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Abstract 


Recent developments in two main phases of cumulus cloud studies are discussed; first 
the relations between the clouds and their energy sources in the subcloud layer and second, 
the interaction between the clouds and their environment in the cloud layer. Under the 
first heading, the various originating impulses for cumulus formation are mentioned, and 
how the character of this impulse affects the spacing, scale and temperature structure of 
the clouds. The possible origins of trade cumuli are considered. The second phase is pri- 
marily concerned with the nature of the resistive forces operating against cloud growth, 
how these arise and how they affect the clouds’ life cycle. The resistive forces are shown 
to be of importance in the transports of momentum, moisture and heat by the cumuli 
and to serve as controls in the feeding of latent energy into the high tropical troposphere 


and hence into the general circulation. 


I. Introduction 


For the past five years or more, the cumulus 
cloud problem has been developing very 
rapidly. A summary of the recent contributions, 
involving an attempt to relate them to each 
other, shows not only that our knowledge of 
convective elements and their detailed structure 
is today far advanced beyond five years ago, 
but that what has been learned about convective 
processes is vitally connected to the most 
important meteorological problems, such as 
the vertical transports of heat, moisture, and 
momentum, the structure of the trade inver- 
sion, and even to the general circulation of 
the atmosphere. 

Many of the key questions and ideas upon 
which the new important developments have 
rested were originally advanced by European 
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or American meteorologists as far back as the 
1930's. The critical observations by which to 
test the hypotheses and from which the orders 
of magnitude of the quantities appearing in 
the equations could be compared, however, 
were largely lacking until the large-scale 
observational programs undertaken by Ameri- 
can agencies in the decade since the war. 
Outstanding among these programs were the 
Woods Hole Oceanographic Institution Ca- 
ribbean Expedition (led by Wyman and 
Woodcock) in 1946, the United States Wea- 
ther Bureau Thunderstorm Project (field 
studies 1946 and 1947), the University of 
Chicago Project Tyrena (Pacific Trade obser- 
vations, 1947), the Hawaiian Pineapple Re- 
search Institute Field Study in 1948, and General 
Electric Company’s Project Cirrus (fight 
programs beginning 1946). The contributions 
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arriving as direct and indirect results of these 
field studies are still accumulating, and form 
the main material for this discussion. Never- 
theless, the time is now approaching when 
many of the newly evolved questions raised 
by analysis of these data must await further 
measurements for their answers. One purpose 
of this summary will be to indicate the direc- 
tions led by several of these questions and some 
possible observational approaches that they 
suggest. 

For the purpose of critical evaluation, 
convection study may be broken down into 
three closely interwoven phases: 


1. The detailed structure of thermals, cumu- 
lus, and thunderclouds; distribution in 
space and time of drafts, gusts, hydro- 
meteors, entrainment, etc. 

. Thermals and convective clouds in rela- 
tion to their energy sources. 

3. Interaction of convective elements with 
their environment; their role in transport 
processes. 


ty 


In the following discussion, results of the 
first phase will be introduced only as they 
contribute to the second and third phases of 
the problem. The physics of the condensation 
process, nuclei, state changes, and cloud 
electrification will be omitted since the 
present summary is almost entirely restricted 
to thermals and smaller cumulus in which 
icing does not occur. For these situations it is 
assumed that only the effects and not the 
details of the latent heat release need be 
considered. 


II. Thermals and Convective Clouds in Rela- 
tion to their Energy Sources 


Because of the dominant role in the energy 
supply of cumulus played by condensation, 
and its regulation by the structure of the 
adjacent air, the major features of individual 
clouds may not differ widely, despite widely 
different initiating impulses. Nevertheless, the 
process which first raises parts of the air to 
saturation may significantly affect several 
features of the cloud layer as a whole, among 
them the spacing and scale of the clouds and 
often the sign and magnitude of the cloud-air 
temperature difference. This originating im- 
pulse may on some occasions be provided by 


organized thermal-convective circulations in 
originally unsaturated air, by unorganized 
thermal-convective bubbles, by occasional 
eddies in a well-mixed layer randomly hitting 
their condensation level, by wave motions 
within a layer or at an interface, or by updrafts 
caused by flow over barriers, convergence in 
the large-scale circulation patterns, sometimes 
due to fronts, seabreezes, or frictional differ- 
ences, and possibly due to other effects not 
yet described. Many times, the effective impulse 
consists of a combination of several of these 
factors. 

Lending itself most easily to theoretical 
analysis is the situation where the impulse for 
cloud formation lies in organized dry convec- 
tion. So far, most of the analytical models 
attempting to describe the initiation of cumulus 
have had as their aim the prediction of the 
vertical motion field in unsaturated air. Two 
extreme physical models have been studied, 
namely, first, the motion field arising from 
uniform heating or cooling of an air layer at 
a horizontal boundary. This is the classical 
“convection cell” problem, first raised experi- 
mentally by BÉNARD (1900a; 1900b; 1901). 
The second model studies the convective 
motions resulting from localized differential 
heating. Despite the different physical situa- 
tions, the underlying premises of all these 
approaches are similar. The basic hypothesis is 
that if convective updrafts are predicted in 
unsaturated, conditionally unstable air, and if 
sufficient moisture is present, cumulus clouds 
may be expected to form at the places of 
maximum lifting. The hydrodynamic equa- 
tions of motion, continuity, state, and energy 
(either the heat conduction equation or the 
first law of thermodynamics) are set up and 
lead, usually by the method of perturbations, 
to a differential equation for the velocity and/or 
temperature field. In such a use of the hydro- 
dynamic equations, an analysis strictly applic- 
able only to laminar motion has been carried 
through and transformed by analogy to apply 
to the mean motions in a turbulent convective 
situation. 


Bénard type convection cells and their relation to 
the cloud problem 


The classical or Bénard-type convection 
cell theories have been well summarized in a 
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review article by STOMMEL (1947a). The 
complexity of these perturbation equations 
is so great that only limited theoretical predic- 
tions are possible. The most famous is the 
criterion for the changeover from molecular 
conduction to laminar cellular convection 
first established by RAYLIEGH (1916). This 
criterion states that when the rate of tempera- 
ture decrease with height (by analogy, poten- 
tial temperature in the atmosphere) exceeds a 
critical value, steady polygonal cells should 
be observed. The critical temperature gradient 
depends directly upon the product of con- 
ductivity and viscosity and inversely upon the 
thermal expansion coefficient and fourth 
power of layer thickness. A preferred ratio 
between horizontal cell dimension and layer 
thickness is also derivable, giving cell diameters 
about 2.5—3.5 times their height. Controlled 
laboratory results (BENARD and AvsEc, 1938; 
SCHMIDT and SAUNDERS, 1938) appear to 
check these criteria in the case of laminar 
canvection. They show steady, honeycomb- 
like, usually hexagonal cells which break 
down into rolls or strips as a horizontal flow 
is superposed. More recent experimental 
studies discussed by SEIDENTOPF (1941) and 
KUETTNER (1949) demonstrate that as tempera- 
ture gradients are further increased to about 
25 times the critical value (or, in practice, 
layer thickness increased by the fourth root of 
this factor, or about 2.2 times), a new change 
takes place in the flow and non-steady, growing 
and dying, convection cells occur. These cells 
have stronger and more localized updrafts, 
diameters approximately equal to their height, 
and a total lifetime given by the ratio of their 
height divided by the mean vertical velocity. 

Intriguing though these results may appear, 
considerable caution must be exercised in 
applying them to the atmosphere and in seeking 
the origin of cumulus clouds in such cells. In 
the first place, uniform heating is assumed in 
the model and very rarely found in nature. 
Second, the models are based on constant 
conductivity and viscosity (laminar motion). 
Their large-scale analogies, the “eddy” viscos- 
ity and conductivity, even if they could be 
measured unambiguously, vary widely with 
space and time, and perhaps by orders of 
magnitude upon the onset of turbulent con- 
vective transport. This raises serious, though 
not necessarily insurmountable, difficulties 


before an observational attempt to determine 
whether the critical condition for the onset of 
Bénard-type convection cells is in real cases 
reached or exceeded. Therefore, in order to 
see if such models have any usefulness in 
application to the real atmosphere, it first must 
be ascertained whether regular polygonal or 
strip-like convective patterns are found. Some 
qualitative evidence has been presented (MAL, 
1930; DURST, 1932; 1933; WALKER, 1933; 
BENARD and AVSEC, 1938; BRUNT, 1938). In 
the case of cumulus formation, Bénard-type 
cells should be sought only in the subcloud 
layer over the open sea, where the heating 
from below can be expected to be the most 
uniform, and where also the lapse rate is 
superadiabatic. Under these conditions, there 
exists some suggestive but not conclusive 
evidence. Among the most recent and pertinent 
is a brilliant observational paper by Woop- 
cock (1940) deducing columnar and strip-like 
convective patterns over the autumnal Atlantic 
from the soaring routines of herring gulls. 
This study was followed by photographic 
investigations of the behavior of smoke trails 
over the sea off Massachusetts, Florida, and 
San Juan, Puerto Rico; the deflections in the 
smoke plumes were consistent with the exi- 
stence of polygonal cells in the lowest hundred 
or so meters above the sea (WYMAN and 
Woopcock, 1947). The hope that similar 
polygonal patterns could be found in trade 
cumulus and that their roots lay in Bénard- 
type convection cells led to the now famous 
Wyman-Woodcock (Woods Hole Oceano- 
graphic Institution) expedition to the San Juan 
area in 1946. Disappointingly enough, not 
only were no polygonal or detectably regular 
cloud patterns found (see Fig. 1), but these 
observations showed that the existence of 
Bénard-type cells throughout any thick por- 
tion of the subcloud layer were almost cer- 
tainly precluded, since the mean lapse rate in 
this area was slightly subadiabatic, becoming 
more stable as the cloud base level was approa- 
ched (see Bunker, HAURWITZ, MALKUS, 
STOMMEL, 1949) thereby shielding the cloud 
layer from such convective motions as might 
exist below. Thus the question of the origin 
of trade cumulus was left open and is still far 
from settled. Due to the great importance of 
the trade cumulus, and their role in the budgets 
of the air over the tropical oceans, this question 
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Fig. 1. Map showing distribution of trade cumuli cal- 
culated from high altitude photographs Nos. 13—26 
made by the Wyman-Woodcock Caribbean expedition 
on the afternoon of April 25, 1946. Scale: 1 cm = 10 
km. Location of center 20°00’ N., 64°25’ W. The wind 
was east-southeast, decreasing slightly with elevation. 
Note the rather irregular spacing of cloudy and clear 
areas. Only the medium and large clouds within the 
53-km circle have been accurately located. The shaded 
areas are distant cloud banks. See Figure 4 for one of 
the original photographs. After J. WyMAN et AL. (1946). 


will be reconsidered in the concluding para- 
graphs of the present section. Meanwhile, 
attention will be turned to a model of cumulus 
origin in organized dry convection which has 
had rather definite observational verification. 


A model of cumulus origin based on differential 
heating; cloud streets 


Following the results of the Wyman-Wood- 
cock Caribbean expedition, in which the 
problem of the origin of trade cumulus 
appeared rather obscure and complicated, the 
attention of the Woods Hole group became 
directed toward the broader question of the 
origin and roots of cumuli in general. It was 
then sought to determine whether any clear- 
cut cases occurred in which convective clouds 
could be related to a specific energy source and 
what role they played in the transports of heat 
and vapor from such a source to a passing air 
stream. The hope was to find a sample situation 
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which could be studied both theoretically and 
observationally, and that the results of such a 
study in which only relatively few factors were 
operating on the airflow might eventually 
provide a starting point for the examination 
of more complex cases. 

A rather striking example of patterned con- 
vective clouds with an identifiable “source” 
are the cloud “streets” formed by islands on 
sunny summer days. These periodically-spaced 
rows of small cumuli (see Figures 2 and 3) are 
frequently seen stretching far downwind of 
small, flat oceanic islands when slightly stable 
air flows across them. 

The basic hypothesis underlying this study 
was that the heating by the island excites one 
or more of the “free periods” of the unsatu- 
rated stable airstream, creating a wave-like 
pattern in the stream lines downstream of the 
“obstacle” and that condensation may some- 
times occur at the wave crests. If the air is also 
conditionally unstable, these may be preferen- 
tial spots for the outbreak of cumuli, and the 
wavelength of the oscillations may then be 
inferred from the cloud spacing. A theoretical 
analysis (MALKUS and STERN, 1951) was under- 
taken, similar in principle to those on convec- 
tion cells, but in which following LANGWELL 
(1951a) and others, the first law of thermo- 
dynamics was the energy equation. Thus a 
distribution of heating as a function of space 
could be assumed, so that the air would be 
heated over the island and not elsewhere. From 
this, it is predicted that the production or non- 
production of lee waves (visible as cloud streets 
if sufficient moisture is present) does not 
depend, within wide limits, upon the character 
or space distribution of the heat source, but 
upon the properties of the air stream before it 


reaches the island. Specifically, no significant | 


downstream oscillations are expected unless 
the atmosphere possesses some non-uniformity 
in the vertical, such as a change in wind speed. 
or a very strong inversion. Second, the 


wavelengths of the lee waves, when they do | 
occur, depend only upon the undisturbed | 


stability, windspeed, and height of the discon- 


tinuity and are independent of the exact form | 


of the heating. For islands about 10 km wide 
under normal insolation, the amplitude of the 
vertical motions in unsaturated air depends 
directly on the heating, inversely on the square 
of the windspeed and may exceed 1/; m/sec. 
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An observational program carried on by 
MALKUS and BUNKER (1952) showed that in 
cases where the island cloud streets were 
observed, a wind change at some upper level 
was found. The spacing checked the predicted 
wavelength in every instance (see Fig. 2 
and 3). In addition to verifying the major 
premises of the simple theory, these observa- 
tions suggested several improvements and new 
lines of 1 inquiry, leading to the removal of the 
assumption of arbitrarily distributed heat 
sources and the inclusion of the effects of 
turbulent eddying and a well-mixed ground 
layer. Even without the refinements, however, 
this investigation has demonstrated that in 
some instances convective clouds do have their 
origin in unsaturated thermal-convective 
motions, and that perturbation methods 
predicting these motions have been successful 
in predicting the occurrence and scale of the 
clouds. 


Other cases of cloud streets have been ob- 
served which may have a similar origin, such 
as those formed by the flatter portions of the 
Hawaiian Islands.1 In the lake country of the 


Fig. 2. Typical cloud street formed over and to the lee 
of SNaatucket Island, Aug. 14, 1950. The low-level wind 
blows from right to left. The first three clouds from the 
right were over the island; those remaining were over 
the downwind waters. The spacing of the clouds was 
measured as 1 km, within experimental error of the 
predicted wavelength of lee waves. After MALKus and 
BUNKER (1952). 


1 Described to the writer by Mr. Wendell A. Mordy 
of the Hawaiian Pineapple Research Institute. 
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Fig. 3. Scale map of cloud street shown in Fig. 2. The 
numbered crosses indicate the positions at which airplane 
soundings were made. After MALKUS and BUNKER (1952). 


far north, especially on the Scandinavian 
peninsula, a common sight is $o—100 km 
long rolls of cumuli extending parallel to the 
wind. These have been described by KuETTNER 
(1949). Since they appear mainly on sunny 
days following cold fronts, it is quite likely 
that their origin also lies in differential heating, 
although wind shear and the earth’s rotation 
must play an important role in their mainte- 
nance. 

Certain common features, however, may be 
sought in all cumuli whose individual origin 
thus lies in any organized dry-convective up- 
drafts and these may provide a way to recognize 
this class of clouds from preliminary observa- 
tions. First, they will usually form in recog- 
nizable patterns more or less simultaneously; 
second, such clouds will possess “roots” or 
detectable updrafts of cloud scale extending 
well below their bases and through a large 
fraction of the subcloud layer; and third, ike 
several other categories of cumuli, they may 
at times be colder than their surroundings 
(discussed by BLEEKER and ANDRE, 1950) “if 
the initial updrafts are sufficiently strong. 


Cumuli formed 


Frequently over heated peninsulas or islands, 
of cumulus clouds appear which can definitely 
be associated with a “‘seabreeze front”. These 
have been studied in Hawaii by LEOPOLD (1949) 
and over Cape Cod, a peninsula on the eastern 


lines or areas of convergence 
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United States seaboard, by Markus, BUNKER, 
and McCastanp (1951). These seabreeze 
clouds owe their existence only indirectly to 
heating, which may not be strong enough to 
release the conditional instability except upon 
the addition of a sharp line of convergence at 
the zone where the seabreeze meets the gradient 
wind or the opposite seabreeze. Whether the 
updrafts caused by the clash of these two wind 
systems serve merely as a trigger to release 
self-sustaining convective clouds with warmer 
temperatures than their surroundings, or 
whether the initial upward velocity is sufficient 
to maintain clouds which consume energy and 
may thereby be colder than their surroundings 
probably varies from instance to instance. In 
any case, the individual clouds would not be 
expected to show deep roots. This expectation 
was supported by the Cape Cod study, in 
which the seabreeze cumuli were noticeably 
warmer than their environment at all levels. 

Another way in which convergence extended 
over wider areas may serve to trigger the 
outbreak of cumuli is by vertical stretching of 
the so-called “mixed” or ground layer. If the 
top of the mixed layer over a wide peninsula, 
for example, lay somewhat below the con- 
densation level and then a seabreeze circulation 
or convergent flow of other origin were 
superposed, so that in places the top of the 
mixed layer approached or reached the con- 
densation level of the air parcels within it, 
random eddying would, under suitable con- 
ditions of instability, cause clusters of cumuli 
to break out in these places. Such spreadout, 
gradual convergence probably plays a signifi- 
cant role in the onset of convective cloudiness 
over the Cape Cod and Florida peninsulas (see 
Byers and RODEBUSH, 1948; Byers and 
BRAHAM, 1949) even when no sharp “‘seabreeze 
front” is present, and convergent flow (perhaps 
of other origin) probably remains as the most 
plausible explanation of the cluster-like forma- 
tions of trade cumuli, as shall be suggested 
presently. Clouds triggered in this manner 
clearly will not show roots in the sense this 
term has been defined in the preceding dis- 
cussion, although warm, moist pockets may 
be detectable in the subcloud layer prior to 
their appearance. 

Probably the most important effect exerted 
by convergence on convective clouds, how- 
ever, lies not just in triggering their outbreak, 
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but in determining to a large extent how far 
they will develop once condensation has 
occurred. The role of convergence in the 
interaction of clouds and environment is 
discussed in the concluding paragraphs of 
this paper. 


On the possible origins of trade cumuli 


The trade cumuli are one of the most 
dramatic and vital features of the subtropical 
atmosphere (see RIEHL, 1950). Day and night 
they grow and decay over a hundred million 
square kilometers of warm oceans, their bases 
at the small stable stratum topping the mixed 
surface layer (— 600 m), their tops pushing 
against, and some of them penetrating, the 
trade inversion (base averages 2 000 m). Their 
degree of development, ranging from scattered 
cumulus humilis to large cumulonimbus, is 
controlled by the height of this inversion, 
which in turn is regulated by the divergence 
field of the passing synoptic systems. In 
relatively undisturbed weather, they have the 
characteristic appearance shown in Fig. 4, 
and form in clusters separated by about 30-km 
wide clear areas, as shown in Fig. 1. As 
these cloud columns rise, eventually mixing 
partially with the inversion air aloft, they 
carry with them vast amounts of energy in the 
form of latent heat, raised initially from eva- 
poration at the sea surface. It was shown b 
Rent, YEH, Markus, and LASEUR (1951) that 
the net heat and vapor transported upward by 
these cumuli is easily sufficient to account for 
the wellknown downstream rise in height of 
the trade inversion, which goes on in the face 
of overall subsidence and divergence. There- 
fore, the relation between these clouds and 
their moisture and energy sources in the sub- 
cloud region warrants study. 

Evidence on these points is still extremely 
scarce, and at first glance, perhaps even contra- 
dictory. However, for one particular locality, 
namely the Caribbean near San Juan, in one 
relatively undisturbed April period, certain 
deductions may be safely drawn. First, these 
individual clouds did not have roots. This is 
not surprising in view of the preceding dis- 
cussion (page 75) which showed that even if 
thermals or convection cells were found at 
low levels above the sea (at ship’s deck level, 
for example), the stabilization of the air as 
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Fig. 4. High altitude photograph No. 15 of trade cumulus made by the Wyman-Woodcock expedition on April 

25, 1946. One of the series of photographs used to construct Fig. 1. The camera was pointed outward from 

the center along line ıs in Fig. I, so that the photograph gives a cross-section exactly in the plane of the east- 

southeast wind. The wind therefore blows from left to right across the photograph and decreases with height. 
After WYMAN et AL. (1946). 


cloud base is approached would damp them 
out long before saturation. Furthermore, 
BUNKER et AL. (1949) showed that the major 
features of the subcloud layer, with maximum 
turbulence at 300 m, decreasing to a minimum 
at cloud base, were incompatible with the 
deep penetration of convective updrafts. 
Supporting evidence was found in horizontal 
accelerometer runs flown below cloud base, 
in which no increased turbulence appeared 
under the clouds. 

A second significant feature of the trade 
cumulus studied by the Wyman-Woodcock 
expedition was their buoyancy. The in-cloud 
temperatures were, on the average, warmer 
than those of the environment, wavering 
between a + 1° C excess and a — 2° C defi- 
ciency near the cloud tops. The buoyancy 
forces were calculated in two cases by STOM- 
MEL (1951) and in each case were positive 
through the lower portion of the cloud, be- 
coming negative in the top few hundred 
meters, indicating overshooting. It seems 
likely that during a major fraction of their 
life cycles, these clouds were self-driving and 
not running on the energy provided by an 
external impulse. 

For these two reasons, it is not reasonable to 


seek the origin of these clouds either in dry ther- 
mals of cloud size or in very sharp convergence 
lines in the flow, but most plausibly in gradual 
variations in the depth of the mixed surface 
layer. On the average, the top of the mixed 
layer was found to be at about 85 % the height 
of the lifting condensation level of the lowest 
air. Variations of several hundred meters from 
day to day and even within a few hours on the 
same day were noted. On some occasions its 
top reached or exceeded the condensation level, 
and it was, in the mean, noticeably higher in 
cloudy than in clear areas. One or more 
mechanisms, then, causing spatial oscillations 
in the thickness of the mixed layer are needed. 
These are probably of the dimensions of the 
cloud clusters, random eddying in these regions 
accounting for the outbreak of individual 
clouds. It is readily showa that fluctuations in 
the thickness of the mixed layer are only very 
weakly coupled with variations in the height 
of the trade inversion and the depth of the 
moist layer as a whole. Attention must there- 
fore be directed to independent, weakly 
convergent flows set up within the subcloud 
layer itself. LANGwELL (1951b) has proposed 
unstable gravitational waves at the interface 
between the nearly adiabatic ground layer and 


78 JOANNE STARR MALKUS 


the isothermal stratum at its top. One should 
also inquire about more random causes such as 
slight inhomogeneities in sea surface roughness, 
radiative properties, or temperatures on a 
10—30 km scale. 

The deductions just made should not be 
generalized without more observations to 
apply to the origin of trade cumuli in all cases. 
In several examples of Pacific clouds, BARRETT 
and RIEHL (1948) have presented evidence 
which may indicicate negative buoyancy 
throughout a cloud. To be sure, their examples 
were cumulonimbus, indicating strong con- 
vergence in the large-scale flow. Under such 
circumstances it is conceivable that an initial 
updraft might provide enough kinetic energy 
to build a cloud against both friction and 
slight stability. The orders of magnitude cannot 
be estimated accurately until more is known 
of the resistive forces on clouds. In conclusion, 
the evidence gathered so far by no means 
precludes the possibility that the origin of trade 
cumuli may be quite different over different 
parts of the subtropical oceans or even over 
the same region at different times. This sugge- 
stion is supported in the next section when 
resistive forces on clouds are discussed. 


III. The Interaction of Convective Elements 
with their Environment 


The preceding section has dealt with the 
relation between the clouds and the subcloud 
layer. Various mechanisms were discussed for 
raising heat and moisture through this layer 
so that they are available for cumulus forma- 
tion. How these properties are further trans- 
ported within and through the cloud layer 
depends upon the dynamics of the individual 
clouds, and particularly upon the interaction 
of the cloud with its environment. 

This approach to the cloud problem as part 
of many larger ones was relatively new with 
the Wyman-Woodcock expedition and was 
phrased clearly in the first written report of 
this group (WyMAN and COLLABORATORS, 
1946). In the classical development of metero- 
logy, cumulus clouds were regarded almost 
entirely as “weather phenomena”, to be 
predicted as part of the synoptic situation. 
The main interest was, therefore, in “stability 
criteria” and the early work on convection 
consisted primarily of the evolution of thermo- 
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dynamic methods to predict the onset and 
probable degree of development of the clouds. 
Their internal structure was regarded as 
interesting mainly as it led to precipitation or 
affected aviation. The interaction of the ascend- 
ing parcel or column with its environment was 
ignored. 

Dissatisfaction with these results soon ap- 
peared, however, when it was frequently found 
that cumuli failed to grow tall even when the 
parcel method stated that a vast amount of 
energy was available to them. Papers by several 
prewar writers pointed to the inhibitory effects 
exerted by the environment, in particular a 
brilliant paper by Namias (1939), who 
described the cutting off of tall cumulus towers 
by mixing. CHRISTIANS (1935) and other 
European authors set up a theoretical frame- 
work for discussion of the several resistance 
forces operating against the growth of cumulus 
clouds, but at that time there were not adequate 
data even to estimate correctly their orders of 
magnitude, and hence their possible importance 
in a cloud’s life cycle. Observations leading 
to the quantitative calculation of one of these 
resistance effects were for the first time 
provided in 1946 by the Wyman-Woodcock 
expedition to the trades. 

This expedition found that temperature 
lapse rates inside clouds were, contrary to 
classical expectations, not moist adiabatic but 
considerably steeper, and that the in-cloud 
temperatures fluctuated with height, as if some 
sort of turbulent mixing process were taking 
place. At this point, STOMMEL (1947b) intro- 
duced, by analogy with the simple fluid jet, 
the idea that environment air was continuously 
“entrained” and mixed with the rising cloud 
air. From comparison of temperatures and 
mixing ratios inside and outside the cloud, 
he devised a method to determine the fractional 
amount of the entrainment, level by level. In 
the tradewind clouds studied, it frequently 
corresponded to a doubling of the mass flux 
of the rising air in 1 km ascent. 

While the purpose of the entrainment paper 
by Stommel was specifically to explain the 
departure from a moist adiabatic lapse rate 
inside trade cumulus, it provided the corner- 
stone for many contributions to follow. From 
consideration of the effect exerted by the 
surrounding air upon a cloud and its logical 
converse, the effect exerted by the cloud upon 


RECENT ADVANCES IN THE STUDY OF CONVECTIVE CLOUDS 79 


the surrounding air, a framework can be set 
up to investigate many previously baffling 
problems. For example, it can now be under- 
stood why observations sometimes showed 
cloud tops which ceased growing while 
warmer than the environment, and why 
“pulsating” cumulus towers are rare. A sug- 
gestion can be made concerning the different 
features observed in trade cumuli contrasted 
with their middle-latitude counterparts. The 
effects of windshear upon convection can be 
discussed analytically, leading to at least 
preliminary models of convective clouds as 
transport elements. In effect, the role of cumuli, 
especially of trade cumuli, as an integral part 
of the layer in which they grow, rather than 
as incidental photogenic objects sometimes 
acting up to produce aircraft hazards and rain- 
all, is now a fruitful topic of investigation. 
Clearly, the question of predicting the onset 
and growth of these clouds from initial air 
conditions becomes considerably more com- 
plex as a result of entrainment and other 
resistive effects. 


Stability criteria and their modifications by resis- 
tive forces 


The model of the air parcel which rises 
without interaction with its environment, 
while now regarded as a drastically over- 
simplified one to describe real cases of convec- 
tion, still serves an important function as a 
clearcut limiting case, just as the undamped 
harmonic oscillator is a clearcut limiting case 
in studying the behavior of real, frictionally 
damped oscillators. The original purpose in 
setting up stability criteria was to establish the 
conditions of thermal stratification in the 
environment under which an air parcel when 
displaced slightly from its equilibrium position 
would acquire an initial acceleration away from 
its origin. Since a rising, saturated, non-inter- 
acting air parcel would cool at the moist adia- 
batic lapse rate, the well-known result was 
obtained that if the environment cooled more 
rapidly with height than this rate, the ascending 
parcel would become less dense than its sur- 
roundings and would accelerate upward. 

Strictly speaking, these criteria regarding 
initial accelerations stand unmodified, since 
frictional and resistive forces come into action 
only after an initial velocity has been acquired. 


The difficulty arises in attempting to apply 
them to finite convective motions. The first 
quantitative recognition that buoyancy forces 
were cut down by the parcel’s environment 
came with the slice method of J. BJERKNES 
(1938). He recognized that the compensatory 
sinking (usually dry adiabatic) of the outside 
air would cause its warming, thus cutting 
down the relative temperature excess of the 
rising parcel and thereby its buoyancy. This 
effect was shown to be greater the more a 
region was filled with updrafts because of the 
requirement thereby of more limited, hence 
stronger downdrafts and consequent intensi- 
fied warming of the sinking air. 

The entrainment of drier environment air, 
however, ordinarily reduces buoyancy forces 
far more than does compensatory sinking. 
Simple calculations show that even under 
exaggerated conditions of 3/10 of the area 
occupied by active updrafts, the effects of 
compensatory sinking reduce the buoyant 
energy predicted by the parcel method only 
by about 40 %. Upon entrainment of a nor- 
mally observed fraction of outside air with 
relative humidity about 70 %, the buoyant 
energy may be reduced by over 90 % from 
that predicted by parcel considerations. 

AUSTIN (1948) has discussed in detail the 
reduction in buoyant energy under various 
assumed rates of entrainment and environment 
humidities. He showed what the lapse rates 
inside clouds would become, and hence what 
environment lapse rates would be necessary 
for cloud growth, thus illuminating the 
inhibitory effects of a dry environment upon 
the development of large cumuli. Such 
methods, however, can only be expected to 
give statistical probabilities for cloud growth, 
rather than clearcut predictions from an initial 
stratification. The reason is that one or more 
assumptions, usually including that of an 
average entrainment rate, is required. Actually, 
the entrainment rate which will prevail at a 
given time, under given conditions is itself 
a complex, adjusting function of the vertical 
motions, cloud size and other factors, and may 
be highly variable. Hence it is probably no 
longer fruitful to discuss convection from the 
point of view of onset criteria, parcels, and 
initial perturbations, except for reference, as 
is the undamped harmonic oscillator in physics. 
The more difficult problem of a cloud’s “life 
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cycle” and the interacting adjustments between 
buoyancy and resistance after the updrafts are 
of finite size must sooner or later be tackled. 


Present knowledge of the resistive effects on cumu- 
lus clouds 


It has been shown that the buoyancy forces 
producing cloud growth are reduced by two 
causes: first, by the compensatory sinking of 
the environment; and second and more drastic- 
ally by the entrainment of drier outside air. 
Clearly, if the buoyancy forces on its air are 
upward, the cloud grows, and hence for 
stability criteria, the net buoyancy is all that 
need be determined. Nevertheless, these two 
effects, compensatory sinking and entrainment, 
further reduce the kinetic energy in an existing 
updraft via ordinary friction or reduction of 
the draft’s upward momentum. Therefore, 
these too must be considered in a “life cycle” 
approach, since anything reducing the vertical 
motions in turn acts back to reduce the resistive 
effects, which: control the buoyant energy 
release, and so forth. The nature of these 
“drag” forces has been set forth in an excellent 
treatise by SCHMIDT (1947) who shows that 
they are exactly the same forces which would 
oppose the motion of a non-buoyant jet of 
water created within a resting fluid. The first 
he calls “form” or “profile drag”, and the 
second “skin friction”, similar to terminology 
in aerodynamics. The form drag is due to the 
pushing aside of the surrounding medium as 
the jet grows and accelerates it back as a 
“countercurrent”. The skin friction arises 
from viscous forces at the moving boundary, 
which in the fluid jet give rise to entrainment. 
In contrast to aerodynamics, however, in 
which solid bodies are moved through fluids, 
the form or profile drag may prove to be 
neglectable in many aspects of the cloud 
problem. Clearly, it will be most important 
as the cloud tower grows rapidly upward and 
the countercurrent is being established, and 
less so when the cloud later approaches a 
steady state. Entrainment and its consequent 
reduction of upward velocity, however, is 
clearly operative throughout the entire life of 
the updraft. This conclusion concerning relative 
drag magnitudes may be invalid, though, in 
the later considerations of the horizontal 
motion of the updraft through the air, since 


a horizontal countercurrent must then con- 
tinually be produced similar to the case of a 
moving solid, thus creating a pressure difference 
in the horizontal across the draft. 

For preliminary models of a cloud’s growth 
and life cycle, it would indeed be fortunate 
were both resistive effects due to counter- 
currents neglectable compared to those due to 
entrainment. In view. of the above discussion 
it seems plausible to proceed upon this tentative 
assumption, especially since both buoyancy 
reduction and momentum reduction by en- 
trainment can be readily dealt with analytically. 
The buoyancy reduction by countercurrents 
was shown to be relatively smaller, and their 
frictional effects seem prohibitively complex 
at present. Intensive observational tests of such 
models will then be called for, to determine 
whether this assumption requires later modi- 
fication. 

Preparatory to any analytic treatment of the 
interaction of entrainment and buoyancy in 
cloud dynamics, however, a further discussion, 
of the entrainment mechanism is required. Two 
points of view prevail at present concerning 
this topic. It is this writer's belief that these 
points of view are by no means incompatible, 
but rather arise from two different operational 
definitions of entrainment. STOMMEL (1947b) 
calculated the proportion of environment air 
of known properties which must have mixed 
with the rising air to produce the observed in- 
cloud distribution of temperature and moisture. 
By his analogy with the spreading jet in fluid 
mechanics he implicitly assumed that all this 
mixed air was dragged along with the updraft, 
permitting its cross section to expand with 
height. Other writers (AUSTIN, 1948; HoucH- 
TON and CRAMER, 1951) have assumed ex- 
plicitly that the draft cross section remains 
nearly constant with height, in more apparent 
agreement with existing sparse observations. 
The latter writers contend that the major 
entrainment is that which is required to 
preserve continuity in a constant cross-sectional 
cloud whose ascent rate increases upward. 
The observed entrainment rates called upon 
by them are based upon direct measurement 
of inflow winds around clouds, such as the 
calculations described by ByErs and BRAHAM 
(1949) and are much smaller than those found 
by Stommel. Clearly, there is no conflict 
between these two viewpoints if the possibility 
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is admitted that all the air which is in or 
becomes mixed with that ascending is not 
retained indefinitely within the updraft, but 
that some is “detrained” during the ascent. 

Stommel’s calculations might be said to give 
the “gross entrainment” or the total amount 
of environment air mixed with cloud air, 
while the Austin, Houghton-Cramer method 
considers “net” or as they call it “dynamic” 
entrainment, the difference between the two 
amounts in a given cloud being detrained. An 
observational test of this hypothesis would be 
provided if a vertical velocity profile of the 
cloud were measured simultaneously with in- 
cloud and environment soundings. The amount 
of net influx required to preserve continuity 
could then be calculated and compared to the 
total amount of outside air mixed with cloud 
air. Although such measurements are still 
wanting, some less direct evidence supporting 
this hypothesis appears in the following 
sections. 


Recent approaches to cloud dynamics 


HOUGHTON and CRAMER (1951) have con- 
structed the first completely analytic model of 
the interaction of buoyancy, vertical motions, 
and entrainment in a cloud. Although their 
work, because of the assumptions necessitated 
by present data, cannot yet be regarded as a 
“life cycle” model, it nevertheless represents a 
significant step in this direction. 

Under the basic assumption that dynamic 
entrainment be the only resistive effect, these 
authors obtained and solved differential equa- 
tions for the height distribution of vertical 
motions, net mass influx, and liquid water. 
These equations were derived from the laws 
of mass, momentum and energy conservation 
under the additional restriction to a steady- 
state cloud of constant cross section. Their 
work has the great advantage over previous 
attempts in that the entrainment rate need not 
be assumed or inferred but appears as one of 
the dependent variables in the system of 
equations and hence can be predicted. The 
vertical velocities so calculated are of the 
observed order of magnitude, but probably 
somewhat too great. Their derived influx rates 
are comparable to but larger than those found 
by the Thunderstorm Project (Bvers and 
BRAHAM, 1949) and smaller than most of those 


found by Stommer (1947b). They found that 
although the liquid water content of their 
clouds was considerably less than that predicted 
by the parcel method, dynamic entrainment 
of even very dry outside air does not suffice 
totally to “dry out” the cloud.’ This appears 
to be in conflict with many observations 
(Markus, 1952). Clearly, then, the greatest 
departure from realism of such a model is that 
it fails to provide a way in which updrafts 
may die out with height, other than by in- 
creasing stabilization of the environment. 
Dynamic entrainment ceases to operate as a 
resistance when the ascent rate has stopped 
increasing upward, since no more influx is 
then required to maintain continuity. The 
authors, however, themselves suggested that 
additional resistances, such as form or profile 
drag, weight of liquid water, or further mixing 
with the environment probably should be 
introduced. 

A second entrainment paper by STOMMEL 
(1951) appearing simultaneously with that of 
Houghton and Cramer provides interesting 
comparison with their work. He computed, 
again using the Wyman-Woodcock soundings, 
the vertical velocity distribution within these 
trade cumuli. Similar to Houghton and 
Cramer, he assumed a steady state and that the 
only resistive effect was the entrainment in- 
ferred from observed soundings. The vertical 
velocities obtained were comparable to those of 
Houghton and Cramer, although the environ- 
ment air of the tradewind clouds was far 
moister. Some more striking differences also 
appeared. Although Stommel’s rate of increase 
of updraft with height was nearly the same as 
that predicted by Houghton and Cramer, the 
measured “gross entrainment’ was far greater 
than that needed to satisfy continuity and 
continued at a high rate above that level where 
the vertical motions began to diminish upward. 

If further observational work corroborates 
these indications, a mechanism to cause such 
additional entrainment or mixing must be 
sought. A factor not included in the Houghton- 
Cramer model which operates in real situations 
is shear in the external wind. Observations of 
clouds in the Woods Hole (eastern U. S.) area 
seem to indicate that those clouds which cease 
to grow despite warmer temperatures at their 
tops occur on days of pronounced wind shear. 
The Thunderstorm Project (Byers and Bra- 
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HAM, 1949) also found that strong shear 
exerted an inhibitory effect upon thunderstorm 
development. 


Some effects of windshear upon convection 


The crux of these investigations is that 
momentum is a vector quantity. While the 
preceding discussion has shown that the vertical 
momentum relations in a cumulus cloud are 
extremely complex, due to the interaction of 
buoyancy and friction, the horizontal compo- 
nent is altered by friction alone. If a parcel or 
updraft (the presence of liquid water is, for the 
moment, unimportant) were to rise without 
interaction through an environment in which 
the external wind changes with height, it would 
acquire a velocity relative to the air, since its 
forward momentum would remain that of its 
original level. If the shear vector between any 
two levels is defined as the upper minus the 
lower wind, the parcel’s relative velocity 
would be in the sense opposite to this. 

It has been shown by Markus (1949a; 1952) 
that frictional forces act to reduce but do not 
eliminate this relative velocity. A striking 
verification of this prediction was found near 
Woods Hole, when on several days of marked 
wind turning with height, small cumulus were 
actually observed moving at an angle of 45° 
to the wind at their midlevel. This case has 
been discussed quantitatively by Markus, 
Bunker and MCCASLAND (1949). 

The entrainment of outside air is one factor 
operating to alter the momentum of the up- 
draft toward that of its surroundings. The total 
amount of outside air mixed with the updraft 
must be included, hence the gross entrainment 
rates found by STOMMEL (1947b) were used in 
these considerations. When this is assumed to 
be the only resistive force, relative horizontal 
velocities of about 2 m/sec should be common 
about 1 km above draft base, for normal wind 
shears. If form or pressure drag creates addi- 
tional resistance of equal magnitude, the relative 
velocities are decreased about 40 %. Direct 
verification of relative cloud-air velocities was 
found in the Thunderstorm Projects tracking 
of clouds by radar (see Byers and BRAHAM, 
1949). At roughly 1 km above cloud base, they 
found average relative speeds of 2 m/sec, 
increasing to an average of 5 m/sec at 6 km 
(fairly large average windshears). 


The relative velocity between draft and. 
environment has been shown by Markus 
(1949a) to produce important asymmetries 
between the location of updrafts and that of 
the visible, liquid cloud. How these may arise 
is easy to see. If the horizontal speed of the 
draft through the air is of the order of 1 m/sec, 
it will be greater than the normal rate of 
inflow of entrained air, which for a stationary 
draft might be fairly symmetrical around the 
periphery. Hence by vector addition, the 
inflow rate on the relative upwind side is 
greatly increased and actual outflow or de- 
trainment of cloudy air from the updraft may 
take place on the opposite side. If verifiable 
observationally, a mechanism has clearly been 
presented whereby the total influx or “gross 
entrainment” is caused to be greater than the 
dynamic or “net entrainment” required to 
preserve continuity. 

Several types of observation verify this 
prediction. The most striking are time-lapse 
motion pictures (motion pictures exposed at 
long intervals, say one every five seconds and 
projected at normal speed of 16 per second). 
A number of sequences made in the Woods 
Hole area showed cloud fragments streaming 
out from one side of the moving cumul, 
becoming caught in the adjacent downdraft 
and dissipating, while the cloud grew vigo- 
rously and put out new towers on the opposite 
side. In each case, the cloud dissipated on the 
downwind side of the shear vector and grew 
on the upwind side (see discussion by Markus, 
1949b; 1952). Nineteen sequences of films 
made by the Thunderstorm Project in Florida 
were also checked and the prediction verified 
in every case that the windshear for the region 
could be determined reliably (17 cases). In 
Fig. 4, showing trade cumulus, the wind 
blows from left to right across the photograph 
and decreases upward so that the shear vector 
points toward the left. Signs of active growth, 
new turrets, and turbulence are apparent on 
the right edges of the clouds, while they 
appear fuzzy and weak on the left edges. — 

An additional verification of these asym- 
metries in trade cumulus was found from the 
accelerometer records made in the horizontal 
airplane traverses by the Wyman-Woodcock 
expedition. Fig. 5 presents an average of 
eight such horizontal traverses through trade 
cumuli, in each of which the ambient wind 
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increased upward. If it can be assumed that 
active updrafts are associated with cloud-scale 
turbulence the data are easily interpretable as 
follows: In the region 6, where cloud air has 
been blown off the updraft, the water droplets 
may be expected to evaporate and the sou 
die. In the region D, where air is being brought 
into the updraft, new cloud may be expected 
to torm. Similar observations were reproduced 
by Markus, BUNKER and MCCASLAND (1949) 
in the Woods Hole area. Thus it is seen that a 
cumulus is not simply an aggregate of droplets, 
drifting passively with the wind, but rather 
is a transient, dynamic balance between rapid 
growth and equally rapid destruction. The 
air is seen to have a component of motion 
through the updraft. An extreme case of air 
motion through a draft is found in clouds 
which appear to remain standing near moun- 
tain tops for hours despite strong winds. 

The asymmetrical distribution of upward 
accelerations within thunderclouds and its 
relation to windshear has been discusséd ana- 
lytically by LANGweLt (1951c). She showed 
that a large cumulus is most likely to begin 
transition to a mature stage on its edge toward 
which the shear vector points (6 region), 
because there the updrafts are weakest or 
absent and yet liquid particles plentiful enough 
to cause a significant downward drag. Her 
predictions concerning the side of convective 
clouds that first produces rain were well 
verified from data taken by the Thunderstorm 
Project. The downward transport of momen- 
tum by the downdraft obeys principles similar 
to those just presented and is discussed by 
Byers and BRAHAM (1949). The role of this 
downward transport of usually very high 
westerly momentum has been shown to be 
important in the cold air outflow from 
thunderstorms, the initiation of new thun- 
derstorm cells, and the formation and move- 
ment of squall lines. 


The slopes of cumuli; trade cumuli vs. those of 


middle latitudes 


The slope of an updraft at any level depends 
only upon the ratio of the vertical speed to 
the departure of the horizontal speed at that 
level from that at its base. Therefore, it is very 
easy to calculate the slopes of the clouds 
discussed in Stommel’s second entrainment 
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Fig. 5. Average of eight horizontal airplane traverses 
with accelerometer through middle levels of trade cumuli. 
Traverses were made by the Wyman-Woodcock 
Caribbean expedition, April 1946. In each of these cases, 
the wind increased with height about 1 mps/km. The 
cloud-scale turbulence was at least an order of magnitude 
in excess of clear-scale turbulence. The time of entrance 
and exit into visible cloud was recorded directly on the 
accelerometer record. The arrow labelled U, is the cloud 
motion seen by an observer moving with the wind; 
an observer riding on the cloud would feel a wind in 
the opposite direction. After MALKUS (1949a). 


paper (1951). The horizontal speed of the 
clouds can be calculated from the entrainment 
rate and observed wind shear using the tech- 
nique evolved by Markus (1949a). Stommel 
already calculated the vertical speeds. The 
result of this calculation is reproduced in 
Fig. 6. This should be compared with the 
photograph presented earlier, in Fig. 4, 
which was made in the same area at the same 
time as the sounding used by Stommel. Due 
to negative shear (easterlies decreasing upward), 
the predicted slope of the cloud is upwind, as 
is the observed slope. Most of the clouds in 
the photograph appear to slope more strongly 
backward than does the theoretical cloud. 
This cannot be due to the 18 000 ft altitude 
from which the photograph was taken, since 
even the nearest clouds were about 70 000 ft 
from the camera. It may be due to additional 
resistive forces such as form drag, not included 
in the theory, or it may be that most of the 
clouds pictured here were in a later stage of 
their life cycles than Stommel’s, thus having 
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Fig. 6. Theoretical slope of the trade cumulus on April 
25, 1946 whose vertical velocity profile was calculated 
by STOMMEL (1951). The dark arrows to the left are the 
observed external wind. The light curved arrows repre- 
sent schematically the predicted air motions into and 
out of the cloud, according to Markus (1949b). The 
cloud should be building on the right and decaying on 
the left side as it moves from left to right more rapidly 
than the external wind. Compare with Fig. 4. 


slower ascent speeds. A more detailed study 
of cloud slopes (Markus, 1952) shows that 
greater horizontal friction forces or weaker 
updrafts produce more slant, as does stronger 
external shear, other conditions remaining 
the same. 

The inhibitory effect of a strong slant upon 
cloud growth has been discussed by Byers and 
BRAHAM (1949). They demonstrated that a 
cloud grows by the production of successive 
turrets at or near its top. If wind shear is weak 
the new turrets grow in air already moistened 
by dissipated predecessors. If external shear is 
strong, the old turrets become strongly slanting 
as their updrafts decay. The new turrets with 
strong drafts grow nearly vertically and hence 
are penetrating drier air. That small cumuli 
as well as thunderclouds grow in this manner 
has been verified by time-lapse movies in the 
Woods Hole area where individual towers have 


been observed to have life spans lasting from 
six to sixteen minutes. 

An additional inhibitory effect of wind shear 
upon cloud growth has been suggested in the 
foregoing discussions. The relative cloud-air 
horizontal velocity increases linearly with 
wind shear. If it can be shown, as proposed 
here, that additional entrainment results from 
this motion of the cloud through the air, then 
the resistive effects and buoyancy reduction 
are substantially augmented by wind shear. A 
suggestion concerning the difference between 
trade cumulus and their extra-tropical counter- 
parts may now be feasible. 

The constancy of the trade wind through the 
subtropical cloud layer is well known. The 
average shear vector through a corresponding 
layer in the Woods Hole area, for example, 
was 2—3 times greater (6 mps/km) in the 
summer months than near San Juan during the 
Wyman-Woodcock expedition (2—3 mps/ 
km). From the data of RieL et AL. (1951), the 
shear in the Pacific Trade appears to be even 
smaller (~ 1 mps/km). Trade cumuli may 
therefore be operating against considerably 
smaller resistive forces than middle latitude 
cumuli. This would mean that the clouds are 
more easily triggered by weaker impulses (see 
concluding paragraphs of Section II) and that 
smaller buoyancy forces are required to pro- 
duce comparable vertical speeds. Turbulence 
within trade cumuli could also be expected to 
be less than that in middle latitude clouds for 
the same ascent rates. This would happen 
because the in-cloud lapse rates need not be 
so much steeper than moist adiabatic. AUSTIN 
(1948) has pointed out that turbulence within 
a cloud should be greater the greater the 
departure of the in-cloud lapse rate from moist 
adiabatic, because an unstable state is thereby 
presented for the small-scale, moist adiabatic 
overturning near cloud center. 

It was mentioned that the observed clouds 
in Fig. 4 may have slanted more than predic- 
ted by the Stommel-Malkus theory because of 
the omission of horizontal form drag from the 
model. Observations from the Thunderstorm 
Project also indicate that form drag may be 
important in reducing the relative horizontal 
cloud-air velocity. Some entirely different 
evidence presented by RIEHL et AL. (1951) re- 
garding the transports in the subtropical cloud 
layer further support this possibility. 
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Transports by convective elements 


Most of the work on cumulus clouds dis- 
cussed so far emphasizes the effects of the 
environment upon the cloud. The data 
gathered by the University of Chicago Project 
Tyrena (reported on by RıeHr, YEH, MALKus 
and La SEUR, 1951) permitted for the first time 
framing the problem in reverse. Clearly, if the 
environment exerts a frictional force on the 
cloud, the cloud exerts an equal and opposite 
force on the environment. Also, according to 
the picture drawn above, the environment not 
only supplies drier air to the cloud, but by 
“detrainment” the cloud supplies moister air 
to the environment. 

From the Tyrena observations it was possible 
to calculate effective large-scale momentum and 
moisture transport coefficients layer by layer 
for a section of the Pacific Trade. In the cloud 
layer the momentum “Austausch” was found 
to be significantly larger than the moisture 
“Austausch”. Since the clouds themselves are 
the sole contributors to these coefficients within 
this layer (BUNKER et AL. 1949, have shown 
the turbulence in between clouds is at least 
an order of magnitude smaller than cloud 
scale) these data suggest that for equal gradients, 
trade cumulus clouds are far more effective 
transporters of momentum than of water 
vapor. If form drag is of roughly equal im- 
portance as entrainment in the horizontal 
momentum relations of the clouds, this 
observation is explainable qualitatively. If by 
means of pushing the air around itself, the 
cloud thereby exerts a force on the environ- 
ment, it imparts some of its horizontal momen- 
tum to the surrounding air. It also exchanges 
horizontal momentum with the surrounding 
air by means of mixing or entrainment. It 
exchanges water vapor with the environment 
only by the latter process. Hence one might 
expect intuitively that the measured “eddy” 
viscosity for a convective layer might be 
considerably in excess of the “eddy” diffusivity. 
It can be shown analytically, even by a rather 
crude model, that if form drag and entrainment 
contribute friction forces (horizontal) of 
roughly equal magnitude on the elements of 
a convective layer, that this ratio comes out 
considerably greater than one. A required 
assumption is that the individual elements 
somehow overshoot their equilibrium positions 
and sink partially back down. 


STOMMEL (1951) has demonstrated the over- 
shooting of trade cumuli in the Caribbean, 
and RIEHL et AL. (1951) have further discussed 
how penetration of even a small fraction of the 
clouds into the trade inversion and later dissipa- 
tion of their tops by mixing can provide 
sufficient moisture to account for the observed 
downstream rise in height and weakening of 
the inversion. 

These latter calculations emphasize the 
tremendous amounts of energy continually 
pumped into the troposphere by the trade 
cumuli. The magnitude may perhaps be 
grasped by the following: If a one square kilo- 
meter area is so % covered by trade cumulus, 
only a small fraction of which need be actively 
rising at any one time, they can easily provide 
a net upward transport of latent energy in 
water vapor equal to that released by a 500 
kg bomb of trinitrotoluene (3 000 cal/gm) 
exploding every minute. If this is summed 
up over all the tropical oceanic area normally 
covered by trade cumuli, a net upward trans- 
port of latent energy is given which is greater 
by a factor of roughly four or five, than the 
rate of kinetic energy dissipation by all the 
wind systems of the globe. 

It has been strikingly pointed out by RıEHL 
(1950) however, that tropical convection does 
not uniformly feed energy into the atmospheric 
heat engine, but mainly in the easterly waves 
and polar troughs which, in passing, superpose 
strong convergence upon the flow. This 
large-scale convergence results in a greatly 
deepened “moist layer” and a tremendous 
diminution in the resistive forces operating 
against the growth of the cumuli. In the first 
place, the stable lid usually provided by the 
trade inversion is raised or removed. Second, 
net ascending motion in the region diminishes 
the required compensatory sinking and hence 
the reduction in buoyancy predicted by the 
slice method (CRESSMAN, 1946, has discussed 
this effect in detail) and finally, the growing 
clouds are surrounded by a deep layer of 
moist environment up to great heights so that 
entrainment is less effective in reducing their 
buoyancy. 

How this energy, once made available in 
these disturbed regions to the upper subtropical 
atmosphere, is then injected into middle 
latitudes is being studied by Richl and other 
investigators of the role of the tropics in the 
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general circulation. It is now apparent, how- 
ever, that cumulus clouds are important 
driving elements, and the forces restraining 
them are important controls, in the basic 
machinery of the atmospheric heat engine. 


IV. Conclusion 


In concluding, it can be seen that many more 
questions concerning cumulus cloud dynamics 
have been raised by this paper than have been 
answered by it. Many of these are framed in 
terms of specific sets of measurements which 
must be carried out. Consequently, the 
Woods Hole Oceanographic Institution has 
planned a new program of meteorological 
observations, to commence in the Spring of 
1952. The principal observing tool is to be a 
large, slow-flying aircraft equipped with 
psychrograph, accelerometer, variometer, 
liquid water collector, drift sight, and cameras. 


By numerous flights through and near both 
trade and middle latitude cumuli, it is hoped 
to answer some of the existing questions and 
upon these answers to base new inquiries which 
are better evolved and more searching. 
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Abstract 


It is suggested that the jet stream at about 12 km-level may be idealized as a system of two 
or three zonal currents with uniform but different values of absolute vertical vorticity £, and 
such that the wind velocity itself is continuous at each boundary. 

The approximation of a horizontal non-divergent motion is made for the jet-stream’s waves 
(chapter II), and a proper system of equations is derived for these waves (by using the meridional 
coordinate y defined by tanh y = sing (g = latitude), the equations are particularly simple 
even when the earth’s curvature is taken into account, and in the case of constant £ the problem 
is reduced to the resolution of the very classical equation V *y = o, where y is the local disturb- 
ance of the stream function. 

The case of a simple jet (two currents), then that of a double jet (three currents) are succes- 
sively considered in the next chapter, the wave motion being assumed to vanish at each pole, 
and it is also shown how some of the results can be generalized for a larger number of currents. 
(In any case the dispersion equation, giving the angular phase-velocity & in terms of the angular 
wave-length À (inverse of the wave number, which must be an integer), is an algebraic one 
in &, and its degree equals the number of internal boundaries.) The waves of a simple jet are 
therefore all stable (their dispersion equation is «= (9 — 094, if Wo is the angular wind- 
velocity at the axis, and 20, denotes the ¢-discontinuity, positive northwards). On the contrary, 
a system of more than two currents may have unstable (amplified) waves, if one discontinuity 
20, is negative while the other ones are all positive. If in addition ©, is small, the wave lengths 
of the unstable waves form a definite series of spectral lines, the maximum number of which 
equals the number of the positive discontinuities, and the phase velocity of these waves is ap- 
proximately the wind velocity at the latitude of the negative discontinuity: in the special case of 
a double jet, there is therefore at most one unstable wave lengthA,, given by Ay = (@1 —@) /oy, 
where the subscripts o and 1 are relative to the negative and positive discontinuity, respecti- 
vely. For all the waves of a simple or double jet, there is a maximum of the amplitude at 
a definite latitude and the wave motion becomes negligible at a distance of the order of the 
wave length. 

The effect of horizontal eddy-viscosity is discussed in chapter IV. If the kinematical coefficient 
is of the order of 10° C. G. S. at most, it is found that the wave motion is generally modified 
only in the vicinity of the boundaries. The conclusions are however different in the special 
case where one discontinuity is very small: some of the stable waves are then changed to damped 
waves, while the unstable waves are practically unmodified. The problem involves Bessel 
functions of order 1/3. 

Chapter V contains a tentative theory of the semi-permanent centers of action, based on the 
properties of the stationary waves of a simple or double jet. For any given wave number, 
stationary waves are found to be possible in a westerly simple jet with a maximum wind- 
velocity Uy along the axis, provided Uy = o9L/27 (resonance condition), where L is the wave- 
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length at the latitude of the axis, and if oy is assumed to have the constant value 3.5 x 1075 
C.G.S., the results do not differ essentially from those derived from Rossby’s theory. But 
this agreement disappears in the case of a double jet, since there are generally at least two very 
different such jets corresponding to the same given wave number, and that may explain why 
the correlation between the zonal index and the number of centers of action is often rather 
loose. It is also possible to have two simultaneous different stationary wave lengths in the same 
jet-stream. 

The last chapter is similarly devoted to a tentative theory of cyclone waves, based on the 
properties of the unstable waves appearing in a simple or double jet, when such a jet is slightly 
deformed by an additional negative ¢-discontinuity. It is suggested that the agent of this defor- 
mation may be the vertical convection created in the lower atmosphere by thermal instability, 
since this convection tends to reduce the wind velocity in a relatively narrow zone above the 
instability area. If the assumption is made that the unstable waves can represent cyclone waves 
provided their wave length is larger than 1,000 km and smaller than 4,000 km, with an ampli- 
fication coefficient larger than 10° C. G. S., it is found that cyclone waves can develop in a 
simple jet if the deformation takes place in either one of two “‘sensitive zones’’ located at less 
than 300 km on each side of the axis, and not larger than soo km about. For a double jet there 
are normally four such sensitive zones. All the properties of these theoretical cyclone waves 
agree quite well with that of the actual cyclone waves of the polar fronts, except for the phase 
velocity which is somewhat too large, but this discrepancy is not surprising since the dynamical 
effect of the lower atmosphere has been disregarded. Finally it is suggested that the same theory 
can also apply to tropical cyclones and easterly waves, if the assumption is made that these 
features are produced when one of the sensitive zones of the jet stream happens to penetrate 
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into a monsoon area for instance. 


I. Introduction 


Au cours des 30 dernières années, deux 
tendances principales et apparemment diver- 
gentes se sont manifestées en météorologie 
théorique. La première, directement dérivée 
des idées norvégiennes, consiste à rattacher les 
propriétés des perturbations atmosphériques à 
la baroclinité de la troposphère moyenne 
(propriété traduisant le fait que les surfaces 
isothermes sont inclinées par rapport aux 
surfaces isobares) et au gradient vertical du 
vent qui lui est nécessairement associé, et à 
faire jouer en outre un rôle essentiel aux 
surfaces de discontinuité de la température. 
Les problèmes de perturbation traduisant cette 
conception se sont révélés particulièrement 
complexes, car les modèles hydrodynamiques 
auxquels on est conduit sont essentiellement à 
trois dimensions. La seconde tendance, déve- 
loppée surtout après les travaux de RossBy 
(1939), attribue au contraire le rôle principal 
à la couche barotrope voisine de la tropopause, 
et qui est celle où la circulation générale possède 
la plus grande énergie aux latitudes moyennes; 
en outre, comme les mouvements dans cette 
couche sont tous très approximativement 
horizontaux et en moyenne zonaux 4 l'échelle 
synoptique, on est ramené à des problèmes 
hydrodynamiques à deux dimensions seulement, 
donc relativement simples. Nous citerons en 
particulier, comme travaux se rattachant à 


cette tendance, ceux d’Haurwitz (1940), de 
Kuo (1949 et 1950), et ceux de PLATZMAN 
(1949). Le dernier auteur est le premier qui 
semble avoir montré l’importance des discon- 
tinuités du tourbillon dans la couche barotrope: 
son modèle de circulation générale comporte, 
en effet, deux courants zonaux où la compo- 
sante verticale du tourbillon absolu est uni- 
forme, avec une discontinuité de ce tourbillon 
sur leur frontière commune, ce modèle étant 
manifestement inspiré de celui considéré 
antérieurement par RossBy (1947). 

L'objet du présent travail est précisément 
l’étude des ondes horizontales associées aux 
discontinuités du tourbillon dans la couche 
barotrope, en vue d’applications aux ondes 
atmosphériques réelles, et en particulier aux 
ondes de cyclone; nous retrouverons l’un des 
résultats de Platzman, mais au lieu de considérer 
seulement le cas d’une discontinuité unique 
(jet-stream simple), qui ne conduit qu’à des 
ondes stables, nous envisagerons celui où l’on 
a deux discontinuités à des latitudes différentes 
(jet-stream double), de sorte que nous arrive- 
rons à des résultats plus généraux et plus 
directement applicables (d’après CRESSMAN 
1950), le cas le plus fréquent est en effet celui 
d’un jet double). 

Bien que le développement de l’aérologie 
semble confirmer de plus en plus que c’est bien 
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la couche barotrope voisine de l'altitude 12 km 
qui joue le rôle principal dans la dynamique de 
la plupart des perturbations atmosphériques 
d'échelle synoptique, il n’en reste pas moins 
certain que l'apparition et même le développe- 
ment de ces perturbations ne peuvent s’expli- 
quer que si l’on tient compte des phénomènes 
ayant leur siège dans les couches baroclines 
situées de part et d’autre, et en particulier des 

hénoménes de la troposphére moyenne et 
inférieure. Cela est évident pour les ondes 
longues quasi-stationnaires associées aux centres 
d'action, et qui sont manifestement comman- 
dées par la répartition des continents et des 
océans, ce qui ne peut se concevoir que par 
l'intervention des basses couches de l’atmos- 
phére. Mais cela n’est pas moins certain pour 
les perturbations cycloniques extra-tropicales, 
lesquelles affaraissent essentiellement là où une 
masse d’air froid d’origine continentale ou 
polaire parvient au-dessus d’une région mari- 
time ou continentale relativement chaude: c’est 
ainsi qu'en Méditerranée les perturbations 
d’hiver se forment toujours lors d’une invasion 
d’air froid venant du Nord, et SVERDRUP (1943) 
a montré de même que la plupart des pertur- 
bations de front polaire nord-atlantique peu- 
vent s'expliquer par un phénomène analogue 
dans la région des eaux tropicales du système du 
Gulf-Stream. On comprend ainsi pourquoi les 
fronts polaires sont surtout actifs en hiver, 
lorsque les masses d’air froid sont particulière- 
ment bien développées, et pourquoi c'est 
également en hiver qu'ils s’approchent en 
moyenne le plus de l'équateur. En été, c’est 
au contraire surtout sur les continents que l’on 
observe la formation de perturbations cyclo- 
niques, et alors il est naturel de faire appel, pour 
les expliquer, à l’arrivée d’une masse d’air 
frais océanique sur le sol relativement chaud 
des continents. 

En somme, l’un des facteurs déterminants 
de la cyclogénèse, ce serait non pas le conflit 
entre deux masses d’air juxtaposées de tem- 
pératures différentes, mais bien plutôt le conflit 
entre une masse d'air froid et un sol chaud, que ce 
sol soit de nature maritime ou continentale; ce 
serait donc un facteur d’instabilité thermique 
ayant son siège dans la troposphère inférieure. 
Nous aurons à nous baser sur cette conclusion 
dans le chapitre sur les ondes de cyclone, et à 
chercher alors sous quelle forme on peut 
concevoir l'interaction nécessaire qui existe 
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entre la troposphére inférieure et la couche 
barotrope du jet-stream. 

Ainsi que l’ont fait la plupart des auteurs dont 
nous avons mentionné ci-dessus les travaux, 
nous considérerons des perturbations de faible 
amplitude (infinitésimales) superposées à un 
courant de base zonal, et telles que le mouve- 
ment perturbé soit non seulement horizontal, 
mais en outre sans divergence, car l’effet de la 
compressibilité de l'air est négligeable à 
l'échelle synoptique s’il n’y a pas de mouvement 
vertical. L'étude du mouvement peut alors 
être faite sur une surface de niveau déterminée, 
sans qu'il soit nécessaire de rien connaître de 
l'atmosphère en dehors de cette surface. Il est 
important de remarquer qu’une telle hypothèse 
relative au mouvement perturbé n’est possible 
que grâce à la barotropie du mouvement de 
base: nous avons en effet montré dans un 
travail récent (QUENEY 1950) que le fait de 
négliger a priori la vitesse verticale dans les 
équations de perturbation conduisait à des 
contradictions mathématiques dans le cas d’une 
atmosphère barocline. 


Nous prendrons comme perturbations élé- 
mentaires ce que nous avons déjà appelé des 
«ondes zonales», c'est-à-dire telles que la per- 
turbations de toute grandeur hydrodynamique 
soit donnée en fonction du temps f et de la 
longitude x par la partie réelle d’une expression 
de la forme 


À (œ) ein (x—œt) (1) 
où # est le nombre d’onde (entier positif), & 
la vitesse de phase angulaire complexe (con- 
stante), i l'unité imaginaire, et A (y) la demi- 
amplitude complexe, supposée infiniment petite 
et fonction seulement de la latitude @ sur la 
surface de niveau considérée. Dans la suite c’est 
toujours de ces ondes que nous nous occupe- 
rons, et quand nous parlerons de la perturba- 
tion d'une grandeur, il sera entendu qu’il 
s'agira de son expression complexe de la forme 
(x), et non de la partie réelle de cette expression. 

En plus des paramètres n et &, nous intro- 
duirons aussi la longueur d’onde angulaire A, 
definie comme l’inverse du mombre d’onde, 


u (2) 


et la pulsation orbitale complexe N, définie en 
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fonction de la vitesse angulaire © du mouve- 
ment de base par la formule 


N =n (© — 2) (3) 


La partie réelle de cette quantité représente la 
pulsation orbitale vraie, tandis que la partie 
imaginaire donne le coefficient d’amortissement 
de onde, compté comme positif ou négatif 
suivant que cette onde est amortie ou amplifiée 
(instable) en fonction du temps, et égal par 
définition au taux de réduction de l'amplitude 
par unité de temps. Dans les formules (1) et 
(3) le systéme de référence peut étre indifférem- 
ment un systéme sans rotation absolu ou un 
système fixe par rapport à la Terre, mais afin 
d’avoir des formules simples nous supposerons 
toujours que c'est un système absolu, de sorte 
que & représentera la vitesse angulaire absolue 
du mouvement de base: si U désigne la vitesse 
linéaire par rapport à la Terre (comptée posi- 
tivement vers l'Est) et Q la rotation terrestre, 
on à ainsi 


U =1(w—2Q) (4) 


où r désigne la distance à la ligne des pôles. 

Nous admettrons que la surface de niveau 
considérée est assimilable à une sphère concen- 
trique à la Terre et de rayon R = 6.368 km, 
donc nous écrirons r = R cos y, et ainsi que 
l'a déjà proposé Kuo (1950) il nous sera 
commode de prendre comme coordonnée 
méridienne, au lieu de la latitude y, la quantité 


y définie par 


ou par les trois relations équivalentes 


I 


= tO) 


th y = sing; sh y = te y; ch y = re > 
Cette variable y est celle qui intervient lors- 
qu'on effectue une représentation conforme 
de la surface de niveau sphérique sur un 
cylindre tangent au cercle équatorial; on voit 
que y varie de — co à + © lorsque @ varie 


7 TT 
de — —A + = ; la correspondance numérique 
2 2 


entre y et p est indiquée sur l’&chelle des 
abscisses de la figure 1. Grace à l’emploi de 
cette coordonnée méridienne, on obtient 
l'expression simple suivante pour ¢, compo- 
sante verticale du tourbillon absolu: 


&=F—y=20thy— — (7) 


où F = 2m sin y est une sorte de paramètre 


Be yar et Ow 
de Coriolis généralisé, et y = — 


a différe peu 


du gradient méridien de U. 


Le mouvement étant supposé sans diver- 
gence, les composantes et v de la perturbation 
locale de la vitesse, suivant les directions Est et 
Nord respectivement, s'expriment en fonction 
d'un potentiel des vitesses y, par les formules 


vcosp = ov sing dcosp= = — (8) 


Ox 


et d’autre part, sil n’y a pas de viscosité, le 
mouvement satisfait à l’&quation du tourbillon 
D¢é/Dt =o. En linéarisant cette équation, et 
en remplaçant ensuite la perturbation locale 
de £ par son expression en fonction de u et v, 
puis u et v par leurs expressions (8) en fonction 
de y et de dy/dy, on arrive à l’équation diffé- 
rentielle suivante pour y: 


n OÙ 
2 5) yp =o (9) 


ou l'opérateur V? a pour expression 


ES LE 0? no ba a? 

HAUTE" iP Noy 
(dans ces formules, r et A¢/Ay sont relatifs au 
mouvement de base). Cette équation (9), déjà 
obtenue par Kuo (1950) avec des notations un 
peu différentes, constitue l'équation différen- 
tielle fondamentale qui détermine les propriétés 
des ondes; elle a d’ailleurs aussi été utilisée par 
plusieurs auteurs dans les cas où la courbure 
terrestre est négligeable, de sorte qu’alors x et y 
représentent des coordonnées cartésiennes ordi- 
naires, mais nous voyons que grâce à l'emploi 
de la variable y l'équation est générale et tient 
compte automatiquement de la courbure 
terrestre. Nous devrons lui adjoindre la for- 


(nv: 


nt) (10) 
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mule suivante, donnant en fonction de y la 
perturbation locale p de la pression, et que 
l’on obtient à partir de l'équation mécanique 
de perturbation relative à la direction du 


parallèle géographique: 


2 ET re 

ZI II 
(o désigne la densité dans le mouvement de 
ase, que l’on peut considérer comme une 
constante puisque ce mouvement est supposé 


barotrope. 


II. Choix d’un mouvement de base simple 
et conditions aux limites 


L’équation (9) indique que les propriétés des 
ondes dépendent essentiellement de la réparti- 
tion de £ sur la surface de niveau considérée, ce 
qui nous conduit à chercher une loi de réparti- 
tion simple de ce paramétre, aussi voisine que 
possible des répartitions réellement observées. 
Mais il est évident que la loi la plus simple est 
celle qui consiste à supposer ¢ uniforme, car 
alors on obtient tout simplement l'équation 
classique de Laplace 


et d’une autre côté, ainsi que l’a déjà montré 
PLATZMAN (1943), c'est précisément à la con- 
sidération de courants où ¢ est uniforme que 
l’on est conduit par les seules données empiri- 
ques. Rappelons à ce sujet que, dans l’une de ses 
premières publications relatives au jet-stream, 
RossBy (1947) avait déjà décrit un modèle de 
jet simple caractérisé par une valeur constante 
de £ du côté polaire de l’axe du jet, et par une 
discontinuité de Ë sur l’axe lui-même, corre- 
spondant à un point anguleux dans le profil 
méridien de la vitesse. Nous sommes ainsi tout 
naturellement conduit, autant par des raisons 
d'ordre mathématique que par des raisons 
d'ordre expérimental, à adopter comme 
modèle du jet-stream le système formé par un 
certain nombre de courants zonaux séparés par 
des lignes de discontinuité de & coincidant 
avec les axes, et tels que € soit uniforme dans 
chacun des courants au niveau considéré: en 
outre nous devons supposer que la vitesse elle- 
même est continue sur les axes: ainsi un jet 
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simple sera représenté par un système de deux 
courants de cette nature, et de même un 
système de trois courants analogues nous four- 
nira une représentation d’un jet double. 

Il ne peut d’ailleurs y avoir aucun incon- 
vénient sérieux à supposer que les courants 
extrêmes s'étendent jusqu'aux pôles, c’est-à-dire 
jusqu’à y= + co et y = — respectivement, 
si on détermine les solutions de l’équation (12) 
par la condition qu’elles sannulent pour ces 
deux valeurs de y: en effet ces solutions ont 
alors une valeur négligeable en dehors de la 
zone centrale du jet, tout comme les ondes 
réelles du jet-stream. 

En plus des conditions à l'infini que nous 
venons d’indiquer, les conditions aux limites 
à exprimer sont la continuité de la vitesse 
méridienne et de la pression sur chacun des 
axes du jet. Par suite de la continuité de w 
ces conditions équivalent à la continuité de w 
et de p, ce qui en vertu de (11) se traduit par 
les deux relations suivantes, pour chacun des 
axes: 

Ay = 0; À (dp/dy) = — (an’/N)y (13) 
en désignant de façon générale par AA la 
discontinuité algébrique d’une grandeur À sur 
l'axe considéré, et en posant 


G = /\C/2 (14) 


On note que, si on tient compte des formules 
(8), la seconde des conditions (13) donne 


as (15) 


d'où l'on déduit que le mouvement ondula- 
toire possède nécessairement une discontinuité 
de glissement sur chacun des axes. 

Dans la suite nous supposerons toujours une 
discontinuité comptée positivement lorsque 
l'accroissement de la grandeur considérée a 
lieu vers les y croissants, c’est-à-dire vers le 
Nord, de sorte que o sera toujours positit pour 
un axe le long duquel w est maximum, comme 
c'est le cas normal. 

Remarque. — Si au lieu de supposer € con- 
stant, on suppose w = const, on a aussi N = 
const et d’autre part Æ/dy = (w/Q) rB + rB, B 
étant le gradient méridien du paramétre de 
Coriolis. Si en outre on considére des ondes 


= 


a a, RE, 


dont l'échelle méridienne est grande par 
rapport à la longueur d’onde, tout en étant 
petite par rapport au rayon terrestre, on peut 
dans l'équation (9) négliger 22/9y? à côté de 
n?, et considérer rß comme une constante. On 
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ch 
© = Wo + Yo vs sh (y — Yo) + 
sh? (y — 0) 
0 ch2 Yo > (17) 


obtient ainsi une relation d’où l’on déduit 
facilement la formule de Rossby 


¢ = U— (8/4 x?) L? (16) 
ou L désigne la longueur d’onde linéaire et c 
la vitesse de propagation, toutes deux assimi- 
lables à des constantes. Bien que cette formule 
ne soit strictement valable que pour des ondes 
très courtes, on sait que Rossby l’a cependant 
appliquée même aux ondes d'échelle plané- 
taire. Nous verrons au chapitre V dans quelle 
mesure les résultats auxquels elle conduit pour 
de telles ondes sont comparables à ceux déduits 
de l'hypothèse d’un système de courants à 
tourbillon vertical uniforme. 


III. Ondes en l’absence de viscosité 


1. Jet simple 


Pour éviter toute complication dans la 
terminologie, nous conviendrons d'employer 
l’expression «et simple» pour tout système de 
deux courants du type précédemment défini, 
et celle d’«axe» pour la ligne de séparation de 
ces courants. Nous généraliserons aussi ces 
expressions pour un système d’un nombre 
quelconque de courants. Suivant que la 
discontinuité du tourbillon est positive ou 
négative, nous dirons que le jet simple est 
positif ou négatif. Dans le cas usuel d’un jet- 
stream d'Ouest, on a affaire à un jet positif, 
mais la théorie que nous allons faire est générale 
(un jet-stream d’Est serait un jet négatif). Nous 
désignerons par y, la valeur de y sur l’axe, et 
de façon générale l’indice zéro sera utilisé pour 
désigner la valeur de toute grandeur sur cet 
axe: en particulier la demi-discontinuité de & 
sera donc représentée par 09. Un profil typique 
de jet simple positif est fourni par la courbe 
A de la figure 1. L’équation d’un tel profil 
dans chacun des deux courants peut être 
obtenue à partir de l'expression (7) de £: 
l'intégration de cette équation par rapport à 
«æ donne en effet, lorsque & est supposé in- 


dépendant de y, 


où Wo et Yq sont les deux constantes d’intégra- 

tion, égales respectivement aux valeurs de w 

et de y sur l’axe et dans le courant considéré 
+ 

(oniaediailleurs)= 2.0, th, 170 —L). Une 

autre expression équivalente de w est 


© =W) + 2 sh 2 (y— yo) + = sh? (y— Yo) 
| (18) 


où Fo’ = (2@ /ch*yo) + 279 th yo désigne la va- 
leur de JF /dy surl’axe, d’ailleurs peu différente de 
Bor, (on a en effet Bory = 22/ch*y,). La formule 
(17) montre que le profil d’un jet simple est 
complètement déterminé lorsqu'on se donne 
le point sur l’axe (point anguleux) et les deux 
tangentes en ce point, ce qui fait au total 4 
paramètres indépendants. Notons aussi que la 


discontinuité de y sur l’axe est opposée à celle 


de ¢: 


NY =—- Ast = — 20 (19) 


comme cela résulte immediatement de la 
formule (7). 

Les solutions de l’&quation d’onde (12) qui 
sannulent à l'infini et qui sont en outre 


continues sur l’axe sont de la forme 


Y = Yo e—nly—yol (20) 
De 1a on tire la discontinuité de y sur I’axe, 
Ao = — 2ny, et par suite la seconde des 
conditions (13) fournit la relation simple 


(21) 


Comme on a d'autre part N, = n (wo — 2), 
on en deduit la formule de dispersion suivante, 
donnant la vitesse de phase angulaire « en 
fonction de n ou À considérés comme pouvant 
varier de façon continue: 


(22) 


Cette expression étant linéaire, toutes les ondes 
sont stables, et leur spectre de propagation 
(donnant la partie réelle de « en fonction de 4, 


& = Wy — (d9/n) = wg — 094 
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pour À > o) est donc une portion de droite, 
ou plus exactement c’est l’ensemble des points 
de cette droite qui correspondent aux valeurs 
entières positives de n = 1/4. On remarque 
que la formule de dispersion ne dépend que 
des deux paramètres wo et 69, et non de quatre 
comme le jet simple lui-même. 

Il y a accord de phase, au signe près de 
l'amplitude, pour tous les éléments du mouve- 
ment ondulatoire le long de chaque méridien, 
et on déduit de la formule (11) l'expression 
suivante pour p: 


(23) 


où & désigne la demi-somme des valeurs de & 
sur l’axe, et où le signe à adopter est celui de 


(Yo — y): 


Pr Ok [do + n (© — &5)] Y 


2. Jet double 


Conformément 4 la convention indiquée 
précédemment, nous désignerons par «et 
double» tout systtme de 3 courants du type 
considéré dans ‘cette étude, et par «axes» de ce 
jet les deux lignes de séparation des courants. 
Nous dirons qu'un axe est positif ou négatif 
suivant le signe de o sur cet axe, et nous dirons 
que le jet double est positif, négatif ou mixte, 
suivant que ses axes sont tous deux positifs, 
tous deux négatifs, ou de signes contraires. 
Nous désignerons par yo et y; les valeurs de y 
sur les axes (ces mémes notations serviront 

our désigner les axes eux-mémes), par d 
fee [Yo — y1l ou «pseudo-largeum du cou- 
rant médian, et de façon générale par les 
indices o et 1 les valeurs d’une grandeur sur 
Yo et yı respectivement, tandis que les indices 
a, b et c seront utilisés pour les trois courants 
respectivement, ces courants étant supposés 
rangés dans le sens de y, vers }1. 

Si on prolonge le courant médian au delà 
de l'axe y; par exemple, on obtient un jet 
simple d’axe yo, et de même son prolongement 
au delà de y donne un jet simple d’axe y,: ces 
deux jets simples sont ce que nous appellerons 
les gets simples composants du jet double 
considéré. On ne peut pas les choisir arbitraire- 
ment, car le profil de chacun d’eux doit 
satisfaire à la condition de passer par le point 
anguleux du profil de l’autre, donc un jet 
double dépend au total de six paramètres 
indépendants. Les spectres de propagation des 
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jets composants, constitués par deux droites 
D, et D,, seront appelés «spectres composants» 
du jet double. 

Cela posé, on établit facilement que toute 
solution de l’équation d’onde qui s’annule à 
l'infini et qui vérifie en outre les deux condi- 
tions (13) sur l’axe y, par exemple, est donnée 
par les formules 


Ya = Vo e »\v— vol Ye = Yi e—*ly—ml 
ee erp ia Chee 
Yo=Y1 (Fhe nly Ms aren Bi ai) 


(24) 


Si nous écrivons ensuite que les conditions sur 
Yo sont satisfaites, nous obtenons deux relations 
que nous pouvons mettre sous la forme 


No¥1 1 


Ni Yo 


Co end ee No — 09 (25) 
N, Fa 0] 0] end 


De là nous tirons la relation suivante, qui 
généralise (21): 


(No — 00) (Ni—91) = 001€ 2" (26) 
et en remplaçant N, et N, par leurs expressions 
en fonction de « on en déduit l'équation de 
dispersion des ondes, qui est cette fois du 
second degré. Il nous sera commode d’écrire 
cette équation sous la forme 


C=O Oa end 
0 x = sn (27) 
Oo O1 ne 


où % et x, représentent les vitesses de propaga- 
tion des ondes des jets composants, données 
par les formules 


(28) 


et qui sont représentées respectivement par 
D, et D, lorsqu'on porte A en abscisse. 

On voit que le spectre de propagation du 
jet double, défini par l'équation (27), est 
complètement déterminé lorsqu'on se donne 
les spectres composants D, et D,, et en outre 
le paramètre d, donc le spectre dépend de cinq 
paramètres au lieu de six. Voici maintenant 
plusieurs propriétés que l’on établit aisément: 

1°) Le système DoD, constitue la forme 
limite du spectre de propagation lorsque, ces 
droites restant fixes, d devient infini. 


% = Wy — (0o/f); &ı = @, — (0,/n) 
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2°) Le spectre est tangent à D, et D, sur 
l’axe des &, avec un contact d’ordre infini. 
3°) Le spectre et ses composantes possèdent 
À ERS: 
le même diamètre D relatif à la direction de 
1: d di \ ta = 
axe des &, ce diamètre ayant pour équation 


I 
oa. (to + &ı) (29) 


4°) Si sans changer d on intervertit D, et Dy, 
le spectre n’est pas modifié, et si on remplace 
le système D,D, par son symétrique par 
rapport à une parallèle à l'axe des A, le spectre 
tout entier subit la même transformation. 

5°) Le spectre admet en général deux 
asymptotes ayant D pour diamètre, et dont 
l’une est parallèle A laxe des À et a pour 
ordonnée 


Ooo © a) 
aes on (Se - 1 


Gees Oye Lo, 24) (30) 
Dans le cas particulier où oy + 0, = 0 (D est 
alors un axe de symétrie), on a deux branches 
paraboliques. 

L’équation (27) étant du second degré, 
fournit deux valeurs de « pour chaque longueur 
d'onde, et dans certains cas ces valeurs peuvent 
être imaginaires conjuguées, ce qui correspond 
à des ondes instables ou amorties. Quatre cas 
sont à distinguer, suivant le signe de o, 0, et 
celui de l’abscisse Ap du point d’intersection 
de D, et D,, donnée par la formule 


O9 —i 


ln = (31) 


Cet 


Lorsqu'on a Ap > 0, les ondes des deux jets 
composants ont la méme vitesse de propagation 
pour À = Ap. 

Ier Cas: 6001 > 0; Ap < 0 (fig. 1, courbe a). 
Les racines de (27) sont manifestement réelles 
quelque soit n et comprennent % et %, dans 
leur intervalle, donc on a deux familles d’ondes 
stables, et le spectre de propagation (fig. 2a) 
se compose de deux courbes distinctes séparées 
Par, D; et IDE EE possédant chacune une 
branche infinie hyperbolique. En outre chaque 
courbe a pour position limite, pour d infini, 
la droite D, ou D, à laquelle elle est tangente, 
et elle est constamment inclinée dans le méme 
sens que cette droite par rapport aux axes A 
et x, ce qu’on peut établir de la façon suivante: 


méridiens de 


Fig. I. Profils systemes de courants 
à tourbillon vertical uniforme: A = jet simple; a, 
b = jets doubles positifs; c, d = jets doubles mixtes. 


L’équation (27) étant écrite sous la forme 


X— Ky X—%] 


Vby by = e-"4/n, où bo Og 


0 O71 


(32) 


sa dérivation par rapport 4 n donne 


2 % + &ı\ 0a 
& 
0,0, 2 on 


+ b—2 bo by (x + nd)end 


a ee 


(33) 


Si les deux quantités de même signe by et b, 
sont négatives, cette relation prouve directe- 
ment que le premier membre est négatif, et si 
elles sont positives on en déduit, en remarquant 
que l’on a 1 + nd < end, 


% + %\ Ja = 
(=) erates (Vbo — 


“a 


2n? 


vb,)? 
(34) 
donc cette fois le premier membre est positif. 


4 I 
Comme les trois quantités bo, bi et «— = 


Oo 01 


(to + 1) ont nécessairement un signe bien 
déterminé pour chacune des courbes du spectre, 
et comme ce signe ne varie pas lorsque d 
varie et que D, et D, demeurent fixes, on en 
déduit que la pente de chaque courbe est de 
même signe que celle de sa position limite pour 
d infini, c’est-à-dire que celle de D, et Dy. 

Il nous sera commode de désigner par les 
termes «ondes supérieures» et «ondes inférieures» 
les ondes relatives respectivement à la courbe 
supérieure et à la courbe inférieure du spectre 
de propagation. 


Fig. 2. Spectres de propagation des ondes des systèmes de courants a, b, c, d de la fig. ı.: 
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«=O 


spectres en 


trait continu fort, spectres composants D, et D, en trait continu faible, asymptotes en trait interrompu, axes 
de symétrie oblique en trait ponctué. 


2ème 


courbe b). 


Le seul changement par rapport au cas précé- 
dent consiste dans le fait que, pour d infini, la 
position limite de chacune des courbes du 
spectre de propagation (fig. 2b) est constituée, 
non plus par Dy ou D,, mais par l’une des 
lignes brisées formées par le système de ces 
deux droites d’un même côté du diamètre D, 
de sorte que lorsque d est grand le spectre 
présente une très forte courbure au voisinage 
de la longueur d’onde Ap. 


CAS 20004 > 0. 4m 10) (tin ts 


ReMe Cast. OeGs 103 Ap Ode 
courbe c) = Jet mixte du type stable. 


Dans ce cas, le profil méridien a une forme 


nettement différente de celle des cas précédents; 
malgré cela on a encore deux familles d’ondes 
stables quelque soit la longueur d’onde, mais 
cette fois les deux racines de (27) sont comprises 
entre & et %. Le spectre de propagation (fig. 
2c) se compose de deux courbes distinctes 
ayant chacune une branche infinie et chacune 
tangente à Dy ou D,. Le fait qu'il n'y a pas 
de racines imaginaires s'établit facilement à 
partir de la propriété sur le signe de la pente 
démontrée dans l’étude du premier cas et qui 
est générale: les deux courbes devant toujours 
être inclinées en sens inverse, puisqu'il en est 
ainsi de D, et Dy, sont toujours nécessairement 
séparées par le système des deux parallèles à 
l'axe des À d'équations « = wo et & = w,, et 


DS) à ND a a a it 


LE d A 


Ale. 
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par suite elles ne peuvent pas se couper; or 
si x pouvait être imaginaire l'équation (27) 
aurait une racine double pour certaines valeurs 


de A. 


HOME cas, 0,0, 0; Ap 7.0 (fig. 5; 

courbe d) = Jet du type instable. 
Ce cas est le plus intéressant car c’est lui qui 
nous conduira à une théorie des ondes de 
cyclone. En effet il existe alors une bande 
spectrale d'ondes instables entourant la lon- 
gueur d'onde Ap, car pour À = Ap onag%y = &ı 
et les racines de (27) sont visiblement imagi- 
naires. 

La propriété sur le signe de la pente du 
spectre permet d’ailleurs d’etablir facilement 
que la bande spectrale d’instabilité est unique 
et que sa largeur croît régulièrement à partir 
d’une valeur nulle lorsque d décroit à partir de 
+ co. Ce même résultat peut aussi s’obtenir 
en discutant le signe du déterminant de 
l'équation (27) qui est alors égal au produit des 
deux facteurs réels &, — & + 24V — 0,0, 
(e—"d/n), et on trouve en outre que la bande 
d'ondes instables a en général une largeur finie, 
sauf si on a og + 6, = 0 et en outre d< Ap. 

Le spectre de propagation (fig. 2d) comprend 
encore en général deux courbes distinctes, 
séparées par l'intervalle d’instabilité, mais c’est 
la même courbe qui est ici tangente à D, et Dj, 
l’autre courbe possédant les deux branches 
infinies (dans le cas particulier que nous 
venons de mentionner, cette dernière branche 
cesse d’exister); en outre le spectre doit être 
complété, dans l'intervalle d’instabilité, par la 
portion du diamètre D joignant les sommets 
des deux courbes précédentes. 

Remarque. — Dans tout ce qui précède nous 
avons supposé que n était un paramètre pouvant 
varier de façon continue de o à + ®, alors 
qu'en fait il ne peut prendre que des valeurs 
entières. Il en résulte en particulier que dans 
le dernier des cas considérés il n’existe effective- 
ment des ondes instables que si la bande 
spectrale d’instabilité est assez large pour 
contenir au moins un nombre d’onde entier. 
Toutefois nous verrons que cette restriction 
est sans conséquence pour la théorie des ondes 
de cyclone. 


Jet quasi-simple. 


Si l’une des deux discontinuités du tourbillon 
est petite, o, par exemple, le jet double est 


voisin du jet simple composant J, d’axe yj, 
et nous dirons alors que c’est un jet «quasi- 
simple», d’axe principal y, et d’axe secondaire 
Yo. Dans ce cas on a approximativement 
Ap = Ag, en posant 


Ay = (@1 —@))/0, (35) 


et si le jet est du type instable, ce qui est le 
cas Si On 212050; < 0° et A, > 0, la bande 
spectrale d’instabilité est étroite, donc si on 
emprunte le langage de la spectroscopie on 
peut dire qu'elle constitue une raie spectrale. 

Dans le cas particulier où o, est suffisamment 
petit pour que l’on puisse écrire 


Voo < Yard (no = 1/20) (36) 


on obtient des formules simples pour la forme 
de cette raie: si on désigne par 64 sa demi- 
largeur et par Sy la valeur du coefficient 
d'amplification en son centre Ay, où il est 
maximum, on a en effet 


OA 35 Euer 
I ee V 01 end; Sy = V— 090, € Mod 
(37) 


La première de ces formules montre que, pour 
toutes les ondes instables, on a A— 4, < À, 
e—2tod, Inversement, si cette dernière condition 
est satisfaite en même temps que la condition 
(36), et si l’on introduit les variables réduites 


fete 
j 


=; 


TER A (38) 


on trouve que, quelque soit le signe de oy 04, 
les deux valeurs de N, et celles de (N, — 0) 
qui correspondent 4 une longueur d’onde 
donnée s’expriment uniquement à l’aide de ces 
deux paramètres par les formules 


donc ces deux quantités sont petites par rapport 
Ao, e2"0d, et on voit en outre qu’au voisinage 
de la longueur d’onde A, le spectre de pro- 
pagation est assimilable à une hyperbole ayant 
D, et D; pour asymptotes, et que la courbe 
donnant le coefficient d'amplification des ondes 
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instables en fonction de la longueur d’onde est 
assimilable à un demi-cercle. Enfin des for- 
mules (15) et (39) on déduit, pour la discon- 
tinuité de glissement A ou sur l’axe secondaire yg 


(40) 


ce qui montre que cette discontinuité de 
glissement est négligeable, et en outre pour 
A =, ona 


Doug = 2 (Nie oc) ra € 29 <1 


(41) 


donc dans ce cas /\çu tend vers zéro en même 
temps que oy. Nous verrons l'importance de 
ce résultat lorsque nous considérerons l’effet de 
la viscosité. Le fait que la discontinuité de 
glissement est très petite sur l'axe secondaire 
permet également de conclure que, lorsque la 
condition (36) est satisfaite, le mouvement 
ondulatoire correspondant à chacune des ondes 
instables d’un jet quasi-simple est très peu 
différent de celui d’une onde du jet simple J, 
ayant la même longueur d’onde. 


Aoulvo = + 24 Vo0lo, emd 


3. Généralisations 


1°) Pour un jet multiple (système d’un 
nombre de courants supérieur A deux) l’équa- 
tion de dispersion des ondes est algébrique en 
& et son degré est égal au nombre des axes. 
Ainsi pour un jet triple ou système de quatre 
courants, comportant par conséquent trois axes 
Yo Yu Yo, on obtient l'équation du troisième 
; 
egré 


RX QE 
EN erzndo + | erzndıt 
O0 Oï 


“= ana) 
up) 


2, 
+ a en (do + di + do) — fe) (42) 


Où dy = |y1 — Val; di = lY2— Vols de = |Yo — 
Yıl, et OU a, a, &, désignent les vitesses de 
propagation des ondes des trois jets composants 
d’axes yo, V1 Ya respectivement. 
2°) Si un jet triple a ses trois discontinuités 
de même signe, toutes ses ondes sont stables. 
Supposons par exemple que, dans le jet 


triple précédent, oo, 61, et G2 soient tous trois 
positifs. On peut toujours, vu l’interchangeabi- 
lité des indices, supposer que l’on a yo < 1 
< yspar exemple, de sorte que si l’on prolonge 
le courant compris entre y, et y, au delà de 
Ya, on obtient un jet double positif d’axes y, et 
y,, dont les ondes vérifient l'équation de 
dispersion (27). Comme il s’agit d’un jet 
positif, cette équation a ses deux racines réelles 
et distinctes, et si «’ est la plus petite et « la 
plus grande, nous savons que l’on a « < «5; 
ae, et 7 0-91 f(a) destene 
le premier membre de (42), on établit facile- 
ment d’autre part, en utilisant la formule (32), 
que l’on a 


3 Ve 


et comme on a f(— ~) <o et f(+ æ) > 0, 
on en déduit que les racines de (42) sont bien 
toutes trois réelles, et qu’en outre elles sont 
alternées avec celles de (27). 

Cette propriété s'étend par récurrence à un 
jet multiple quelconque dont toutes les dis- 
continuités du tourbillon sont de même signe. 

3°) Si Pune des discontinuités o, d’un jet 
triple d’axes yo, Yı, Ya est supposée petite, le jet 
triple est voisin du jet double J;, défini par 
les six paramètres yı, Yo, ©1, Oo, 01, Cat NOUS 
dirons que c'est un jet «quasi-double» d’axes 
principaux y, et y, et d’axe secondaire yy. Dans 
le cas où o, et oy sont de même signe, et où 
Gy est de signe contraire, il peut exister des 
raies spectrales d’ondes instables, au nombre 
de deux au maximum, et dont les longueurs 
d'onde centrales sont fournies par l'intersection 
du spectre de propagation des ondes du jet 
Je avec la parallèle à l'axe des À d’équation 
% — @g. Pour chacune des ondes instables la 
vitesse de propagation est voisine de la vitesse 
du vent à la latitude de l'axe secondaire yo, 
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et en outre la discontinuité de glissement sur 
cet axe est très petite, de sorte que le mouve- 
ment ondulatoire est peu différent de celui 
d'une onde du jet Jıs ayant la même longueur 


d'onde. 


Ces résultats peuvent être étendus À un jet 
multiple quelconque dont l’une des disconti- 
nuités est petite et de signe contraire à toutes 
les autres supposées de même signe. 


IV. Effet de la viscosité 


Avant d'appliquer les résultats précédents, 
il est indispensable d'examiner quel peut être 
l'effet de la viscosité turbulente sur les ondes 
obtenues. En effet, même si cette viscosité est 
faible, ce qu'il est permis d'admettre en ce qui 
concerne la couche barotrope, il est bien comme 
que son rôle peut être néanmoins prépondérant 
lorsque la vitesse de propagation des ondes est 
voisine de la vitesse du vent au voisinage d’une 
discontinuité de glissement, où les forces de 
viscosité sont nécessairement très grandes. 
Notons d’ailleurs que si « =«, c’est-à-dire 
si N =o, l'équation générale (9) comporte 
un point singulier, donc pour avoir une 
solution déterminée et physiquement possible 
il est nécessaire de faire intervenir la viscosité 
(ou tout autre facteur de dissipation) au 
voisinage de la valeur correspondante de y, 
ainsi que l’a montré Kuo (1949 et 1950). 
Toutefois si 9&/dy = 0, comme nous l'avons 
supposé, de tels points singuliers cessent 
d'exister, de sorte que seul est à considérer le 
voisinage des axes (c’est A encore un avantage 
de l'hypothèse que & est uniforme). 

Afin de rendre le problème accessible au 
calcul, nous allons supposer que la viscosité 
turbulente est du même type que la viscosité 
moléculaire, caractérisée par un certain coeffi- 
cient cinématique v. On trouve alors que la 
formule (11) donnant la perturbation de la 
pression doit être remplacée par 


ae «| y (44) 


et que si » est une constante, ce que l’on peut 
supposer au voisinage d’un axe du jet-stream, 
l'équation du tourbillon devient D¢/Dt = 
= v/\2C, ce qui conduit à l'équation différen- 


tielle suivante pour y, à substituer à (9), et qui 
est cette fois du quatrième ordre: 


à 9 
[W+ivvy v4 55] v=o (as) 


Comme nous supposons & indépendant de y, 
cette équation se réduit à 


NST PAS ETATS (46) 


D'autre part les conditions aux limites sur 
les axes doivent être également modifiées: on 
doit exprimer d’une part la continuité de u 
et », et d'autre part celle des deux composantes 
horizontales de la tension de viscosité exercée 
le long de l’axe considéré par l’un des courants 
sur l’autre, ce qui fournit quatre conditions au 
lieu de deux. Si on tient compte de (8) et de 
(44), on trouve que ces conditions sont 
équivalentes aux suivantes: 


ae CHE DN N, 
Ava) nal) ne 


a 


Py Cae 
= (i) Ace 


1. Jet simple. — Lorsque v est très petit, la 
modification du mouvement par la viscosité 
est localisée dans une «zone de frottement très 
étroite entourant l’axe yo, donc dans cette zone 
on peut considérer v, r et .N comme des 
constantes (v= vo, f= 19) N = No): Legu- 
tion (46) devient ainsi une équation à coefficients 
constants, qui avec les conditions (47) redonne 
la formule N, =o obtenue pour v =a, 
moyennant la condition 


(47) 


(48) 


ou lon-pose ky =mir, = 2r/L,, Ly étant la 
longueur d’onde linéaire à la latitude de l’axe. 
En outre les solutions obtenues indiquent que 
dans la zone de frottement la quantité n(y — yo) 
est de l’ordre du premier membre de (48), donc 
la largeur de cette zone augmente si on fait 
croître v, sans faire varier ky et Go, et alors il 
arrive forcément un moment où N ne peut 
plus être considéré comme une constante dans 
la zone de frottement. Pour déterminer à quel 
moment cela se produit, nous remarquons que, 


Fig. 3. Effet de la viscosité sur les ondes d’un jet 

simple lorsque la discontinuité du tourbillon est pe- 

tite: courbes de propagation (graphique supérieur) et 
d’amortissement (graphique inférieur). 


d’après (18), ou peut écrire au voisinage de 
l'axe, dans chacun des deux courants, 


I 
N = No + Yot (Y — Yo) + > Fon (y— yo)? 
(49) 


Yo étant la valeur de y sur l’axe dans le courant 
considéré, de sorte que nous avons à chercher 
quand le second ou le troisième terme de ce 
développement devient comparable à 64. 
Comme n(y — yo) est de l’ordre de (ko? ¥9/a9)!, 
et comme d’autre part F,’ différe peu de Boro, 
nous voyons ainsi que c'est lorsque v, devient 
comparable à l’une des deux quantités o%/y02k? 
ou Og2/Poko. Or, avec les ordres de grandeur 
usuels suivants (en c. g. s.): 


Gg™ 10>; y, mw 10°; By m 1015: ky ~ 1078 


ces deux quantités sont de l’ordre de 10", ce 
qui est trés grand pour une viscosité turbulente 
dans la couche barotrope. En conséquence, 
nous pouvons dire que dans le cas d’un jet simple 
normal leffet de la viscosité sur l'équation de 
dispersion des ondes est négligeable, et que cet effet 
se réduit a une modification du mouvement dans 
une couche de frottement très étroite (comparée à 
la longueur d’onde). 

Naturellement cette conclusion cesse d’être 
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correcte si 09 OU Yo est nettement inférieur 
à l’ordre de grandeur 10”, et si 66 devient 
très petit l’effet de la viscosité doit au contraire 
étre prépondérant et produire un amortisse- 
ment des ondes qui augmente indéfiniment 
lorsque o, tend vers zéro: en effet à limite on 
a un courant unique où & est uniforme, donc 
l'équation différentielle (46), qui est du 
quatrième ordre, n’admet aucune solution non 
nulle s’annulant à l’infini et qui soit partout 
continue, en d’autre tèrmes il ne peut y avoir 
aucune onde du type considéré. Pour préciser 
ce résultat, supposons que y, conserve l’ordre 
de grandeur 105, et que seul o, soit petit; 
alors, comme les deux valeurs de y, diffèrent 
de 20, on peut considérer que ces deux valeurs 
sont égales, et on trouve que la viscosité 
intervient dans l'équation de dispersion des 
ondes lorsque o, est au plus comparable à la 
quantité 


(so) 


qui est seulement de l’ordre de 10° au plus 
pour vy = 10°, si on suppose Ly > 1.000 km. 
En outre le développement (49) peut étre 
limité 4 ses deux premiers termes, ce qui 
conduit pour (46) 4 une équation d’Orr- 
Sommerfeld, dont les solutions s’expriment 
à l’aide des fonctions de Bessel d’ordre 1/3. 
En écrivant que les conditions aux limites sont 
satisfaites, on arrive finalement, compte tenu 
de l’ordre de grandeur des paramètres, à une 
équation qui ne dépend que des deux variables 
réduites 


Uo = (Yo? ko? vo)s 


g(—iN) =— 


i 
9 
où g est la fonction définie par 


oo Bp 
J exp (2 =) dt (53) 


ce} 


Cette équation (52) détermine N, en fonction 


de o, par l’intermédiaire de Ng; on trouve 

, A er Ë ‘ 
qu'elle admet une infinité de solutions, qui 
correspondent toutes à des ondes infiniment 
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Fig. 4. Effet de la viscosité sur les ondes d’un jet quasi-simple 


IOI 


1ère solution principale 
2ème 


” ” 


pour les longeurs d’onde voisines de Ay: 


courbes d'amortissement & = const, chiffrées suivant les valeurs de £ (paramètre réduit proportionnel à l’écart 


oe Way. 


amorties lorsque 6, = 0 et dont l’une d’elles 
(solution «principale) est asymptotique à 
No =% pour 05 > fo. Certaines de ces 
solutions sont représentées en fonction de 
Go par les courbes de la fig. 3, dont la partie 


supérieure donne la partie réelle R(No), et la 
partie inférieure la partie imaginaire i (No), 
fournissant une valeur réduite du coefficient 
d’amortissement; la solution principale est 
figurée en trait continu, les autres solutions en 
traits diversement interrompus suivant la solu- 
tion dont il s’agit. On voit que pour la solu- 
tion principale les ondes sont amplifiées si l’on 
a Gy > 1,65, mais que le coefficient d’amplifica- 
tion est très petit ct pratiquement négligeable. 


2. Jet double — Tout comme pour un jet 
simple, on établit que pour un jet double 
normal, où do, 01, Yo et y, sont tous de l’ordre 
de grandeur 10 c. g.s., l’effet de la viscosité 


sur l'équation de dispersion des ondes est 
négligeable; au contraire cet effet consiste en 
général en un amortissement si par exemple 
l’une des demi-discontinuités o, du tourbillon 
est très petite, auquel cas on a affaire à ce que 
nous avons appelé un jet quasi-simple. Toute- 
fois ce dernier résultat comporte une exception 
importante, dans le cas où la longueur d’onde 
À est voisine de la valeur À, définie par (35) 
et qui, en l’absence de viscosité, constitue le 
centre de la raie spectrale des ondes instables 
éventuelles: on trouve en effet que, dans ce 
cas, l'équation qui détermine N, est d’un type 
autre que (52) et que, moyennant certaines 
conditions, elle peut se mettre sous la forme 
- (£— No) 


g (—iNo) = (54) 


GUESS 


dépendant des deux paramètres & et 7 défini par 
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Fig. 5. Effet de la viscosité sur les ondes d’un jet 
quasi-simple, pour les longeurs d’onde voisines de Ap: 


courbes d’amplification I(N,) = const, chiffrées suivant 


les valeurs de — I(N),). 

Les hachures continues indiquent les domaines des 
ondes instables fournies par la première solution prin- 
cipale, les hachures interrompues ceux des ondes four- 
nies par la seconde solution principale, et la parabole 
en pointillé représente la frontière du domaine des 

ondes instables pour une viscosité nulle. 


a we Vente oy es T __ Go 
Mo Ho À 


(ss) 
(€ et t sont définis par (38)). 


Cette équation (54) admet encore une infinité 
de solutions, mais parmi elles il en existe une 
(«première solution principale) qui devient 
identique à celle que nous avons obtenue en 
négligeant la viscosité dans le cas où odo est 
petit par rapport 4 u,, Contrairement à ce qui 
se passe dans le cas général où l’on a alors, au 
contraire, un amortissement trés grand. Ce 
résultat était d’ailleurs 4 prévoir, puisque nous 
avons établi dans le chapitre précédent (formule 
(41)) que la discontinuité de glissement tendait 
vers Zéro en même temps que Cp pour À = Ay. 
Donc, en particulier, l'effet de la viscosité sur 
les ondes instables d'un jet quasi-simple est 
pratiquement négligeable, résultat qui nous 
permettra d'appliquer sans aucune modification 
les formules du chapitre précédent A une 
théorie des ondes de cyclone. En plus de la 
première solution principale; il existe une 
autre solution de l'équation (54) («seconde 
solution principale) qui possède des pro- 
priétés asymptotiques intéressantes, sur les- 
quelles nous ne croyons pas utile d’insister ici. 
La fig. 4 donne les variations du coefficient 


d'amortissement réduit I (No) en fonction du 
paramètre VT pour 0, 0, > 0 (partie droite de 
la figure), et en fonction de V/-T pour 0,0, <0 
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(partie gauche), pour diverses valeurs de & et 
pour les deux solutions principales, et la 
fig. 5 montre quel est l’effet de la viscosité 
sur la largeur de la raie spectrale d'ondes 
instables pour 09 01 < 0. (On suppose 0, > 0). 


V. Application aux centres d’action 
semi-permanents 


Dans les chapitres précédents il n'a été 
question que des ondes libres du jet-stream, 
puisque nous avons supposé le mouvement 
strictement horizontal; mais il est bien connu 
que la théorie des ondes forcées d’un milieu 
fluide se rattache étroitement à celle de ses 
ondes libres: ainsi lorsque le fond d’une rivière 
peu profonde et de vitesse uniforme présente 
une petite dénivellation, il en résulte un 
système d’ondes de gravité forcées apparaissant 
à l’aval de la dénivellation et dont la longueur 
d’onde est égale à celle des ondes libres station- 
naires, si toutefois la vitesse du courant est telle 

re. 
que de telles ondes existent. Il s’agit A d’un 
phénomène de résonance, que l’on retrouve à 
plus grande échelle lorsqu'on considère par 
exemple les ondes produites dans un courant 
atmosphérique par une chaîne de montagnes 
peu accidentée. 

Si nous nous plaçons à l'échelle de la circula- 
tion générale, nous pouvons envisager de même 
des ondes forcées produites dans le jet-stream 
par les actions géographiques possédant cette 
échelle, ce qui nous conduit à une explication 
des ondes quasi-stationnaires associées aux 
centres d'action semi-permanents: c’est 1a le 
principe même de la théorie formulée par 
RossBy (1939), et nous voyons que le problème 
se ramène à la recherche des ondes libres du 
jet-stream qui sont stationnaires par rapport à 
la Terre, c'est-à-dire telles que l’on ait « = Q, 
ce qui nous oblige à supposer que le jet-stream 
lui-même est stationnaire. Toutefois la re- 
sonance dont il s’agit ici n’est certainement pas 
aussi stricte que dans le cas des ondes forcées 
À la surface d’une rivière, de sorte qu'il suffit 
que & soit voisin de Q, et non egal exactement 
à 2, pour que l’on puisse observer dans 
l'atmosphère des ondes longues de grande 
amplitude; en d’autres termes, il suffit que les 
ondes soient très lentes, et non pas exactement 
stationnaires, et on peut supposer aussi que le 
jet-stream lui-même, au lieu d’être rigoureu- 
sement stationnaire, possède une certaine 
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évolution pourvu qu’elle ne soit pas trop ra- 
pide. 

De façon générale, on obtient la condition 
pour qu'un jet-stream soit capable d’ondes 
stationnaires en remplaçant simplement « par 
Q dans l'équation de dispersion des ondes de 
ce jet-stream, ce qui fournit par conséquent 
une certaine relation entre les différents 
paramètres qui le définissent, pour chaque 
valeur de n. 

Avant d'appliquer ces considérations aux 
modèles de jet-stream étudiés précédemment, 
rappelons d’abord que, si on fait l'hypothèse 

"un courant de vitesse uniforme et d’une 
courbure terrestre négligeable, de sorte que 
l'équation de propagation des ondes est la 
formule (16), on obtient pour les ondes 
stationnaires les formules classiques de Rossby 


Xu PE 


n? 4 T°? 


(56) 


où ß et r sont assimilables à deux constantes, 
et où L = 2ar/n désigne la longueur d’onde 
linéaire. En donnant à n les valeurs entières 
successives, on obtient ainsi la série des valeurs 
de U pour les courants capables d’ondes 
stationnaires, et on voit en outre que pour ces 
ondes U est proportionnel à L?, de sorte que 
la relation entre U et L est représentée par une 
parabole (A). Cette parabole a été construite 
sur la fig. 6a en supposant B = 1,62 x 1071 


c. g.s., valeur relative à la latitude de 45°. 


1. Jet simple. — Si nous supposons le jet 
simple formé de deux courants du type 
considéré dans les chapitres antérieurs, nous 
devons appliquer la formule de propagation 
(22), d’où nous déduisons la condition de 
résonance suivante pour les ondes stationnaires: 


o 
op — © = = (57) 


déjà donnée par PLATZMAN (1949) avec d’autres 
notations; si nous introduisons la vitesse 
linéaire sur l'axe, Uy = ro(wo — 2), et la 
longueur d’onde linéaire mesurée à la latitude 
de cet axe, L = 2x ro/n, nous déduisons de 
la des formules correspondant à (56), 

Oot Gil 

NL 


n 27 


(58) 
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Fig. 6. Ondes stationnaires d’un jet simple (a) et 

d’un jet double (b), et comparaison avec la formule 

de Rossby: spectre de vitesse du jet-stream (portions 

de droites en trait continu) et relation entre la lon- 

geur d’onde et la vitesse du jet-stream (courbes en 
trait interrompu). 


Ces formules indiquent que pour un jet 
positif (6,9 > 0) on doit avoir U, > 0, autre- 
ment dit on doit avoir des vents d’Ouest au 
moins dans la partie centrale du jet, tout 
comme l’exige la formule de Rossby. D’autre 
part les coupes méridiennes publiées par 
CRESSMAN (1950) par exemple indiquent que 
pour un jet-stream dont l’évolution est lente 
la demi-discontinuité du tourbillon est toujours 
peu différente de la valeur moyenne 3.5 x 10 ° 
c.g.s. sur chacun des axes, de sorte que le 
produit 6, ry varie lui-même très peu et peu 
être considéré comme une constante en 
première approximation. Les formules (58) 
fournissent alors une valeur bien déterminée 
de U, pour chaque nombre d'onde n, et 
indiquent en outre que pour les ondes station- 
naires U, est proportionnel à la longueur 
d'onde L, et non plus à L?, donc la relation 
entre U, et L est représentée par une droite (B). 
Sur la figure 6a cette droite a été construite 
en supposant l’axe à la latitude de 45°, et en 
prenant pour o, la valeur précédemment 
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Fig. 7. Ondes stationnaires d’un jet simple ne com- 
portant que des vents d'Ouest (n = 3). 


indiquée; en outre les valeurs de U, corre- 
spondant aux valeurs entières successives de 
n sont mises en évidence sous la forme d’un 
«spectre de vitesse», constitué par des segments 
de parallèles à l’axe des L et se terminant sur 
la droite (B). 

Les formules (58), associées à la formule de 
propagation (22), conduisent à des résultats 
qualitativement identiques à ceux obtenus par 
Rossby: 

(a) Les ondes stationnaires sont d’autant plus 
longues que le jet-stream est plus rapide, 
donc que l’indice zonal est plus élevé. 

(b) Pour un jet-stream donné, les ondes plus 
courtes que les ondes stationnaires se meuvent 
vers l'Est, les ondes plus longues vers l'Ouest. 

En outre il est remarquable que la droite 
(B) est dans l’ensemble assez voisine de la 
parabole (A), ce qui signifie que les résultats 
déduits de l'hypothèse de Rossby song égale- 
ment comparables quantitativement à ceux que 
donnent les formules (58). 

L’aspect des lignes de courant stationnaires 
est donné par les fig. 7 et 8, pour deux jets 
simples particuliers: dans le premier cas on 
an = 3, et le jet-stream ne comporte que des 
Ventsi cd Ouestin(dy = 3.5.0% 10 Foie 
d'où U, = 52.6 m/sec), tandis que dans le 
second cas, qui est celui du profil Afdetla 
fig. I, on am —4, et une zone de vents 
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Fig. 8. Ondes stationnaires d’un jet simple com- 
portant une zone tropicale de vents d’Est (n = 4). 


d’Est de part et d’autre du tropique (04 = 3,5 X 
10%; yo. =.0,9, dot. U, =.38,9,m/sec) Dass 
les deux cas on a supposé que l’axe déformé 
qui coincide avec une ligne de courant puisque 
le mouvement est stationnaire, a une demi- 
amplitude A, égale à 8 degrés de latitude, donc 
petite par rapport à la longueur d’onde con- 
formément à l'hypothèse toujours faite jusqu'ici 
(cet axe déformé est indiqué par une suite de 
pointes de flèches sur les figures); en outre on 
a représenté, sur les parties supérieures, d’une 
part le profil du jet simple en coordonnées 
p et U (courbe indiquée par «U»), et d’autre 
part les variations méridiennes de la demi- 
amplitude A, des lignes de courant et de celle 
A, des isobares, exprimées en pourcentage de 
Ao, et dans l'hypothèse où ces demi-amplitudes 
sont infiniment petites. 

Les formules donnant A, et A, sont, lorsque 
ces quantités peuvent être considérées comme 
comparables à Ao, 


À, _ Og tn(w—2) 


a (© + Q) sin o 
(59) 


où le signe à prendre est celui de (yp — y), 
d'où l’on déduit que A, et A, deviennent 
relativement grands lorsque U est petit, et cela 
explique pourquoi on trouve le maximum de 


U 
—9 on] v—yol: 


AY’ a Pr 
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ces quantités, dans le cas de la fig. 7, non 
pas au voisinage de l'axe, mais plutôt au 
voisinage de la latitude où U est minimum, 
c’est-à-dire vers 30°. En outre dans le cas de 
la fig. 8, le fait que U s’annule à deux 
latitudes entraîne que A, et A, sont très grands 
à ces latitudes (théoriquement infinis), donc 
les lignes de courant et les isobares cessent 
d’être des sinusoides dans les régions tropicales, 
et prennent l'aspect de courbes fermées 
formant un chapelet d’anticyclones alternant 
avec des dépressions. Tout cela est bien con- 
forme à ce qu'on observe, et par ces deux 
exemples on se rend compte aussi que, pour 
une même action perturbatrice, l'amplitude 
des ondes stationnaires est en moyenne d’autant 
plus grande que U est lui-même plus petit, 
donc que l'indice zonal est moins élevé, 
résultat déja énoncé par Rossby sur une base 
purement empirique. 


Enfin on voit sur les figures que À, et À, 
sont toujours assez peu differents, ce qui 
signifie que le mouvement est approximative 
ment géostropique (on a d’ailleurs exactement 
A} — A; à la latitude où A, est maximum, 
dans le cas où U est toujours positif). 


2. Jet double. — Si on fait « = Q dans 
l'équation de dispersion (27), et si on introduit 
les vitesse U, et U,, sur les axes, on trouve 
comme condition pour qu'un jet double 
possède des ondes stationnaires 
end 


U, I 
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donc si n est donné on obtient une relation 
entre les trois paramètres Up/o9 ro Ui/01 rı 
et d. Ainsi que nous l’avons déjà dit, les 
produits o, ro et 61 r\ sont peu variables pour 
un jet-stream qui évolue lentement, de sorte 
qu'en première approximation nous pouvons 
les considérer comme égaux à une même 
constante. D'autre part d’après Cressman la 
distance entre les axes est peu variable et de 
l'ordre de 15 à 20 degrés de latitude, ce qui 
nous conduit à assimiler également d à une 
constante. Enfin, comme Uy, et U, ne sont 
jamais très différents, nous pouvons aussi 
admettre que les deux paramètres sont égaux, 
et nous voyons alors que le nombre de para- 
mètres se crouve réduit à un seul Up, de sorte 
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Fig. 9. Ondes stationnaires d’un jet double ne com- 

portant que des vents d'Ouest, et capable de deux 

longeurs d’onde stationnaires différentes (n = 4 et 
iO) 


que la condition (60) établit, tout comme pour 
un jet simple, une relation bien déterminée 
entrerU, eC 71, à savolt 


d’où nous tirons 


o lo (1 + ena) 


n 


Nous obtenons donc deux groupes de solutions, 
et il est facile de voir que les ondes stationnaires 
correspondantes sont des ondes supérieures si 
on prend le signe + devant l’exponentielle, 
et des ondes inférieures si on prend le signe —. 
Si on désigne encore par L la longueur d’onde 
mesurée à la latitude de l’axe yo, on déduit 
de (62) une relation entre U, et L, qui est 
représentée par un système de deux courbes 
(C) et (C’), correspondant respectivement aux 
ondes supérieures et aux ondes inférieures; sur 
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Fig. 10. Superposition des deux ondes stationnaires 
de la fig. 9. 


la figure 6b ces courbes ont été construites en 
prenant pour 09 fo la même valeur que pour 
la figure 6a, et pour d la valeur 0,42 (valeur 
moyenne déduite des données de Cressman) ; 
le spectre de vitesse correspondant (ensemble 
des valeurs de U, pour les valeurs entières 
positives de n) est représenté de la méme 
facon, mais il est constitué cette fois par la 
superposition de deux séries de raies. 

Cet exemple particulier suffit à établir le 
résultat important suivant: Dans le cas d’un 
jet double rn général d’après Cressman), 
l'indice zonal ne suffit plus à déterminer, même 
approximativement, le type ondulatoire de la 
circulation générale, puisqu'on peut avoir deux 
valeurs totalement différentes de la vitesse 
axiale pour une même valeur du nombre 
d'onde; toutefois il reste vrai que les ondes 
stationnaires les plus longues (n —2) ne 
peuvent apparaître que dans des jet-streams 
rapides. De même il n’est plus exact que, pour 
un jet-stream donné, les ondes plus courtes que 
les ondes stationnaires se meuvent vers l'Est 
et les ondes plus longues vers l'Ouest, et on 
peut aussi avoir des ondes de même longueur 
que les ondes stationnaires et qui se meuvent 
soit vers l'Est, soit vers l'Ouest; des études 
comme celles de Namras-Crapp (1944) pour- 
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raient donner d’intéressantes confirmations à ce 
sujet. 

Du chevauchement des deux séries de raies 
du spectre de vitesse de la fig. 6b, il résulte 
que certaines des raies de l’une des séries sont 
très voisines de l’une des raies de l’autre série, 
comme c’est le cas pour les raies n = 6 et 
n = 7 voisines de U, = 23 m/sec, et pour les 
raies n = 5 et n =.6 voisines de U, = 28 
m/sec, et ce résultat signifie que pour un même 
jet double on peut avoir simultanément des 
ondes pratiquement stationnaires ayant deux 
longueurs différentes. D’ailleurs, il est toujours 
possible de trouver des jets doubles tels que la 
condition (60) soit exactement vérifiée pour 
deux valeurs différentes de n, car si on suppose 
par exemple que d est donné on obtient ainsi un 
système de deux équations qui déterminent 
U,jlozer, tet: Uileit rm Lestfhengtetaeisons 
précisément relatives à un tel jet double, 
capable à la fois d'ondes stationnaires n = 6 
et d’ondes stationnaires n = 4 (les premières 
sont des ondes supérieures, et les secondes des 
ondes inférieures). Sur la fig. 9 on a repré- 
senté, en plus du profil du jet, des ondes 
stationnaires n = 6 et n = 4 dans l'hypothèse 
où celles-ci existent seules, et où l’axe le plus 
déformé a une demi-amplitude égale à 6 
degrés de latitude: on voit que pour n = 6 
c'est l’axe méridional qui est le plus déformé, 
l’autre ne l’étant pratiquement pas, tandis que 
pour n = 4 c'est l'inverse qui se produit, et 
comme dans les deux cas le maximum de A, 
ou de À, est voisin du milieu du courant 
médian, il en résulte que l’amplitude des lignes 
de courant et des isobares décroit très brutale- 
ment à partir du maximum, vers le Nord 
pour n = 6, et vers le Sud pour n = 4. Cela 
se traduit par un phénomène de confluence- 
diffuence très prononcé en certains points du 
courant médian. 

Dans la fig. 10 on a supposé que les deux 
ondes stationnaire précédentes sont superposées, 
et on voit qu'elles interfèrent seulement dans 
le courant médian, et que le phénomène de 
confluence-diffluence est encore accentué en 
certains endroits. 

Les divers résultats peuvent permettre de 
s'expliquer des effets tout à fait comparables 
décrits en particulier par NAMrAs-CLAPP (1944), 
sans qu'il soit nullement nécessaire de faire 
appel à des mouvements verticaux comme 
certains auteurs l’ont fait, et il serait important 
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de vérifier que ces effets sont effectivement 
associés à une structure double du jet-stream. 

En somme, la théorie que nous venons de 
faire rend compte d’une façon très satisfaisante 
des principaux phénomènes en rapport avec 
les centres d'action semi-permanents. 


VI. Application aux ondes de cyclone 


Nous avons montré dans l’Introduction que 
les seuls faits d'observation conduisaient à 
considérer la cyclogénése comme due à l’asso- 
ciation d’un facteur thermique troposphérique 
et d’un facteur dynamique propre au jet-stream, 
du moins en général. D’autre part nous avons 
trouvé dans les chapitres II et IV qu’un 
système de courants à tourbillon vertical 
uniforme pouvait être le siège d’ondes instables, 
ayant un maximum d'amplitude à une latitude 
bien déterminée, et une longueur d’onde 
voisine d'une ou plusieurs valeurs également 
bien déterminées, lorsque l’une des disconti- 
nuités du tourbillon est petite et négative alors 
que toutes les autres sont au contraire positives. 
Ces propriétés des ondes instables rappellent 
suffisamment celles qui caractérisent les ondes 
de cyclone dans le cas d’un courant de per- 
turbations bien établi, pour qu'il soit permis 
de penser que le facteur dynamique ci-dessus 
invoqué n’est autre que l'instabilité dynamique 
acquise par le jet-stream sous l’effet d’une petite 
déformation susceptible d’être schématisée par 
l'apparition d’une discontinuité négative du 
tourbillon.1 

Il n’est d’ailleurs pas difficile d'imaginer une 
suite de phénomènes simples conduisant à une 
déformation de ce type. Considérons par 
exemple le cas des cyclones du front polaire 
nord-atlantique qui, ainsi que nous l'avons 
suggéré, ont très probablement pour cause 
initiale l’arrivée d’une masse d’air froid con- 
tinentale ou arctique sur les eaux tièdes de la 
région du Gulf-Stream. Le premier phénomène 
qui doit logiquement résulter de ce contraste 
thermique est la naissance et le développement 

de courants de convection verticale dans la 
masse d’air froid, puisque celle-ci tend à être 


1 Depuis la rédaction de cet article, nous avons pu 
montrer que des ondes quasi-instables, de longeur 
d’onde bien déterminée, peuvent aussi résultér de la 
modification du tourbillon vertical dans une zone 
étroite du jet-stream, ce qui paraît plus admissible que 
l'introduction d’une discontinuité du tourbillon. 
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instabilisée par le bas lorsqu'elle atteint les 
eaux océaniques plus chaudes qu’elle. Ces 
mouvements convectifs, dont l’existence est 
d’ailleurs révélée par l'observation directe, 
doivent s'étendre vers le haut d’autant plus 
facilement qu'ils sont accompagnés de forma- 
tions nuageuses, donc rien ne semble s’opposer 
à ce qu'ils atteignent le jet-stream lui-même, 
et comme ils sont équivalents à une viscosité 
turbulente lorsqu'on les considère à l’échelle 
synoptique, ils doivent avoir pour effet de 
réduire la vitesse du jet-stream par rapport à 
la Terre, cela dans une zone relativement 
étroite débutant un peu à l’aval de leur lieu 
de formation et allant plus ou moins loin vers 
l'Est. En première approximation on a donc 
une déformation zonale du jet-stream, et 
comme il y a tout lieu de penser que la réduc- 
tion de vitesse est maximum à une certaine 
latitude, on voit qu’elle est assez bien sché- 
matisée sous la forme de l’apparition d’une 
petite discontinuité négative du tourbillon, ce 
qui conduit à un jet quasi-simple ou quasi- 
double, suivant que le jet-stream initial était 
simple ou double: on obtient bien ainsi la 
déformation du type voulu pour avoir des 
ondes instables. Ainsi se trouverait précisé 
à son tour en quoi consiste l'association 
des facteurs thermique et dynamique dont 
nous avons parlé: ce seraient les mouvements 
convectifs qui établiraient le lien nécessaire 
entre les basses couches de l'atmosphère, siège 
de l'instabilité thermique, et les couches 
voisines de 12 km, siège de l'instabilité dyna- 
mique à l'échelle de la circulation générale. Il 
nous reste maintenant à examiner de plus près 
cette théorie au point de vue quantitatif. 


1. Ondes de cyclone d’un jet simple 


Pour fixer les idées nous raisonnerons sur le 
jet simple ayant pour profil la courbe À de la 
fig. 1. L'indice 1 étant affecté à la latitude de 
Vaxe;jon ayalors: 20,99 =45 75.0, = 
OST MIO SECTE 0 à SAN MIO COTES! 
nous supposons ce jet transformé en un jet 
quasi-simple d’axe principal y, et d’axe 
secondaire yp, nous savons que le nouveau 
système est du type instable si l’on à la fois 
By < 0 et My <w,, et que la longueur d'onde 
angulaire moyenne des ondes instables est alors 
Ay = (© — ©0)/01, d’où l’on déduit la longueur 
d'onde linéaire correspondante par la formule 
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Tableau 1. Les variations de L, et S, en fonction de ® 


27°6 20°0 13° Or ea 6°o 


ENGSOO NA OOOE ENFE TOO OX o 


Lo (km) oe 1000 7 4000 7 7537 
S 10° N 4,7 
is sy 5,09, 0A 2,35 DORE Co $ ne 
alt zZ — > 
(110 km) (440 km) | (2300 km) 
Po 45 7 48°3 53°7 63°1 72°8 74°0 
Ly (km) 0 7 100007 2750760074 2318700 2520000 ex Oo 
Spee x ay 1,0 
(c. 8. 5.) 1,78 1,070 I 0,22 on [e) 
<— Z—> 
(290 km) (600 km) 


Ly =27 1 Ay. On voit ainsi qu'il ne peut y 
avoir des ondes instables que si la latitude , 
de l’axe secondaire est dans l’un des deux 
intervalles (6°0, 45°7), (45°7, 74°0) limités par 
les latitudes des points du profil A où l’on a 
@ = @,; en outre dans chacun de ces intervalles 
L, est nul aux latitudes limites et maximum à 
la latitude où est minimum sur le profil A. 
D’autre part le coefficient Sy d’amplification 
maximum des ondes instables est donné par 
la formule (37) avec une précision absolue 
suffisante tant que |o,| est au plus égal au 1/10 
environ de o,: nous adopterons la valeur 
Gy = —0,3 X 10° sec, dont l’ordre de 
grandeur paraît correct pour caractériser une 
petite déformation du jet simple considéré. Les 
variations de Ly et S, en fonction de wy dans 
chacun des intervalles ci-dessus mentionnés 
sont données dans le tableau 1. 

Pour que les ondes instables ainsi déter- 
minées soient susceptibles de représenter des 
ondes de cyclone, il faut d’une part que leur 
longueur d’onde ne soit ni trop courte ni trop 
longue, d’autre part que leur coefficient 
d'amplification ait une valeur significative A 
l'échelle synoptique. 

a) Si la longueur d’onde trouvée est tro 
courte, il n’est pas possible de négliger la 
composante verticale du mouvement ondula- 
toire comme nous l’avons fait jusqu'ici, donc 
la stabilité hydrostatique de l'atmosphère 
intervient pour réduire et même supprimer 
l'effet de l'instabilité dynamique. Notre théorie 
est donc incapable de déterminer quelle doit 


être la longueur d'onde minimum des ondes 
de cyclone, et nous devons à cet effet nous 
adresser à l’observation: ainsi que nous l'avons 
déjà admis, la valeur Ly = 1.000 km semble 
raisonnable comme minimum. 

b) Si la longueur d’onde est trop longue, 
l'assimilation de la déformation du jet-stream 
4 une déformation zonale n’est plus admissible, 
puisque la zone de déformation a toujours une 
longueur finie en raison de l'existence des 
phénomènes de dissipation. Ainsi dans le cas 
du front polaire nord-atlantique, dont la 
longueur ne dépasse guère 8.000 km, on ne 
peut à proprement parler considérer les 
perturbations associées à ce front comme des 
ondes de cyclone que si L, est au plus de l’ordre 
de 4.000 km environ, de façon à avoir au 
moins deux ondulations le long de la ligne de 
déformation. Si la valeur obtenue pour L, est 
de l’ordre de 4.000 km il est clair que la 
stabilité dynamique de la portion non dé- 
formée du jet-stream doit interférer avec 
l'instabilité due à la zone de déformation, ce 
qui peut encore conduire à des perturbations 
instables si Lo n’est pas trop grand, mais ce 
seraient alors des perturbations dépourvues de 
tout caractère ondulatoire, du type des cyclones 
tropicaux par exemple: nous y reviendrons 
dans un moment. Enfin si L, est nettement 
supérieur à 4.000 km, toute instabilité dyna- 
mique doit disparaître, donc on ne peut plus 
avoir de perturbations cycloniques correspon- 
dantes. Ce raisonnement simple, qu’il y aurait 
évidemment lieu de préciser, nous conduit 
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donc à adopter Ly = 4.000 km comme maxi- 
mum admissible pour les ondes de cyclone. 

c) L’échelle de temps synoptique étant de 
l’ordre de 1 à 4 jours, le coefficient Sy n’a une 
valeur significative que s’il est de l’ordre de 
hOmo°useGe> au Mons. 

Cette discussion nous conduit donc à ne 
retenir finalement que les intervalles de varia- 
tion dep, où l’on a à la fois 1.000 km < Ly < 
4.000 km et Sy > 10% sec-!, comme définis- 
sant approximativement les «zones sensibles du 
jet simple où doit se trouver l’axe secondaire 
pour que des ondes de cyclone puissent se 
développer. Le tableau I montre qu'il en 
existe deux, désignées par les symboles Z et Z’, 
et situées respectivement au Nord et au Sud 
de l’axe principal. On voit que ces zones sont 
relativement étroites (la plus large n’a que 
600 km), et que leurs distances à l’axe principal, 
également indiquées sur le tableau I, sont 
aussi relativement faibles (moins de 300 km). 
On note d'autre part que la limite septen- 
trionale de la zone Z est déterminée par la 
condition que doit remplir So, tandis que la 
limite méridionale de Z’ est au contraire 
déterminée par la condition que doit remplir 
Lo; de 1a il résulte que les ondes de cyclone 
dues à l’excitation de la zone Z ont une valeur 
maximum nettement inférieure à 4.000 km 
(de l’ordre de 2.500 km seulement), et comme 
c'est surtout cette zone qui a des chances 
dintervenir s'il s’agit des cyclones de front 
polaire (c’est en effet la plus large et la plus 
proche des sources d’air froid), nous arrivons 
au résultat important suivant: Dans le cas d’un 
jet simple, la longueur d’onde maximum des 
cyclones de front polaire est normalement inde- 
pendante de la longueur de la zone déformée et de 
l’ordre de 2.500 km: or c’est bien là, en fait, la 
limite supérieure approximative constatée pour 
la longueur d’onde de ces cyclones. Toutefois 
cela n'exclut pas la possibilité d’une excitation 
de la zone Z’, et par suite celle d’ondes de 
cyclone plus longues et à caractère ondularoire 
moins marqué: ce pourrait être le cas, par 
exemple, pour les perturbations méditerranéen 
nes. 

Voici en outre d’autres propriétés des ondes 
de cyclone théoriques résultant de l’examen 
du tableau I complété par celui du profil A: 

1°) Les ondes sont d'autant plus longues que 
la déformation du jet-stream a lieu plus loin 
de l’axe. 
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2°) La vitesse de propagation d’une onde, 
voisine de la vitesse du vent sur l’axe secondaire, 
est toujours dirigée vers l’Est, et d’autant plus 
grande que le jet est plus rapide et que l’onde 
est plus courte. Ce résultat est en bon accord 
avec l’observation. 

3°) Pour un jet-stream de vitesse moyenne 
(U, de l’ordre de 40 m/sec, cas du profil A), 
cette vitesse de propagation varie, suivant la 
longueur d’onde, de 20 à 32 m/sec si l’excita- 
tion a lieu dans la zone Z, et de 20 à 37 m/sec 
si elle a lieu dans Z’. Cette vitesse est un peu 
forte comparée à celle que l’on observe, et 
qui ne dépasse guère 20 m/sec; toutefois le 
désaccord n’est pas très important, et au surplus 
il n’a rien de surprenant, vu que l’on a totale- 
ment négligé l’action dynamique des couches 
d'air, plus lentes que la jet-stream, situées au- 
dessus et au-dessous de ce dernier. 


2. Ondes de cyclone d’un jet double 


Nous avons vu que la distance entre les axes 
d’un jet-stream double était normalement de 
l'ordre de 15 degrés de latitude au moins, ce 
qui correspond à une pseudo-distance d > 0,26; 
d'autre part la valeur 4.000 km que nous avons 
admise comme longueur d’onde maximum des 
ondes de cyclone correspond à un nombre 
donde n> 5 si on suppose la latitude in- 
ferieure à 60°. Il en résulte que l’on a e~ 24 < 
1/13, valeur petite comparée à l’unité, ce qui 
permet de conclure que dans le problème des 
ondes de cyclone on peut pratiquement remplacer 
un jet double par ses deux jets simples composants. 
Il suffit d’ailleurs pour s’en convaincre de 
considérer par exemple les deux jets doubles 
particuliers ayant pour profils les courbes a et b 
de la figure 1, et pour lesquels la distance 
entre les axes est environ 18 degrés de latitude: 
on voit clairement que les spectres de pro- 
pagation correspondants (fig. 2a et 2b) sont 
confondus avec leurs composants D, et D; 
pour toutes les valeurs de n supérieures à 5 
(A << 0,2). 

Nous voyons ainsi qu’un jet double com- 
porte normalement quatre zones sensibles, 
formant deux systèmes rattachés chacun à un 
axe, et tels que si l’un d’eux est excité il peut 
apparaître des ondes de cyclone ayant leur 
intensité maximum au voisinage de l’axe cor- 
respondant, alors que sur l’autre axe cette 
intensité est au contraire très faible. Toutes les 
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propriétés énoncées pour les ondes de cyclone 
d’un jet simple demeurent valables (à condition 
de remplacer le terme «axe principal» par «axe 
principal le plus proche»); en outre, vu que la 
distance entre les deux systèmes de zones 
sensibles est relativement grande, l'éventualité 
d’une excitation simultanée et indépendante de 
ces deux systèmes a de fortes chances de se 
présenter, d’où la propriété additionnelle 
suivante : 
Dans le cas d’un jet double, on peut avoir simul- 
tanément deux familles distinctes d’ondes de 
cyclone à des latitudes différentes. Or c'est un 
fait bien connu que deux courants de pertur- 
bations peuvent exister en même temps dans 
la zone des latitudes moyennes, et il serait 
important de vérifier que cela se produit lorsque 
le jet-stream possède une structure double. 
Remarque. — On notera que nous n'avons 
fait intervenir nulle part, dans cette théorie 
des ondes de cyclone, la condition que le 
nombre d’ondes n doit être entier. En effet, le 
fait que la zone de déformation du jet-stream 
a une longueur finie ôte toute signification 
physique à cette condition, puisqu'elle corre- 
spond à l'hypothèse d’une répartition zonale; 
il n’y a d’ailleurs pas lieu d’introduire d’autre 
condition de résonance analogue, car la zone 
de déformation va certaimenent en s’estompant 
progressivement à ses deux extrémités. 


3. Cyclones tropicaux et vagues d’Est 


Nous n'avons considéré en somme, jusqu'ici, 
que les ondes de cyclone des fronts polaires, 
mais la théorie précédente peut aussi s’appliquer 
aux perturbations des basses latitudes. A cet 
effet, reprenons le cas particulier du jet simple 
de profil A. En plus de la zone sensible Z’ déjà 
considérée, le tableau I met en évidence 
lexistence d’une autre zone contigüe Z”, 
beaucoup plus large (elle s’étend vers l'équateur 
jusqu’à la latitude 20°), et caractérisée par le 
fait que l’on a à la fois Sy > 10% sec-1 et 
Ly) > 4.000 km, en d’autres termes les pertur- 
bations instables dues 4 l’excitation de cette 
zone ont une instabilité de valeur significative, 
mais leur longueur d’onde théorique est trop 
grande pour des ondes de cyclone, donc 
d'après ce que nous avons dit plus haut elles 
peuvent représenter des perturbations cycloni- 
ques sans caractère ondulatoire, du type des 
cyclones tropicaux. Toutefois nous remarquons 
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que la zone en question se trouve dans une 
region normalement occupée par la ceinture 
des anticyclones subtropicaux, où des mouve- 
ments convectifs 4 grand développement 
vertical ont par suite très peu de chances de 
se produire, et en conséquence la formation 
des perturbations dont nous venons de parler 
est trés peu probable. 

Mais supposons maintenant que le jet-stream 
soit, pour une raison ou une autre, soumis à 
un déplacement général vers l'équateur, sans 
changement appréciable de son profil méridien. 
La zone Z", suivie dans son mouvement par 
la zone Z’, pourra alors atteindre par exemple 
une région de moussons, où les mouvements 
convectifs se développent avec grande facilité, 
donc elle pourra être excitée avec une intensité 
suffisante pour que des perturbations instables 
se développent, et ainsi se trouverait expliquée 
la naissance des cyclones tropicaux. En même 
temps la zone Z’, se trouvant à proximité de la 
même région d’instabilité thermique, peut elle 


D . 
aussi devenir le siège de perturbations instables 


ayant cette fois le are d’ondes de cyclone: 
ce pourrait être là l’explication de l’origine de 
ce que plusieurs auteurs ont décrit sous le nom 
de «vagues d’Est ou «easterly waves», et qui 
apparaissent précisément en bordure septen- 
trionale (dans l'hémisphère Nord) des régions 
de moussons, par exemple en Afrique Occi- 
dentale Francaise. 


4. Rôle possible de ozone 


Nous n'avons considéré jusqu’ici que l’in- 
stabilité thermique de la basse atmosphére 
comme cause initiale de la cyclogénése, mais il 
est évident que tout autre phénoméne suscep- 
tible de ralentir localement le jet-stream peut 
jouer le même rôle. Ainsi la déformation des 
lignes de courant associée aux centres d’action 
semi-permanents peut fort bien produire un 
tel résultat si son amplitude est suffisante, et 
dans ce cas on aurait affaire à une cyclogénèse 
purement dynamique; il semble bien, en effet, 
que certaines perturbations cycloniques ne 
peuvent pas s'expliquer autrement. Mais il 
est également possible que la couche d’ozone, 
qui est probablement le siège de mouvements 
convectifs importants par suite de son fort 
pouvoir absorbant dans l’ultra-violet, inter- 
vienne dans certains cas au même titre que la 
troposphère inférieure, par exemple si un 
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maximum de concentration en ozone se trouve 
créé le long d'un certain parallèle comme 
conséquence de la circulation générale: cela 


permettrait peut-être d’éclaircir le problème des 
relations entre l’activité cyclonique et l’activité 
solaire. 
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Die Konfiguration mittlerer Höhenströmungsfelder 


und ihre Ursachen 


Von F. WIPPERMANN, Mainz 


(Manuscript received 1 Feb. 1952) 


Abstract 


Some computations are presented in order to estimate the relative importance of dif- 
ferent factors for the generation of the quasi-stationary perturbations in the westerlies, such 
as the effect of the large mountain barriers, the heat exchange with the surface of the 
earth and surface friction. For that purpose the vorticity advection of the mean field of 
motion has been computed. It is concluded that the three factors mentioned seem to be 


of the same order of magnitude. 


Selbst über Jahreszeiten gemittelte Karten 
der Höhenströmung zeigen bekanntlich noch 
ausgeprägte Abweichungen von einer reinen 
Zonalzirkulation. So sind die betreffenden 
Winterkarten charakterisiert durch zwei mar- 
kante Hochdruckrücken an der amerikanischen 
und europäischen Westküste, sowie zwei ent- 
sprechende grosse Tröge an der asiatischen und 
amerikanischen Ostküste. Ein dritter, weniger 
stark ausgeprägter Trog findet sich über dem 
europäischen Russland. Die Ursachen dieser 
Störungen können keineswegs als geklärt an- 
gesehen werden, vielmehr sind die einzelnen 
Meinungen zum Teil sehr widersprechend. 
So versucht man auf der einen Seite in einer 
orographisch-dynamischen Betrachtungsweise 
die Existenz dieser Störungen einer Wir- 
kung der grossen Kettengebirge zuzuschreiben, 
die in ihrem Lee grossräumige horizontale 
Wellen auslösen. Errechnete Ergebnisse lassen 
sich mit den Beobachtungen in Einklang brin- 
gen. (CHARNEY, ELIASSEN 1949; BOLIN, 1951.) 
Demgegeniiber glaubt man auf der anderen 
Seite die Ursache vornehmlich im thermischen 


Einfluss suchen zu miissen, der sich (wie hier 
im stationären oder doch quasistationären Zu- 
stand) in der Baroklinität des Feldes ausweist. 
Wegen der mathematisch weitaus schwieri- 
geren Behandlung des Problemes im barokli- 
nen Fall ist man genötigt sich bei der soge- 
nannten thermisch-synoptischen Betrachtungs- 
weise auf qualitative Erklärungen (SUTCLIFEE, 
1951) zu beschränken. Man vermag jedoch 
vom beiden Seiten Argumente für und gegen 
vorzubringen, die nicht übersehen werden 
dürfen. 

Es erscheint aus diesem Grunde lohnend, die 
Ursachen für die quasistationären Störungen 
der Zonalzirkulation zu untersuchen und zu 
prüfen, in welchem Grössenverhältnis die ein- 
zelnen Einflüsse zu einander stehen. 


Zu diesem Zwecke leiten wir uns eine Glei- 
chung her, die die horizontale Advektion der 
auf den isobaren Flächen senkrecht stehenden 
Vorticity-Komponente (als Charakteristikum 
für die Anordnung des Strömungsfeldes) auf 
der einen Seite mit den verschiedenen Ein- 
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flüssen auf der andern Seite verbindet. Im x,, 
Y» p-System, einem rechtwinkligen Koordi- 
natensystem, in welchem die x,, y,-Ebene 
mit den Flichen p = const identisch ist, 
lautet die Vorticity-Gleichung: 


d 
EP C+N=--—-( + f) div, v— 

x ie V il 
—k [Vim x 5 |] +R (1) 


wobei w, die Geschwindigkeit in Richtung der 
p-Koordinate ist, v den Geschwindigkeitsvek- 
tor auf den isobaren Flächen bedeutet, & die 
relative Vorticity, f den Coriolisparameter und 
R den Reibungsterm symbolisiert. Unter Ver- 
nachlassigung des 2. Termes der rechten Seite 
als klein gegenüber den andern, sowie für € < f 
folgt bei Verwendung der Kontinuitätsglei- 


chung 


Po Po P, 
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Hieraus ergibt sich nach Ausführung der Inte- 
gration durch Mittelbildung über den Bereich 
pP =Ppo bis p =o der äquivalent-barotrope 
Betrag fiir die individuelle Anderung der auf 
den isobaren Flächen senkrecht stehenden 
absoluten Vorticity-Komponente: 


FR Po 
RE+P=L m+ ET Rd 6) 


Wegen der Anwendung auf zeitlich gemittelte 
Karten können wir mit quasistationären Ver- 
hältnissen rechnen. Für solche erhält unter 
Überführung der Geschwindigkeit in Richtung 
der p-Koordinate in die (metrische )Vertikal- 
geschwindigkeit, sowie unter Verwendung 
des geostrophischen Windes v, an Stelle des 
Windvektors auf isobaren Flichen v, also mit 


ae, 


Gl. (3) die Form 
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Hierin bedeuten J, den Jacobi’schen und A, 
den Laplace’schen Operator, jeweils auf Flä- 
chen p = const. geltend, h die Höhe der ent- 
sprechenden Druckflache im metrischen Sy- 
stem x, y, z, und H die Höhe der homogenen 
Atmosphäre. 

Im G. (6) gibt der erste Term der rechten 
Seite den Einfluss der Topographie der Erd- 
oberfläche wieder, da eine Vertikalgeschwindig- 
keit im Niveau z =o nur dann nicht ver- 
schwindet, wenn die Erdoberfläche als Be- 
erenzunesflache geneigt ist. Nach dem Ansatz 
von CHARNEY u. ELIASSEN (I) können wir für 
diesen Term auch schreiben 


= £k[Vzz x Vol (7) 


0 


worin 2g (x, y) die Höhe der Erdoberfläche 
iiber dem Niveau z = 0 angibt. Der 2. Term 
der rechten Seite enthält die Reibungseinfliisse ; 
die Kompliziertheit der verschiedenen mög- 
lichen Ansätze gestatten keine quantitative 
Erfassung. Um jedoch eine spätere überschlags- 
mässige Abschätzung der Grössenordnung zu 
ermöglichen, sei der sich durch seine Einfach- 
heit auszeichnende, jedoch nur in der Boden- 
reibungsschicht geltende Ansatz nach Guldberg- 
Mohn verwandt. Es wird dann in Gl. (1) 


R=—uö+kl[vx Von] (8) 
worin # die Guldberg-Mohn’sche Reibungs- 


konstante bedeutet. Entsprechend heisst der 2. 
Term dann in Gl. (6): 
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Der 3. Term schliesslich stellt den iiber die 
gesamte Atmosphäre gemittelten Vorticity- 
transport senkrecht durch die isobaren Flächen 
hindurch dar. Er ist mit demjenigen im äquiva- 
lent-barotropen Niveau nur dann identisch, 
wenn letzteres definiert ist unter der Voraus- 
setzung 


I) =al/ rious lr, top 


barotrop. 


aequiv. 
(10) 


und nicht etwa unter der iiblichen Definition 
fiir ‘dieses Niveau: 


d 


SEIN=S CTP (a) 


Da er vom Mass der Baroklinität abhängt, gibt 
dieser Term den Einfluss der Warmefunktion 
wieder, d. h. all der nichtadiabatischen Effekte 
der Wärmeaufnahme oder des Entzuges in- 
folge Strahlung, Kondensation, Verdunstung, 
Austausch von der Unterlage her usw. Fiir die 
Richtigkeit dieses Schlusses mag angeführt 
werden, dass sich im nichtstationären Fall dieser 
Term rasch zu dem Korrekturglied der zeitlich 
lokalen Änderung umformen lässt, das CHAR- 
NEY, FJORTOFT u. VON NEUMANN bei Ver- 
wendung des vereinfachten baroklinen Mo- 
delles gebrauchen. Der Zusammenhang zwi- 
schen Baroklinitäit und Wärmefunktion ist 
andererseits im stationären Falle klar. 


Es soll zunächst der Versuch unternommen 
werden diese 3 Einflüsse in ihrer Grössenord- 
nung gegeneinander abzuschätzen. 

In Gl. (7) erhalten wir mit einem Bodenwind 
von 10 m/s und einer Neigung der Erdober- 
fläche in Richtung dieses Windes von 2 500 
m : 1000 km (maximale Neigung in einem 
ausgeglichenen Profil) für die Grösse des 
ersten Termes in Gl. (6) den Wert von 
u à 10 Ce 

Der Reibungsterm ist einerseits vom Strö- 
mungsfeld selbst abhängig, (in stark verwirbel- 
ten Feldern — ohne Unterschied ob zyklonal 
oder antizyklonal — ist die Reibung grösser) 
zum andern aber auch vom Untergrund, näm- 
lich der räumlichen Änderung des Reibungs- 
koeffizienten. Der erste Anteil lässt erkennen, 
dass Hoch- u. Tiefdruckgebiete zerstört werden, 
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der zweite Anteil wird besonders gross an den 
Kiisten erwartet werden diirfen. Die Mohn- 
sche Konstante ist von STÜVE (1931) mit 
2—8 : 10 5 sec! angegeben; die Werte von 
Baur u. Parties (HANN-SÜRING, 1938) für den 
Reibungsansatz von Hesselberg und Sverdrup 
lassen sich umrechnen auf 5 - 10°? für Meer und 
9: 10° für Land. Um die ausschliessliche Gül- 
tigkeit in der Bodenreibungsschicht zu berück- 
sichtigen, sei angenommen, dass die Reibungs- 
konstante bis zur Reibungshöhe p = 850 mb 
auf den Wert o nach einem irgendwie anzu- 
setzenden, hier linearen Gesetz abklinge. Den 
für diese Schicht berechneten Einfiusswert 
denken wir uns dann über die ganze Atmo- 
sphäre verteilt wirksam. Auf diese Weise er- 
gibt sich mit u = 8. 10° und ¢ = 3 - 10 für 
den ersten Reibungsterm 1—2 10719 sec~?. 
Zur Berechnung des 2. Reibungstermes mögen 
die oben genannten für Meer und Land ver- 
schiedenen Werte eine Verwendung finden, 
indem wir den Übergang auf einem ca 300 km 
breiten Küstenstreifen annehmen. Es ist also 
Vu='1,3%109% an sect. Vertshrensee 
dann wie beim ersten Reibungsterm so folgt 
für den zweiten der Wert ı - 10-10 sec2. 

Die Abschätzung des Einflusses der Wärme- 
funktion in Form des Vorticitytransportes 
senkrecht zu den isobaren Flächen ist nicht 
minder schwierig. Die Verteilungen & = £ (2) 
zeigen, dass das Maximum in der Regel an der 
Tropopause (p= 225 mb) liegt, während sich 
bei p = 700 mb & meist dem Werte Null 
nähert. Mit einem Gefälle von 3: 10-5 auf 300 
mb und einer in die Druckkoordinate umge- 
rechneten Vertikalgeschwindigheit von s cm/s 
(in gemittelten Druckfeldern) folgt für diesen 
Term 2—3 : 10710 sect. 

Es berechtigen die vorgenommenen Ab- 
schätzungen infolge ihrer Grobheit zu keinerlei 
Zuerkenntnis des Dominierens für einen der 
Einflüsse. Es ist jedoch die Aussage berechtigt, 
dass von den drei Einflüssen keiner gegenüber 
den andern an allen Stellen des Feldes vernach- 
lässigbar klein ist. 


Eine andere Betrachtungsweise der Verhält- 
nisse ergibt sich aus Gl. (6). Gelingt es aus einer 
zeitlich gemittelten Karte die Grösse der lin- 
ken Seite zu entnehmen und synoptisch darzu- 
stellen, so ist damit eine räumliche Verteilung 
der Summe der auf der rechten Seite stehenden 
Einflüsse gegeben. Aus dieser wird man unter 


DIE KONFIGURATION MITTLERER HOHENSTROMUNGSFELDER 


Berücksichtigung der geographischen und 
orographischen Verhältnisse Schlüsse auf die 
Existenz des einen oder anderen der Einflüsse 
ziehen können. Die mangelnde Möglichkeit 
die genannte Grösse durch Mittelbildung über 
Druckkoordinate für eine äquivalent-barotrope 
Atmosphäre zu ermitteln, nötigt zur Anwen- 
dung auf die vorliegende Karte des soo mb- 
Niveau unter der Voraussetzung, dass dieses 
nach den Überlegungen von CHARNEY u. 
ELtassEN (1949) dem äquivalent-barotropen 
Niveau entspricht. Hiergegen lässt sich der 
nicht unberechtigte Einwand erheben, dass die 
500 mb-Fläche nur annähernd dem äquivalent- 
barotropen Niveau entspricht, dass sich viel- 
mehr gerade in der Differenz 
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die Wärmefunktion wiederspiegele und sich 
deren räumliche Verteilung wahrscheinlich 
darstellen lasse durch eine Topographie des 
an einer räumlich gemittelten Vertikalverteil- 
ung von J, festgestellten äquivalent-barotropen 
Niveaus. Noch schwerwiegender wird dieser 
Einwand bei der Feststellung, dass die Fest- 
legung der soo mb-Fläche als äquivalent-baro- 
tropes Niveau nicht unter der Bedingung (10) 
sondern (11) erfolgt. Die vorhandenen mitt- 
leren Höhenkarten schliessen aber alle anderen 
Möglichkeiten als diejenige der Bestimmung 
von J, unter der Annahme 
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aus. Das Ergebnis scheint diese, notgedrungen 
alleine in Frage kommende Verfahrensweise 
trotz aller möglichen Einwände zu rechtfer- 
tigen. 

Als Grundlage diente die von FLOHN (1949) 
entworfene Mittelkarte der absoluten Topo- 
graphie der soo mb-Fläche für den Winter im 
Bereich von 120° W—s50° E und 35° N—70°N. 
Ermittelt wurde also die geostrophische Advek- 
tion der geostrophischen Vorticity, wozu 
mithin eine dreifache Differentiation der topo- 
graphischen Werte nach den Ortskoordinaten 
erforderlich war. Bei der verhältnismässig 
geringen Genauigkeit der Ausgangskarte muss 
dieses Moment sofort bedenklich stimmen, 
doch scheint das Ergebnis den Versuch zu 
rechtfertigen. Um eine möglichst exakte 
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Berechnung zu gewährleisten, wurde das Feld 
zunächst in einzelne ausgeglichene Profile 
(meridional von 5 zu 5 Grad Länge, zonal von 
2,5 zu 2,5 Grad Breite) zerlegt. Jedoch wurde 
nach jedem Differentationsgang wiederum 
eine kartenmässige Darstellung angefertigt, um 
eine synoptische Ausgleichung zu ermöglichen. 
Erst den resultierenden Karten wurden dann 
wieder die neuen Profile entnommen. Durch 
das dreifache Differenzieren einerseits und die 
Behandlung in zwei Koordinaten zum andern 
bleibt das Ergebnis bis zum Schluss derart un- 
übersichtlich, dass es unmöglich ist, bei den 
notwendigen Ausgleichungen durch geschickte 
Durchführung derselben irgendwelche gewoll- 
ten Effekte zu erhalten; vielmehr erscheint 
hierdurch eine gewisse Objektivität bei der 
Bearbeitung gesichert. : 

Das Ergebnis dieses Versuches ist recht be- 
friedigend. Das resultierende Vorticity-Trans- 
portfeld ist in seiner Anordnung durchaus 
sinnvoll, ja es lässt sogar eine recht gute Uber- 
einstimmung mit dem Feld von Einzelein- 
fliissen in den Gebieten erkennen, in denen 
diese dominierend vorhanden sind. Die Dis- 
kussion des Feldes sei mit einem Vergleich der 
Karte von MÖLLER (1950) verbunden, die die 
räumliche Verteilung der Wärmeaufnahme 
und des Entzuges in den untersten 3 km der 
Atmosphäre im Winter darstellt. MöLıer hat 
diese Verteilung aus der mittleren Advektion 
der Temperatur sowie den Austauschwärme- 
strömen erschlossen, also auf gänzlich andere 
Weise als hier. Sieht man von den verschie- 
denen Geltungsbereichen der beiden Karten 
in der Vertikalen ab, so stellt die MOLLER’sche 
Karte den dritten Term der rechten Seite in 
Gl. (6) dar. Sie müsste also mit der unsern über 
dem Atlantik in bester Übereinstimmung sein, 
da hier die von uns noch mitberücksichtigten 
Gebirgseinflüsse ganz verschwinden, die Rei- 
bungseinflüsse aber minimal und in Bezug auf 
die Mohn’sche Konstante überall gleich sind. 
Die Übereinstimmung ist in grossen Zügen 
hier auch recht gut, Auf beiden Karten zieht 
sich über den Atlantik von SW nach NE ein 
Gebiet mit Wärmeaufnahme — etwa dem 
Golfstrom folgend. Der Unterschied an der 
skandinavischen Küste darf vielleicht auf das 
Wirksamwerden der skandinavischen Gebirge 
gewertet werden, der im Luv negative, im Lee 
positive Werte zeitigen muss. Letzterer würde 
dann mindestens auch zum Teil das positive 
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Abb. ı a: Mittlere Wärmeaufnahme und -entzug 

zwischen o und 3 km Höhe im Januar (in je o,1 cal 

pro cm? min 3 km) erschlossen aus der mittleren 

Advektion und den Austauschwärmeströmen nach 
MÖLLER. 


Areal über Ostsee, Baltikum und Finnland 
erklären können. Ein zweiter grösserer Un- 
terschied zeigt sich vor der amerikanischen 
Ostküste südlich von 45° N, wo in einem 
Gebiet ausgesprochener Wärmeaufnahme das 
isobare Transportfeld negative Werte auf- 
weist. Letzteres lässt sich wenigstens qualitativ 
durch das Wirksamwerden des 2. Reibungs- 
termes erklären, wofür auch die Anpassung 
dieser Abweichung an die Küste spricht. Wie 
in der Möurer’schen Karte finden wir ebenfalls 
ein Gebiet mit Erwärmung über der Davis- 
strasse, hier noch ausgedehnt bis nach Neu- 
fundland und über die Hudsonbai. Die nega- 
tiven Werte im Luv und die positiven im Lee 
des amerikanischen Randgebirges werden er- 
wartet, erweitert durch die mit der MÖLLER- 
schen Darstellung übereinstimmenden, von W 
nach E in den Kontinent hineinreichenden 
Zunge positiver Werte. Gerade die Überein- 
stimmung des letztgennanten Gebietes in beiden 
Karten zwingt jedoch nicht dazu, diesen Ein- 
fluss notwendig auf eine Erwärmung zurück- 
zuführen. Es sei nämlich vermerkt, dass die 
Mörzer’schen Berechnungen auf der not- 
wendigen Vereinfachung horizontaler Advek- 
tion beruhen und deshalb auch noch dyna- 
mische Eftekte enthalten, die aus der Vernach- 
lässigung der Vertikalbewegungen resultieren. 
Dies wird vorallem auch für das starke Abküh- 
lungsgebiet im Nordwesten Amerikas gelten. 


Abb. ıb: Wärmeaufnahme und -entzug, Reibungs- 

effekte und Einfluss der Gebirge dargestellt durch das 

Feld der Advektion absoluter Vorticity im “äquivalent- 

homogenen’ Niveau (nahe soo mb) in je 2,5 x 10710 
cm/sec? im Winter. 


Auch die Übereinstimmung beider Karten in 
den mediterranen Ländern ist befriedigend, 
selbst die Leewirkung der Alpen scheint sich 
bemerkbar zu machen. — Auf dem Atlantik 
selbst kann nur die Wärmefunktion und der 
Reibungseinfluss wirksam werden. Die gute 
Übereinstimmung beider Karten in diesem 
Raume lässt einen dominierende Einfluss 
der Wärmefunktion in diesem Gebiet ver- 
muten. 


Den vorangehenden Betrachtungen glaubt 
der Verfasser insofern einen Wert beimessen zu 
dürfen, als sie erkennen lassen, dass die das 
quasistationäre Höhenfeld gestaltenden Ein- 
flüsse im allgemeinen von der gleichen Grössen- 
ordnung sind, und dass eine dominierende 
Wirksamkeit einzelner Effekte auf bestimmte 
Areale beschränkt bleibt. Bewerkenswert ist 
weiterhin, dass der Abstand der Haupteinfluss- 
gebiete gleichen Vorzeichens in der Karte des 
geostrophischen Transportes der vertikalen 
Komponente absoluter Vorticity im Bereich 
der Wellenlänge der freien planetarischen Wel- 
len (RossBy, 1939) liegt, und dass somit der 
Faktor der Resonanz auch bei Berücksich- 
tigung des Einflusses gänzlich verschiedener 
Effekte seine Bedeutung behält. 

Herrn Dipl. Met. H. Zimmermann bin ich 
für die Unterstützung bei der Berechnung der 
Vorticity-Transportkarte zu Dank verpflichtet. 
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Meridional Flux of Angular Momentum in the Tropics 


By V. P. STARR, Massachusetts Institute of Technology 
and R. M. WHITE, U.S.A.F. Cambridge Research Center, Mass. 


(Manuscript received 26 March 1952) 


Abstract 


The results of an attempt to measure the meridional flux of angular momentum across the 


vicinity of 13° N latitude from actual wind data are presented. It is found that the average 
flux for one year is directed northward but is much smaller than the corresponding flux at 
31° N as obtained from a similar hemispheric study made earlier for that latitude. According 
to the data the flux at 13° is due mainly to horizontal exchange processes which are most vigorous 


at an elevation of about 35,000 feet. 


I. Introduction 


During the course of an extended program 
of research concerning the angular momentum 
balance of the northern hemisphere, the 
writers have reported elsewhere (STARR and 
WHITE 1951, 1952) the results of flux com- 
putations from actual-wind data in the sub- 
tropics. In the study described in these two 
previous papers, henceforth denoted as (A) 
and (B) respectively, the meridional transports 
were evaluated from a string of upper-wind 
stations encircling the earth in the vicinity of 
31° N latitude for a period of one year. In the 
region of the subtropics and at more northerly 
latitudes the character of the angular momen- 
tum balance has also received observational 
elucidation from the geostrophic-wind studies 
of WipGER (1949) Mintz (1951), WHITE and 
CooLEY (1952) and Lorenz (1952), and like- 
wise from the actual-wind studies of PRIESTLEY 
(1949, 1951a, b), STARR (1950) and NYBERG and 
SCHMACKE (1951). 

However, much interest is centered in regard 
to the character of the transport mechanisms 
(vertical and meridional) which may be 
operative in the more tropical regions where 


geostrophic wind measurements are of doubtful 
reliability and are more difficult to obtain. For 
these reasons, and also because of other con- 
siderations, recourse must be made to prc- 
cedures involving the use of actual-wind data 
in order to throw light upon the problem in 
this zone. In the present paper an account is 
given of an attempt to duplicate the studies 
contained in (A) and (B) for a string of upper- 
wind stations in the vicinity of 13° N, to the 
extent that this is possible with existing data. 
Although an endeavor was made in (B) to 
show the seasonal variation of the flux, no 
corresponding attempt was made for 13° N, 
due to the fact that a somewhat smaller total 
amount of data was available in the tropics 
and also for other reasons. 

From general physical considerations one is 
led to certain over-all expectations concerning 
the total meridional flux across the latitude 
circle in the vicinity of 13° N. The more 
important of these are the following. (a) Since 
the storage capacity of the atmosphere for 
angular momentum is small when a period of 
the length of one year is considered. the total 
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Table 1. List of key stations (numbered) and alternate stations (italics) 


SS ee ee eee ee 


Station Latitude (N.) Longitude Altitude (ft) Type 

Tel BETS ay, Seats Les ae ae ea ae 702 06 COS OS E 899 pilot balloon 
NET no CO re 31 O2 20.18 HSS) » 
MIA SIN ween hls 13 48 09 oo E 1604 » 
NAO DE ER Le res e 14 30 04 Tr) OW 906 » 

ID AMAR ENS ER ee con 12 38 08 or W 1086 » 

Dep Piethashensne MS MAR re 13 37 25 DOME 2395 » 

DCR DO PRES TS OS cee TNT IS 28 36 24 E 1644 » 
JEL MONAT SSS a nee 13 10 30 14 E 1887 » 
NOS RES SSSR ard Se ges 13 Io SAONE 1253 » 
MialaRal ER sea cee ee 09 33 31 39 E 1276 » 
ae CS PRO Adler. 07 42 28 org B 1440 » 

NPA Khormaksar) nen... 12 so 45 Oni 13 3 radio wind 
Aden Scheikh Othman ...... 12 54 44 58 E 33 pilot balloon 
Royale 2 satis SAW kx sts 140530 io NE 46 » 
Komarani Island 4e 2" 2.0 IS 20 42 a7 18 20 » 

mM IIG NM ODOMrcieyschors te ae, as os Io 49 Visa Ae 5 256 » 
TO MES SRE D IS 55 73 40 E — » 
SULA AD TOR Ed espe? ele 14 41 77 SRE 1148 » 
IMAGES Re kN ea are ave 1 04 80 SRE 52 » 
Mangalorezee 2... June ES T'AS NE) 72 » 

TONER GOERIA Eten hex, cs 09 58 76 14 E Io » 
DEST OUT Se Srnec oe 08 29 76 Sy7 18 200 » 
NEN AMEN 08 PT Gg} XS 1S) 10192 » 

BESASON, Reese ce usb os IO 49 106 39 E 33 » 

Mingaladopnre a ee loan ehe, 16 56 060" 07 EB 92 » 
BE AVOVE en oi Sis ches in: her 14 05 98 Ta) 112 » 
Kretonas Bonn. as ce 09 59 98 ag JE 122 » 
IGT 68 DC ee Seb Oo BOS oo Be 05 18 100.0. 10WE 13 » 
INCU CT TS Sie OS Cee 16 02 108 ign dz, 16 » 
FAT SERRE SEE CORRE 16 23 III ay) 3 20 » 
DEN ee OMS Eee 16 40 105$ oo E 502 » 
Phnom Beni ea here II 33 104 51 E 33 » 

Gam Clarke Hel CRETE 15 Io 120 34 E 644 radio wind 
TOGO ee ae ae TS it 12 One ae 12 pilot balloon 

Th Na ae EL C0 WO) 138 08 E o radio wind 

Sy Hlarmongkielder sanken Ties ii AAO 176 » 

OL RANBIR Vor 08 43 TOMATE Io » 
AMEL ES RTE ARE 0 19 IS 166 30 E o » 

TEs OAR mo On oO RON Soto 16 44 169 31 W 20 » 

EEE PR eee deu à à TON 155 04 W ag » 

25. NÉGOCE EE EE 16 20 95 14 W 75 pilot balloon 
TG CTU Zo = fabs to ole delete 19 2} CE A Io » 
Ciudad del Carmer ........ 18 39 OI so W 5 » 

RON ETC OU E Er denen 12) 08 SCHEN 178 radio wind 

13. Plato Magdalena #72. 09 48 74 48 W 190 » 
PAROI AN AE ES EAS Er 08 58 79 2 N 21 » 
ISM: CORPS tate Gann Oc 10 55 74 45 W 39 » 

ie Waller Bielden ne es 10 36 61 12 W 0 » 

SI TOS Nore dt cine 18 28 66 07 W 82 » 

TS BESTER balay Ata) Aco cach CREER 17 nr Ke) 17020 Wi) 13T » 

Sis PIS 0 mo Lo on Re ete Toy 204 Gy ye)» AYE 10 pilot balloon 
IDE NEP occ cod ogoe gor T4) Ad 17 25 W 62 » 
Gai). Ree ee a aE ees 16 44 22 57 W 180 » 


amount transported into the polar cap north winds near the surface, the northward flux at 
of a given latitude should be equal to the 13° N should be smaller (algebraically) than 
amount removed by surface torques acting on the flux at 31° N. (c) Whether or not the 
the cap. (b) Since the zone between 13° and flux across 13° N is actually negative (directed 
31° N is a region of predominant easterly southward) cannot be determined from general 
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considerations except as noted below, since 
one may not assume complete symmetry of 
the two hemispheres about the equator in this 
connection. (d) Since at 31° N angular momen- 
tum is transported strongly northward by the 
very large disturbances in the upper westerlies, 
one might perhaps expect that this effect would 
not be suppressed completely in a distance so 
short as 18° to the south. 

The degree to which such anticipations as 
those listed above are not in conflict with the 
numerical findings obtained, constitutes one 
rather general criterion for the plausibility of 
the results. 


2. Observational material 


As in the study reported in A) and (B) a 
string of key stations was used. These 15 
stations are listed in table 1. Since reports were 
often missing from the key stations, alternative 
stations in the vicinity of the key stations were 
used in order to increase the amount of obser- 
vational material. The alternative stations are 
entered in italics below the key stations in the 
table. Although it would have been advanta- 
geous to use the same period of time as was 
used for the study at 31° N, various factors 
made this impractical. The average latitude of 
the key stations is 13° N. 

Wind reports for the levels 2, 6, 10, 14, 20, 
25, 30, 35, 40, 45, 50, and 55 thousand feet 
were taken from the data tabulations of the 
Daily Series Synoptic Weather Maps prepared 
by the U.S. Weather Bureau, in cooperation 
with the Army, Navy and Air Force for the 
period from 1 July 1949 to 30 June 1950 for 
the hour 0300 GMT for each day. For the 6, 
10, 20, 30, 40, and 55 thousand-foot levels 
supplementary data as reported for the 850, 
700, 500, 300, 200, and 100 millibar levels, 
respectively, were added, although preference 
was given to constant-height data where 
available. No attempt was made to correct 
original data, although in a few cases garbled 
reports were omitted. The frequency of 
available observations at the various levels is 
given in table 2. This table is analogous to table 
5 in (B) in that it gives the percentages of the 
total possible observations, i.e., 365. Each 
wind observation was resolved into the east- 
ward component # and northward component 
v, tabulated in m sec-t. 
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3. Computation of transports 


In order to avoid repetition, only a brief 
statement is made here of the methods used to 
measure the angular-momentum flux, the 
reader desiring more details being referred to 
(A) and (B). The notation used may be 
summarized as follows: 


Table 2. Percentage of total possible observa- 
tions at each level for each station 


Altitude in thousands of feet 


Station SS 

22100 (10/14/20 25| 30) 35|40\45| 50/55 
Gao yo ee 63165159/46|15| I| O} 0) O} o| o| © 
El Fasher ....| 04lo2|05/87|60/13|11|11| 8] 4] 1| I 
INCH 44154153128147122/45123|38|19|15123 
Trichinopoly ..| 99] 1|99| 0|90)7873| I] ©} oO} o| o 
SAISON ec 100/95|95|75|6ılıı| 6| 3| o| o| OF O 
Clarke econ 84186183175173152144132133123|16|11 
ap ote 30/2s|19|12|10| 4] 2) I] I] I} o| o 


Harmon Field} 77|85|83|68)80|63|73|61|68 50/39/28 


Kwajalein ....| 92!95]92}86]93|79|82]63]70/49]3 8134 
Johnson...... 82188186/74183166/73158164147134128 
Tehuantepec..| 78175 Ela 3| ol ol o| oO} ol o 
Managua 2 © 3 | Du ea) = 9| 0,6) a} 2 TE 
Plato Magda- 

lena Beets 74|84183]70|83|63|65|50 47/28] 9] o 
Waller Field..| 93|98|97|86193|78|88|74[81|65|50|37 
Dakar worden 65166162153155147150/44143135125|20 


u eastward component of the wind. 

v northward component of the wind. 

[x] space average of a quantity « over the 
length of the complete latitude circle. 

& deviation of a quantity & from its space 
average [x]. 

% time average of a quantity &. Time and 
space averaging processes are not neces- 
sarily commutative. 

[x] deviation of the space average [x] from 
the space-time average [a]. 

{x} arithmetic mean of all individual observa- 
tions of & at a given level during entire 
time period considered. 

n number of days for which observations 
were available at a given level. 

N number of individual observations present 
at a given level during entire time-period 
considered. 

r coefficient of linear correlation for the N 
pairs of u and v. 


In (A) it was shown that the quantity [u v] 
which represents the northward transport of 
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linear momentum per unit mass at a given 
level may be resolved into three components, 
namely: 


(i) Transport of the first species, [ul vl. 
This component represents the contribution 
due to a net meridional mass flux [v] during 
the entire period at the level in question. To 


the extent that a non-zero value of |v] at one 
level is usually compensated by values 
having an opposite sign at other levels, this 
component depends upon the existence of so- 
called mean meridional circulations. 


! 


(ii) Transport of the second species, [u]’ [v]’. 
This component represents the contribution 
from fluctuations in [v] due to its possible 
correlation with [u] in time. 


(iii) Transport of the third species, [u’ v’. 
This component represents the contribution 
due to a correlation between # and v along 
the length of the latitude circle. It does not 
require the presence of net air motions [v] 


or [v]. 


Since one basic question involved in the 
subject relates to the role of mean meridional 
circulations, it is useful to form the quantity 
{uv} — {u} {v} which is analogous to the sum 
of the transports of the second and third 
species although the mode of averaging used 
is slightly different as explained in (A). 

The results are given in table 3 which is to 
be compared with table 4 of (B), both being 
of the same basic form. The table gives the 
transports of each of the three species in 
columns 6, 7, and 8 together with various 
other quantities defined in terms of the symbols 
already recapitulated here and explained more 
fully in (A). À number of the quantities carry 
confidence limits defined as twice the standard 
error indicating approximately the 95 per cent 
confidence level. Vertical integrals of several 
quantities with respect to mass (standard 
atmosphere pressure) are given at the foot. 
These figures are to be interpreted as the 
northward transport of linear momentum per 
unit length of the latitude circle for the layer 
from 2 to 55 thousand feet. For comparison 
the corresponding integrals for 31° N as 
reported in (B) are also given at the foot of 
table 3. 


ran 
4. Discussion of numerical results 


The appraisal of the numerical results 
obtained depends upon a simultaneous con- 
sideration of various factors which might limit 
the accuracy of the measurements of the flux 
and other computations. By and large these 
factors are similar to those enumerated in (A) 
and to some extent in (B), although certain 
points are worthy of especial note in view of 
the circumstances of the present study, speci- 
fically the following: 


(a) There exists a rather serious gap in the 
station network extending from Johnson 
Island eastward to Central America. There 
appears to be no source of reports from this 
region which is known to the writers. 


(b) When the twelve values of [v] given 
in column 3 of table 3 are averaged vertically 
with respect to pressure one obtains the net 
meridional air flow for the entire layer. This 
quantity turns out to be — 0.16 m sec-!. Since 
it is quite improbable that a return flow exists 
below 2,000 ft., it would be necessary that an 
average northward flow of rather large 
magnitude should exist in the remainder of 
the atmosphere above 55,000 feet, if this figure 
represents actual conditions. This appears to be 
improbable. The most likely inference to be 
drawn is that it represents evidence of a slight 
lack of randomness of the observations with 
respect to the troughs and ridges of the stream- 
line patterns, which in turn might be related 
to the gap in the observational network 
mentioned in (a) above. In any event the 
figure is rather excessive, the corresponding 
value at 31° N being only + 0.04 m sec! as 
obtained from data in (B). 

(c) The total number of individual wind 
observations is 27,985 at 13° N as compared 
with a total of 38,301 at 31° N. This is a re- 
flection of a poorer station network and more 
missing reports. This and other shortcomings 
such as the one discussed in (b) might be made 
less serious by using more than one year o 
data. 

(d) Resort had to be made to the use of 
relatively more pilot balloon stations than were 
used at 31°. This fact contributes to a paucity 
of observations at high levels and may be a 
factor contributing to a slight lack of random- 
ness discussed in .(b). 


De 
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(e) The range of latitude of the individual 
key stations is only 8 degrees which is much 
less than for the network at 31° N. 

(f) A large number of alternative stations 
was used. This in itself should not be neces- 
sarily detrimental to the results, except to the 
extent that a lack of randomness might be 
introduced in the data. 


(g) The confidence limits for the flux [u v] 


and the component [u’ v’] are smaller at 13° 
than at 31° N, but the values of the quantities 
are also smaller so that on the whole less 
confidence is indicated. On the other hand, 
the alternative method of averaging given in 
columns 10—14 of table 3 gives results which 
compare very well with the first technique. 

Generally speaking the data fulfill plausible 
expectations such as those mentioned in the 
introduction. The following points may be 
noted: 


(t) The vertical integral of the flux [u »| 
may be converted into the total flux of 
angular momentum by multiplying by the 
length of the torque arm and the length of 
the latitude circle. The result is + 8.5 x 105 
gm cm? sec~!, which is considerably smaller 
than the corresponding value at 31° N, i. e., 
+ 31.5 x 10%. The difference represents essen- 
tially the rate at which angular momentum is 
acquired from the surface, since probably only 
small contributions to the fluxes are made by 
the portions of the atmosphere not considered. 

(2) The largest contribution to the angular 
momentum flow across 13° N is due to the 
transport of the third species, the value being 
+ 5.2 x 10% gm cm? sec-1. The next largest 
is due to the second species, the value being 
+ 2.2 x 10%. The smallest is due to the first 
species of transport, i.e., + 1.2 x 10%. Thus 
about 14 per cent of the total is brought about 
by mean meridional circulations, according to 
the data. 

(3) If use is made of parentheses to indicate 
vertical averaging with respect to standard 


atmosphere pressure, the quantity ([v]) has a 
value of — 0.16 m sec”! as already stated. 


Similarly ([u]) turns out to be + 0.50 m sec”. 


Since one may write that [4] = ([u]) + [ul 
with a similar relation involving v in place of 
u, it follows that 


(Ed fol) = (fal) (el) + (el EN) 
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The term on the left is a measure of the 
transport of the first species integrated ver- 
tically. The first term on the right represents a 


portion of this quantity depending on ({v]) 
which as mentioned before may be spurious. 
Although this “correction” procedure is not 
recommended as a substitute for better obser- 
vational material, it is of interest to note that 
in the present case the correction is of negligible 
magnitude compared to the last term which 
contains only departures from the vertical 
averages. 


(4) In the process of calculation, individual 


monthly averages of [u’ v’] were obtained for 
the twelve levels used. Of these 144 monthly 
values 98 were positive in sign, whereas at 
31° N 138 were positive. There appears to be 
a tendency for this quantity to be negative at 
high levels as is also shown in table 3. However, 
caution must be used in accepting this feature in 
view of the large confidence limits associated 
with these data. 


(s) The correlation coefficients r for the N 
pairs of # and v given in column 14 of table 
3 are smaller than the corresponding ones for 
31° N except at 2,000 feet, and actually become 
slightly negative at the highest three levels. 


5. Some meteorological implications 


(a) From time to time various estimates have 
been made of the strength of mean meridional 
circulations which would be needed in order 
to transport the net amount of angular momen- 
tum across given latitude circles prescribed by 
mean surface torques, on the assumption that 
the entire flux is of the first species. Two 
difficulties attend this process. In the first place 
accurate estimates of the requirements due to 
surface torques are difficult to make. In the 
second place there is no objective method for 


specifying the functional dependence of [v] 
on: elevation (pressure). It is of interest to 
substitute the total flux as measured by studies 
of the present kind for the surface torque 
requirements, and to employ a variational 


method to specify [v] as a function of pressure. 
In effect one thus seeks the least intense merid- 
ional circulation which could bring about the 
stated transport T. This procedure was sug- 
gested to the writers by Messrs. E. N. Lorenz 
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Fig. 1. Mean meridional circulations in m/sec at 13° N 

and at 31° N as given by data (full curves). Weakest 

mean meridional circulations capable of transporting 

total amount of angular momentum given by dashed 
curves. 


and H. L. Kuo of the Massachusetts Institute 
of Technology. 


The process takes the form of finding the 
solution of the variational problem 


6S Tt dp = 0 


subject to the side conditions that 


a 


J [v] dp = 0; 


b 


[u] [v] dd = KT 


oS 
a 


Here standard atmosphere pressuresar e used, 
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a and b being the values at 2,000 and 55,000 


feet respectively. The dependence of [u] on p 
is that given by the data and K is a constant. 
The result is simple to state qualitatively. The 
“least square” cell so obtained has the same 


shape as the vertical profile of [u], but displaced 
relative to the origin in order to satisfy the 
first side condition which stipulates that there 
shall be no net transport of air. The amplitude 
is of such a magnitude as to satisfy the second 
side condition, i.e., to provide the specified 
momentum transport. 


The profiles of [v] so obtained are given 
by the dashed curves in Fig. 1, together with 
the actually observed profiles (full curves), at 
13° and 31° N. The contrast between the two 
sets of curves is obvious and probably cannot 
be ascribed to poor cbservational material. 


It would of course be possible to impose 
additional side conditions upon the problem 
stipulating, for example, that the mean merid- 
ional circulation should transport stated 
amounts of energy of several forms. How- 
ever, such added restrictions would not lead 
to less vigorous circulations in the least square 
sense. 


(b) According to what has been said under 
(1) of the previous section about 23 x 10% 
gm cm? sec! represents the average rate at 
which angular momentum is supplied to the 
zone between 13° and 31° N by surface 
torques. Recent estimates of such surface 
torques made by PRIESTLEY (1951b) from 
surface stresses yield a value of about 15 x 10% 
gm. cin? sect 


The average zonal component of the surface 
stress for the belt as derived from the present 
figures is approximately 0.53 dynes cm-?, 
which is somewhat larger than the estimate 
used by Priestley. 
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A Comparison between Observed and Computed Winds 
with Respect to their Applicability for 


Vorticity Computations’ 


By ERNEST HOVMOLLER, University of Stockholm 
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Abstract 


The insufficiency of computed height values for computations of derivatives such as vorticity 
and convergence, at least for the 300-mb and higher levels, is illustrated by theoretical and 
practical examples, and it is recommended to make more extensive use of observed winds. 
This gives a possibility to take into account the cyclostrophic and other non-geostrophic wind 
components for the vorticity computation which is the first step in the numerical forecasting 


methods (CHARNEY 1951). 


I. Introduction 


The approach towards numerical forecast- 
ing methods during the last few years has 
emphasized the need for ample and exact 
observations from the upper air, for a detailed 
and skilful map analysis, and for empirical 
support of the theoretical assumptions on 
which the methods are based. 

One of the assumptions usually made is that 
the absolute (or, in three-dimensional models, 
the potential) vorticity of an air parcel may be 
considered as a constant, or at least as a pro- 
perty which is sufficiently conservative to 
allow a forecast of 24 hours ahead or more. 
Even if this assumption holds true, there are 
several reasons why we cannot expect the 
resulting prognostic charts — c.g., maps 
showing the contours of the soo-mb surface 
— to be absolutely correct: 

(a) The aerological observations are always 
too few to permit a really detailed analysis 
of the whole map, and the errors in the obser- 


1 The inyestigations reported in this paper have been 
supported by a grant from the U.S. Weather Bureau. 


vations cannot all be corrected with a sufficient 
degree of certainty. 

(b) The use of grid points means that the 
vorticity is computed only at a relatively small 
number of points, and for these points a rather 
crude method of computing the relative vor- 
ticity (by means of finite differences) is used. 

(c) Additional theoretical assumptions con- 
cerning the transport of absolute vorticity 
are necessary to permit the construction of a 
soo-mb prognostic chart, namely: The vor- 
ticity is transported with the horizontal com- 
ponent of the geostrophic wind; the vertical 
transport of vorticity may be neglected; the 
horizontal divergence is assumed to be zero, 
which means that the effect of a three-dimen- 
sional field of deformation on the horizontal 
distribution of vorticity is neglected. 

The complex character of the whole prob- 
lem makes it difficult to separate the effect 
of any individual source of error. Thus, for 
instance, it is hardly possible to draw con- 
clusions from the failure of one or a few 
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forecasts as to the validity or non-validity of 
the basic theoretical assumptions. It has 
been tried, however, (NEWTON et AL., 1951) 
to investigate the variation of the abso- 
lute vorticity with time in developing kine- 
matic systems, and it has been indicated that a 
mechanism producing absolute vorticity may 
be of importance in certain situations. 

The present study deals with the problem 
of computing the vorticity from actual ob- 
servations. It arose from an attempt to con- 
struct vorticity charts for a selected period 
over a limited area (the British Isles) where a 
dense network of aerological stations provide 
for a material which is, qualitatively as well 
as quantitatively, probably superior to any 
other material now available for an area of 
equal size. 


2. Computation of the geostrophic wind 
field by means of temperature soundings only 


If the topography of a constant-pressure 
surface is completely known, the geostrophic 
wind can be derived for any point of this 
surface by determination of the slope, taking 
into account the Coriolis parameter as a func- 
tion of latitude. By means of observations 
made at three adjacent points A, B, and C 
(not in one line), we may determine the mean 
slope and, approximately, the mean geo- 
strophic wind within the triangle ABC. By 
means of observations at A and B only, we 
may determine a mean value of the geo- 
strophic wind component perpendicular to 
the direction of AB. Occasionally such cal- 
culations are made, mainly in connection 
with the construction of vertical sections. 

Let us assume that the height of the 1000-mb 
surface at stations A and B is known with 
sufficient accuracy. The actual (harmonic) 
mean of the absolute temperature of the layer 
between 1 000 and p mb may be T4 and Tz 
respectively, but, due to errors in the sound- 
ings, we have obtained instead T4’=T4+AT4 
and Tr’ = Tg-+ ATs. This gives us for the 
height of the isobaric surface p above 1000 mb 


AT. 
eo di ZA (: ae 7), 
oh Sui (a 


Z = ZB + AZB = Zp (+7 
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where Z,’ and Zp’ are computed, Z4 and Zz 
actual values. 


For the difference Z4’ — Zp’, we get 


EEG (1 + oe) + 


A 


(2) 


Sy ps (: + +), 


Ts 


which deviates from the actual difference Z4 — 


ZB by 


A (Za Zp) = AZ 4 AZ = 
er Tg ATs () 


EZ ih 
and as the quotients 2% and — are not very 
ZB TB 


different from 1, this may be approximated by 


7) 
A(Z4— Zp) © = (A T,—A Ts) — 
A 
Za 
= A(T4 — Th) 7, (2) 
the relative error is 
A (ZA = Zp) I (A T4— A T3) 24 es 
ZA LE (Za— Zp) Ta 
_ A(Ta—Ts)Za 6) 
(Z.4— Zp) Ta ~ 


If the distance between A and B is d, we 
have for the geostrophic. wind component 
(v,)„ perpendicular to AB 


CS ae (6) 


by using Z4’ and Zp’ we get, however, 
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Fig. 1. Error in height difference Z 4 — Zp between 
two radiosonde stations, À and B, caused by an error 
of 1° C in mean temperature difference T4 — TB. 
The height error is given as a function of the height 
ZA(x Zp) and the actual height difference Z 4 — Zp, 
and expressed as a percentage of Z4 — Zp. The figure 
also gives the error in the computed geostrophic wind 
component normal to AB in per cent of the corre- 
sponding component of the actual geostrophic wind. 


and the difference 
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Fig. I gives an evaluation of this expres- 
sion under the following assumptions: 

A(T4—Ts) is taken to be 1° C. As the error 
in v, is proportional to the error in tempera- 
ture difference, the diagram may of course 
be used for any value of A(T4—T;) if the 
percentage found from the diagram is multi- 

For T4 the value 260° has been chosen. The 
deviation from this value is in most practical 
cases less than 10%, hence the values in the 
diagram are reasonably correct for practically 
any temperature. 

The figure shows, to take a realistic example, 
that at the 300-mb level (z4 approximately 
9,000 m) the error amounts to 25 % of the 
comparable component of the actual wind, if 
the true height difference between the 300-mb 
level at the two stations is 140 m and the 
error in temperature difference is 1°. 

Evidently, as far as the geostrophic wind 
field of the upper troposphere and the strato- 
sphere is concerned, it is not advisable to 
draw any conclusions from height differences 
based on unsmoothed height values. On the 
other hand, any kind of smoothing will tend 
to efface details which may be both real and 
significant. As an incorporation of observed 
winds would introduce non-geostrophic wind 
components of which our knowledge is very 
deficient, it must be stated that no general 
method exists by which it is possible to deter- 
mine the detailed distribution of geostrophic wind 
in the upper troposphere and the stratosphere. 
Consequently, it is not possible to obtain 
(from daily maps) more than a very general 
picture of the geographical distribution of 
non-geostrophic wind components in these 
levels. Even a considerable increase in the 
number of observations would not bring about 
any material change in this state of affairs, 
whereas a further improvement of the quality 
of the soundings would be more important 
from this particular point of view. 

It may be noted, that besides temperature 
errors, pressure errors too occur in the sound- 
ings. In the troposphere these are almost 
equivalent to the temperature errors, the ‘‘scale- 
value” being dependent upon the actual lapse- 
rate but, at first approximation, independent 
of the height. In the stratosphere, the relative 
importance of pressure errors for computations 
of height (and geostrophic winds) decreases 
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with altitude, as the absolute height error 
remains virtually constant above the tropo- 
pause level. The existence of pressure errors 
thus tends to modify, for levels above 300 mb, 
the normal increase of errors with height 
shown in Fig. 1. On the other hand, tem- 
perature errors often increase, and seldom 
decrease, with height, which accentuates the 
tendency shown in the figure. 

In order to test the accuracy of geostrophic 
winds which may actually be obtained by 
use of uncorrected height values, a random 
example (May 29, 1951, 152) was chosen 
among a number of cases in which complete 
observations of temperature and wind were 
available from all existing British radiosonde 
stations up to 100 mb. A triangular net was 
introduced (Fig. 2), and the geostrophic wind 
was computed for each triangle, using un- 
corrected height values for 500, 300, 200, 
and 100 mb respectively. In the polar diagrams 
shown in Fig. 3a-d each of the computed 
winds is represented by a cross, and each 
actual wind observation by a dot. The figures 
show, for 200 and above all for 100 mb, in a 
striking manner the erratic distribution of 
the crosses and the denser ‘‘clouds’ formed 
by the dots. It should be stressed that the 
scattering of the dots is not necessarily due to 
lack of exactness or representativity of the 
wind data, as there is no reason to assume 
that the actual wind, or the representative 
wind, i.e., the actual wind corrected for the 
(unknown) effect of turbulence is exactly the 
same at all stations. The analagous statement 
might be made as far as the computed winds 
are concerned, though the real differences in 
this case might be supposed to be smaller 
because the centers of gravity of the triangles 
are distributed over a smaller area than the 
stations. In fact, some of the ‘‘absurdly” 
placed crosses (c and g on Figs. 3 c—d) be- 
long to the only triangles without “‘outer 
walls”. The real differences in geostrophic 
winds are, in the stratosphere at least, one 
order of magnitude less than the fictitious diffe- 
rences which result from a computation based 
on height differences. This result, which has 
been derived from a situation with relatively 
quiet conditions in the stratosphere, is not 
directly applicable to situations where a strong 
horizontal shear is present, but it is seen that 
the actual shear must be extremely large to 
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Fig. 2. Triangular net used in comparing computed 


geostrophic winds and observed winds. 


match the “false shear” normally introduced 
by temperature errors. 

For levels above 300 mb the observed wind (if 
obtained by radar or an approximately equally 
good method) normally is a better approximation 
to the geostrophic wind than a value derived from 
computed heights, in spite of the systematic 
difference existing between geostrophic and 
actual wind. For the 300 and soo mb levels, 
the same statement probably is true in many 
cases, but not generally. 


3. Computations of vorticity by means of 
geostrophic winds based on temperature 
soundings only 


If height values computed from temperature 
soundings are not sufficiently exact to allow one 
differentiation, giving the geostrophic wind, 
they must, evidently, be quite unsuitable for 
computations based on a second differentiation, 
such as computations of horizontal shear or 
vorticity. (In the case of horizontal conver- 
gence, the value obtained by means of com- 
puted heights necessarily must be equal to 
zero, which is obviously wrong. Nevertheless, 
the absolute error in the case of this compu- 
tation may be considerably smaller than in the 
case of vorticity computations.) 
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Figs. 3a 


d. Polar diagrams showing the computed geostrophic winds and the observed winds in a case (May 


29, 1951, 152) where complete temperature and wind observations were available for 9 British stations as far 
up as to 100 mb. Each cross denotes a geostrophic wind computed for one of the triangles in Fig. 2, each dot 
an observed wind, the letters and numbers referring to Fig. 2. 


To give an example showing the magnitude 
of the errors, the “fictitious vorticity’’ emer- 
ging from one erroneous temperature sound- 
ing surrounded by correct soundings shall be 
discussed. 

Suppose that the computed height at the 
station Ais Z’=Z+AZ, where Z is the actual 


height and AZ the error. The height values 
at the distance r are supposed to be correct. 
Whether the vorticity in the point A is com- 
puted by comparing Z’ with the average 
height along the circle or at a few fixed points 
at the distance r, the error AZ will introduce a 
fictitious vorticity AË which is added to the 
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actual vorticity. To compute A¢ we have to 
determine the “false geostrophic wind” along 
the circle, which is 


Ayes ee. (12) 


this gives for the fictitious angular velocity 


dar + => AS, (13) 


jane 


and for the fictitious vorticity 


(14) 


Here we may insert the value of AZ which 
can be obtained from (4) by putting A (Z4—Z5) 
=AZ and, consequently, A (T4— Tr) =AT 4: 
AT: 


IB ee 
AFS T, a 


ZA 


(15) 


As AË is proportional to Z,, it is normally 
much greater in the stratosphere than in the 
middle troposphere, even when the tempera- 
ture error does not increase with height. The 
variation with T 4 is unimportant; the variation 
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Fig. 4. Fictitious vorticity resulting from an error of 
1° in the mean temperature at A in the model de- 
scribed in text to Fig. 4. The vorticity is expressed as 
a function of the radius r, the height Z4, and the 
latitude. For instance, if the height is 10 km, the lat- 
itude 45°, and the radius (repesenting the distance be- 
tween adjacent radiosonde stations) 300 km, the ficti- 
tious vorticity is 8:1075 sec" per degree Centigrade. 
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for computations of the vorticity field. 


Diagram illustrating the polygon method 


with latitude is, of course, the same as for the 
geostrophic wind itself, which means that 
the error is particularly large at low latitudes. 

Fig. 4 shows the numerical value of AC 
at 30—60° latitude as a function of Z4 and 
HARAS 1°; 


4. Computation of vorticity by means of 
observed winds only 


On a few occasions, observed winds have 
been used for a computation of horizontal 
convergence (Byers and RODEBUSH, 1948, 
Byers and HULL, 1949), or vorticity (BELLAMY, 
1949). This may be done by combining the 
observed winds in groups of three or four 
and determining the radial, respectively tangen- 
tial, component of the wind relative to the 
center of gravity of the polygon (Fig. 5). 
The mean vorticity, e. g., for the quadrangle 
A,4o4,d, is then, approximately, 


2 Vt 
= 25 == N 
n t)» 


v, being the tangential component (positive, 
if the motion is counterclockwise), r the 
radial distance, and p the index number of the 
point (1 S p Sn). 

This method, of course, is rather crude, as 
each computation is based on a very limited 
number of observations without any warrant 
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Fig. 6. Diagram illustrating an alternative method of 
computing vorticity values, using an analysed field of 
wind components and a grid of semi-squares. 


that these characterize sufficiently well the 
wind distribution within the polygon. Little, 
if anything, could be gained by giving the 


: UES ae : 
different terms — different weights (to com- 
r 


pensate for the non-regular form of the poly- 
gon). A minor difficulty arises from the finite 
and unequal size of the polygons: each vor- 
ticity value computed in the manner just 
described is an approximate mean value for 
the polygon; in the final analysis, the field 
distribution of vorticity therefore should be 
such that the mean value for each polygon 
should be equal to that found by the compu- 
tation. This would be of importance mainly 
in the vicinity of marked vorticity extremes. 

The polygon method is inconvenient if the 
total number of stations reporting upper 
winds is considerable. This is true partly 
because of the arbitrariness of the structure 
of the polygonal net, but mainly because of 
the large work involved in measuring angles 
and distances and making trigonometrical com- 
putations. An arbitrary method which is, 
from this point of view, slightly more con- 
venient, is the following: 

As the relative vorticity is defined by the 
formula 


(16) 


it can be derived by means of maps showing 
the horizontal and meridional wind compo- 
nent separately. Such maps can easily be con- 
structed, if a reasonably large number of well- 
distributed observations is available within 
the area to be investigated. Probably it is 
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justified, at least as a general rule, to use only 
plain interpolation principles for the construc- 
tion of u- and v-isotachs, but nothing will 
prevent the use of more refined methods of 
interpolation in special cases. 

Although the picture which the wind com- 
ponent maps give of the wind distribution is 
based on rather few observations and a not too 
well-founded hypothesis regarding the 
“smoothness” of the wind field, it is thought 
that vorticity computations based on such 
maps are at least as reliable as computations 
made by any other method. As to the question 
whether the wind observations should be 
accepted, without any smoothing, in drawing 
the curves, no general answer can be given. 
If the quality of the observations is very 
good, it is probable that nothing is gained 
by a smoothing unless the network of stations 
is unusually dense. 

The procedure by which the vorticity it- 
self is computed is shown by Fig. 6, repre- 
senting a section of a map for an arbitrary 
pressure level. The solid lines show the west 
component, the broken lines the south com- 
ponent as obtained by the interpolation 
described above. The mean value of vorticity 
within the square shown on the figure can 
now be approximated by 


¢ =k{[(vs—vp)—(uc—ua)], (17) 


the value of the proportionality factor k 
depending upon the size of the square (but 
not on the latitude). 

In practical analysis, it is convenient to use a 
network of squares of the type shown in 
Fig. 6. Because of the curvature of the earth, 
this involves making a compromise with 
respect to the form of each square. If the 
dimensions of the area are of the order of 
magnitude of 10° or 20° of latitude and longi- 
tude, a reasonable compromise may be seen 
from Fig. 7. At first, a meridian AB dividing 
the area in nearly equal parts is chosen as a 
base-line. Along this base-line, division marks 
are made at equal intervals, a convenient inter- 
val being 1° of latitude or 100 km. Next, the 
latitude circles through the division marks 
are drawn, and along each of these circles 
new division marks are made starting from 
the original base line and using the same 
length unit as a base. Of the two systems of 
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main diagonal lines which may be constructed 
in this network, every second diagonal is 
actually drawn; the points which are not on 
any of these diagonals are used as grid points 
for which the vorticity is computed in the 
manner described above (see Fig. 6). The 
fact that the figures are not ideal squares is of 
little or no importance. 

The method outlined above has been used 
in an investigation of the life history of areas 
of excessive vorticity over the British Isles. 
The result of this investigation was meager, 
mainly because the area having a sufficiently 
dense network was not sufficiently extensive, 
but it is thought that the method of computing 
the vorticity proved to be advantageous, or 
at least practicable. 


5. Computation of vorticity by use of maps 
based on all available material 


In most cases the vorticity is actually com- 
puted from contour maps of the ordinary 
type. These maps are constructed by means 
of both computed heights and observed winds, 
and it is generally accepted as a rule for the 
analysis that at least for the 300-mb and higher 
levels the observed winds give a better idea 
of the gradient than the height differences. 
However, if we use the contour values to 
compute, by differerentiating twice, the vorti- 
city, we use the geostrophic wind components 
only, as this is what can be derived from the 
maps. It is not possible, generally, to construct 
a chart of the ‘‘potential of the wind”, from 
which the actual wind might be taken in a 
manner analagous to that which is used by 
determining geostrophic winds from a con- 
tour map. Even the simple transition from 
geostrophic to gradient wind cannot be achie- 
ved by replacing the contour lines by another, 
similar set of curves. 

Thus, the possible (or, rather, probable) im- 
portance of non-geostrophic wind components 
for the distribution of vorticity is left out of 
consideration in any case where the vorticity 
is computed by differentiating twice the 
contour values. As the additional vorticity 
which depends on non-geostrophic winds is 
likely to mix up with the geostrophic vorticity, 
we may see in this limitation of the standard 
methods one reason, though possibly not 
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Fig. 7. Model of a network which is convenient for 
computation of vorticity values from maps showing 
distribution of zonal and meridional wind components. 


a very important one, why the result of 
numerical forecasting cannot be absolutely 
correct. 

It might, however, be worth while to try 
numerical computation based on a vorticity 
map which is constructed by using observed 
winds wherever available, and geostrophic 
winds only where the wind data are insufli- 
cient. The procedure might be this: from an 
ordinary contour map, approximate compo- 
nents of the geostrophic wind are determined 
and plotted for a sufficiently large number 
of points (not necessarily the same points 
for each component) — the zonal compo- 
nents by means of height differences along the 
meridians, the meridional components by 
means of height differences along the latitude 
circles. It will be advantageous to plot the 
numbers on two sheets of transparent paper, 
one for each component. On the same sheets 
the components of the actual winds are plotted, 
most conveniently in another color, and by 
the analysis more stress is laid upon these 
values than on the geostrophic components 
(which are necessary, however, to get a rea- 
sonable picture of the wind distribution over 
areas where wind observations are missing or 
too scanty). The two components of relative 
vorticity are then readily computed, and by 
numerical or graphical addition of these com- 
ponents plus the Coriolis parameter the vor- 
ticity field is determined. 
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Some Remarks on the Angular Momentum Balance 


in the Atmosphere 


By J. van MIEGHEM, University of Brussels! 


(Manuscript received 22 April, 1952) 


1. By substitution of the equation of con- 
tinuity in the equation of the zonal motion, we 
obtain (vAN MIEGHEM 1950) 


E 


) ; 
; (ou) - — [r? cos p (uov' + 


r cos @ dx’ 
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where (x! = A, x? = q, x? =r) are the spher- 


> : dx’ 
ical coordinates, v! = SE the contravariant com- 
[A 
= F dÀ 
ponents of the air velocity v, (4 = r cos p Tr 
= t 
d dr pee 
Dir + = 7 the zonal, meridional and 
t E 


vertical components of v), u = r cosq(r cos: 
w+ u) the absolute angular momentum, & 
the earth’s angular speed, o the specific mass, p 
the pressure, Ti the contravariant stresses along 
the latitude circles, due to internal friction, 
ef =1 when i =1 and ¢,/ =o when i +1. 

Equation (1) shows that there is no pro- 
duction nor destruction of absolute angular mo- 
mentum in the atmosphere (VAN MIEGHEM, 1950), 
but only redistribution of momentum as a con- 
sequence of inflow and outflow of west-momen- 
tum in the surface easterlies and in the surface 
westerlies respectively (skin-friction and moun- 
tain-effect at the earth’s surface). 

2. Averaging equation (1) along the latitude 


circles, one obtains 


1 At present at the Institute of Meteorology of the 
University of Stockholm. 
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a on eo aS 2 % 
dt (om r2 cos œ 0) oz a Az a 
+ r cos p guy" — Ti] = 0, (x = 2, 3) 


(2) 
where X designates the zonal mean value of X 
and X the corresponding weighted mean value 
defined by ox = = 0X. Let X’ and X” be the 
fluctuations of X with respect to X and X. It 
is obvious that X’ =o, HN = = 0; hence 


NE STE 0 0 
= 9X = @X (3) 
and 
DOSE (4) 


For the zonally averaged instantaneous absolute 
angular momentum fi, three modes of redistribu- 
tion are possible: 


a) redistribution of momentum due to mass 
transport ov" in meridional planes; 

ou"w"); 

c) a non convective flux due to the mean 
zonal components of the viscosity and 


small-scale eddy stresses. This flux may be 
disregarded outside the boundary layer. 


b) an eddy flux (ou”’v”, 


3. Replacing in (3) X by v, one finds 


ev=ov 


(3a) 


ie LA 7 ed: 
=ov ov 


and identical relations for 1. 
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The vector v” (v”, w”) may be considered 
as defining the organization of the mass ex- 
change in meridional planes as a consequence of 
the large-scale eddying motion (vAN MIEGHEM 
1949). When, along the latitude circles, the 
fluctuations 0’ and v”* (or v’*) are always of 
the same or the opposite sign, that is to say 
when the mass exchange in meridional planes 
due to the large-scale eddies is well organized, 


the velocity components y”* reach a high abso- 
lute value. 

For instance, in most cases, warm air (o’ < 0) 
ascends (w > o) and comes from the south 
(v > 0), while cold air (0° > 0) subsides (w < 0) 
and comes from the north (v < 0). Hence, along 
the latitude circles, w” > o and, below 250 mb, 
v" >> o. However, along the latitude circles above 
the 250 mb-level, v” < 0. Fromg@ X=0X+0'X’ 
and (3), it follows that: 


ov =o —) and 


ow = p(w— w) with w" > o. (s) 


The instantaneous mean meridional circulation 
is generally defined by (v, w) although it would 
have been better to define this mean circulation 
by (?, w). In order to avoid confusion, we shall 
use the definition (v, w) introduced previously. 

Now, along a meridian, for continuity reasons, 
we must have 

TI 
t+r en 
TETE ow dp = 0, 
t 7 


2 


if the time interval 7 is sufficiently large. Whence, 
by virtue of (s), 

t+r dur t4r 

[4 f owdp>o and | owdt +. 

t a t 
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This inequality demonstrates the existence of a 
mean meridional circulation. Moreover, the hor- 
izontal gradient of 9 being small, the lifting 
in equatorial regions (# > 0) must overcom- 
pensate the sinking in higher latitudes (w < o). 
For the same reasons, along a vertical, 
t+ co 


foie =o 
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On account of the exponential decrease of p 
with altitude z, the sum of the positive ov’, 
above about 10 km, presumably does not com- 
pensate the sum of the negative gv’, below 
10 km, so that 


t+r oo i+r 
f dt fovdz>o and f ovdt +0. 
t fo) t 


Hence, the upper layer with v > o must be 
much deeper than the lowest layer with v < 0, 
as a consequence of the rapid decrease of 9 
with height. 

It should be emphazised that this existence 
proof of a mean meridional circulation defined 
by the time mean value of pv and pw does 
not give any information about the order of 
magnitude of the instantaneous mean meridional 
velocities (v, w), except that they are very small. 

From (5) it follows that, as a result of the 
large scale eddy exchange of mass in meridional 
planes: a) the vertical branch of the instantane- 
ous mean meridional circulation (v, w) is weak- 


ened in equatorial latitudes (w > 0, w” > 0) 
and strengthened in higher latitudes (w < 0, 


w” > 0); b) the horizontal branch is strength- 


ened in the low troposphere (v < 0, v” > o) 


and in the stratosphere (v > 0, v” < o), but 
weakened in the middle and high troposphere 


(> 0, v"> o). 
4. The flux of @-momentum assumes the form: 


| r? co poor = r? cos? po (r — v"), 6) 
| r cos? po o W = r? cos? mao (w— w”). 


In virtue of (3), one obtains for the eddy flux 
of angular momentum 


r cos you" v" = rcos@p + 


(Gr a) + owe), 
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: [5 (u w’ — u"w") + o ww]. 
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Hence, the eddy flux of angular momentum 
may be decomposed into: 


a) a flux (ow’v’, ou’w’) due to the correla- 
tion between the velocity fluctuations only. It 


seems now well established that the flux ou'v’ 


ANGULAR MOMENTUM BALANCE IN THE ATMOSPHERE 


represents the bulk of the horizontal poleward 
transport of momentum, repuired in order to 
compensate the outflow in the surface westerlies 
and to maintain the strong upper westerly flow 
(Mintz 1951, STARR and WHITE 1951, 1952). The 


flux Qu'w' is presumably downward (Wire 
1950), so that it counterbalances the upward flux 
of angular momentum due to the mean meridi- 
onal circulation (BJERKNES 1951). 

b) a flux (— ou" Dv", — ou" w") due to the 
organization of the large-scale eddy exchange of 
mass in meridional planes. In the case of a me- 
ridional tilt to the east of throughs and ridges, 


== ou” Zu ou < 0; consequently this 
flux is directed to the south and to the earth’s 
surface. 

c) a flux (o’u’v’, p’u’w') due to the cor- 
relations between both the velocity and density 
fluctuations. 


The fluxes b and c are, as a rule, disregarded 


and therefore the flux a is generally called the. 


eddy flux. 


On the other hand, the horizontal and ver- 
tical components of the convective flux of rel- 
ative angular momentum may be expressed in 
the following way; 


[ rcos POUV = rcos po 
| [7 + PES Er va 
ae 2 (8) 
rcospouw =rcosgo: 
| : Kt + u" w" — u" w — uw”), 


account being taken of the general formula 
ox = oX +o X' — ox lat (a) This 
convective flux has four constituent fluxes: 

a) a flux (ouv, ouw) due to the instantane- 
ous mean circulation (y, wv) in the meridional 


planes. The time mean value of #v for a large 
period seems to be very small (Starr and WHITE, 
1952) ; 

b) a flux (ou v”, ou" w") due to the orga- 
nization of the large-scale eddy exchange of 
mass in meridional planes. It is remarkable that 
this flux drops out in the sum of the fluxes 
(7) and (8); 2 

(c) a flux (—ou"v, — ou" w) proportional to 


the mean meridional mass transport (0v, 0) 


$*— 202040 
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and the correlation between the fluctuations of 
the density of the zonal wind; 

d) a flux (— our”, — Guw") proportional to 
the mean zonal wind and the “organization” 


vector (v”, w”). 
This last flux is presumably more important 
than the former. 


5. Finally, disregarding viscosity and small- 
scale turbulence, the instantaneous total flux of 
absolute angular momentum is defined by the 
horizontal (positive to the north) and the ver- 
tical (positive in: zenith direction) components 
given below: 


er = wr cos* pov 
+ rcos (on? + ou"v") = 
= rceospo(wrcosmv + uv + uv) =. 
— r cos g X 6 [or cos pv" be phe Ae 
+ who] — wor} (9) 
and 
TOW = or co pam + 
+ rcosp (On + ou"w") = 
= r cos po(wrcos pw + uw Ne 
— rcos g { 5 [or cos pw” + uw” + 


= u'w| — o'u'w' N (10) 


with the conditions 


= | v” > 0, below ~ 250 mb, 
w d 


= Cc) — 
| v" <0, above ~ 250 mb, and 


u” So, depending upon the sign of the me- 
ridional tilt of troughs and ridges. 

Hitherto, only the fluxes in  parenthesis 
in the last member of (9) and (10) have 
been considered and the time mean values of 
v, uv and w’v’ computed (STARR and WHITE, 
1951, 1952). We believe that the instantaneous 
zonal mean value v” may not be disregarded 
a priori and therefore should be computed. 
Moreover, the triple correlation o’u’v’ should 
be investigated numerically in order to establish 
its order of magnitude. 
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The Variation of Gravity within the Earth 


By SVEND SAXOV, Danish Geodetic Institute, Copenhagen 


(Manuscript received 29 March 1952) 


Abstract 


The different methods to determine the values of gravity within the Earth are re- 


viewed. 
By using the equation 


fe" 3 


r 2 
ge So (®) 


r? 


and substituting the density values tabulated by Ramsay the variation of gravity has been 


computed (table ı and fig. 1). 


The g value of 982 gal at the surface of the Earth is 


attained without any correction to the density distribution. 


Ramsay (1949) has recently developed a 
theory in which the Earth’s core and mantle 
are assumed to be chemically identical. The 
discontinuity of seismic waves at a depth of 
2,900 km. is interpreted in the usual way as a 
jump in density. But the jump in density is 
attributed to phase transition under pressure 
and not to the appearance of a new material 
such as iron. Olivine has a ionic structure in 
the mantle; in the core it is in a dense metallic 
state. The density figures given by Ramsay 
are slightly different from figures given by 
BULLEN (1940 and 1942) as the density values 
vary in the mantle from 3.29 to 6.02 and in the 
core from 8.92 to 18.5. The corresponding 
figures in Bullen’s papers are 3.32 to 5.68 and 


9.43 to 17.20. No gravity values have been 
computed from the density variation given by 
Ramsay, as far as it is known, and it is therefore 
the purpose of this paper to compute new 
gravity values. 

The value of gravity at a distance r from the 
center of the Earth is given by 


4% nf 
2 Gs 0, rèdf (1) 
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where G is the constant of gravitation and eve 
the density at the point considered. 

BENHELD (1937) used this formula regarding 
the Earth as a non-rotating sphere, an assump- 
tion which was justified by the uncertainty 


THE VARIATION OF GRAVITY WITHIN THE EARTH 


of the density values. From r = 0 to 1,070 km. 
gr was evaluated by using a power series of 
approximation for @,. Beyond this distance g, 
was obtained by numerical integration in steps 
of 200 km or less. 

STEWART (1949) has investigated theoretically 
the conditions under which maximum and 
minimum values of gravity occur. By means of 
equation (1), and by successive differentiation 
finally he obtained the equation which defines 
the critical values of g 


d ee 
ht) 6 


where 9 is the average density within the sphere 
of radius r. It is further pointed out, that if the 
depths and densities of the various zones are 
known, it is possible to evaluate exactly the 
positions of all maxima and minima, as well 
as the resultant values of gravity. 

Another method was used by BULLEN (1940 
and 1942). The values of gravitational force at 
various points are given, also, using the rela- 
ton 


c= GM () 


where M is the mass within a sphere of radius r 
concentric with the Earth, and G the constant 
of gravitation. Bullen proceeded this way 
obtaining the expression on the right-hand 
side in equation (3) from the following 


dor _ ol ee 
R=-— cm 5 (a— 4 pt) W 


where 9, is the density at a distance r from the 
Earth’s center, « and f the velocities of P and S 
seismic waves there. 

We have, however, preferred not to use 
equation (1) or (3) but to proceed the following 
way; considering the Earth is a sphere and 
made of concentrical strata of equal density it 
is shown (OLczaK 1938) that the gravity value 
g has the following expression 


where o is assumed to be the constant density 
in the concentric layer with the radii r and ro, 
and G the constant of gravitation. Having g on 
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Fig. 1. The distribution of density with depth (the upper 
part) and the variation of gravity (the lower part). 


the outer face of any layer we can compute g 
everywhere in its interior down to its inner 
face. Since g is a continuous function of r, the 
value obtained can be used in the same way 
for calculation in the following layer. 

In table ı the depths and densities according 
to Ramsay are listed in the two first columns, 
and the following column contains the cor- 
responding gravity values obtained from equa- 
tion (s). Fig. 1 illustrates the distribution of 
density with depth (the upper part) and the 
variation of gravity (the lower part). 

The value of g rises to a maximum of 999 
gal at a depth of about 6—700 km. The func- 
tion then falls to a minimum of 986 gal at a 
depth of the order of 1 800 to 2 000 km. rising 
to maximum (the greatest value) of 1 017 gal 
at the boundary between the Earth’s mantle 
and the central core at a depth of 2 898 km. 


I40 
Table 1 


Depth | Depth | Density | Gravity 


Density|Gravity 


SINEINID ES AS OM I 


The gravity values in the central core fall 
continuously from the boundary to the mantle 
and to a depth of 5 125 km. where the den- 
sity value has a jump, and the gravity continues 


km |g/cm?| gal km g/cm? gal at another rate to the center where is has its 

o | 2.72 | 982 || 2600 5.78 999 zero-value. 

33 | 3-29 | 985 || 2800 39 A910 The property used by BırcH (1939) that 
100 3:39 987 808 ff 6.02 nad ‘chi + i 
Ne tac Me as 1017 g is always within ı per cent of the value 990 
300 3.59 993 3000 9.09 1000 gal down to a depth of about 2,500 km re= 
400 | 3.68 | 996 3200 9.43 965 mains true. 

42315370 + ee a While the gravity values given by Benfield 
ae ae ae 3 > 846 were adjusted by 4, % to bring g at the sur- 
700 | 4.27 | 999 4000 10.63 804 face of the Earth to the value of 982 gal and 
800 | 4.38 | 998 || 4200 | 10.89 762 Olczak had to multiply his results by the 
DOOM Dee Ce ae Iu14 Ter factor 1.00572 in order to obtain 982 gal at 
1000 4.60 | 996 4600 11.38 682 iE £ ns ] RE 
DE ol a 648 the surface, we are attaining the value 
1400 | 4.90 | 990 || 4982 11.8 626 982 gal without any correction to the density 
OO M RCE ECC ER { 12.0 \| 618 values given by Ramsay. 
1800 | 5.20 | 986 T7230) 
2000 5.36 986 5700 18.0 342 
2200 | 5.49 | 988 || 6171 18.5 ° Thank is due to Mr. K. Nygaard, M. Sc., 
2400) 15.091992 for assistance in the tabulation. 
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Report on a conference on the application of numerical methods 


in forecasting atmospheric flow patterns. 


During the week of 12 to 14 May, 1952 a con- 
ference was held at the Institute of Meteorology, 
University of Stockholm, on the various aspects of 
numerical methods for forecasting atmospheric flow 
patterns. The following papers were presented. 
(Some of the papers will be published in forthcom- 
ing issues of Tellus). 

May 12 

Dr. A. Eliassen (Norwegian Meteorological 
Institute): A comparison of different dynamic models 
designed for computational purposes. The barotropic 
model developed by Charney and collaborators at 
Princeton (Tellus, 2, 237—254) as well as various 
simplified models developed more recently in an 
attempt to incorporate the baroclinic character of 
the atmosphere (i.e. Fjortoft, Tellus, 2, 237—254; 
Phillips, J. Meteor., 8, 381394; Eady, Eliassen and 
Bolin) were discussed and compared with specific 
attention to basic assumptions made and the ease of 
intrepreting the models in terms of the atmosphere. 

Fil. lic. B. Bolin (University of Stockholm): 
A series of tests of the barotropic model. The several 
independent series of tests of this model were discus- 
sed and results presented. (Tellus, 3, 248—257, 
Tellus, 4, 21—30, Tellus, 2, 237—254). 

Civiling. G. Neovius (Swedish Computer Pro- 
ject): Demonstration of the Swedish relay-computor 
(Bark). Bark is a Swedish realy computer similar to 
Harvard University Mark I and has been in use 
since spring 1950. An explanation of its operation 
was followed by a demonstration. 

May 13 

Dr. J. Bellamy (Cook Electric Company, U.S.A.): 
Comments on the communication problem. The speed 


x Deputy director of the Institute of Meteorology; 
on leave from Norwegian Meteorological Institute. 


of the modern computing machine points out the 
urgent necessity for speeding up the transmission of 
weather data from the observer or instrument to 
the forecaster. Since most of the time involved is 
spent in coding or decoding for purposes of trans- 
mission, accessibility and utility, it was suggested 
that a method of representation be developed which 
would serve all purposes and eliminate the necessity 
for more than one code. A unitary decimal digitary 
system which can be read automatically or visually 
was suggested as a possible solution. 

Prof. R. Fjortoft (University of Copenhagen): 
Preliminary tests of a simplified graphical method for 
integrating the barotropic vorticity equation. By ad- 
vecting the vorticity field with a spacially-averaged 
wind field a 24-hour forecast may be obtained in 
one time step. Since this whole method is done 
graphically a forecast may be completed in approx- 
imately three hours by one person. Results look 
encouraging although no statistics are available yet. 

Fil. kand. G. Kjellberg (Swedish Computer 
Project): The Swedish Electronic Computer (BESK). 
This computer is of the Princeton type and is ex- 
pected to be completed during 1953. It has an 
internal memory capacity of 256 40-binary digit 
numbers. The capacity will be increased by magnetic 
drum storage. The speed and capacity of this 
machine will permit its use for meteorological 
purposes. 

Prof. E. T. Eady (Imperial College of Science, 
London): Methods of representation and techniques of 
integration of Poissons and Helmholtz’s equations. Dr. 
Eady indicated some of the main points to be con- 
sidered when planning the continued integration by 


2 Temporarily attached to the Institute of Mete- 
orology, University of Stockholm. 
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mechanical methods of the Poisson and Helmholtz 
equations appearing in problems of largescale 
atmospheric flow. For a single integration, three 
main methods have been used, namely Green’s 
function, harmonic (or spherical harmonic) analysis 
and mechanical “relaxation”; the last is probably 
both the quickest and the simplest to programme 
for large networks, being particularly suitable when 
spherical coordinates must be used. However, the 
use of finite difference approximations to differential 
coefficients leads to errors which may mount up 
and become serious when the integration is repeated 
many times. 

May 14 

Dr. A. Eliassen and Lt. W. Hubert? (Office of 
Navel Research); A three dimensional study of the 
vertical motion and the vorticity balance in a blocking 
case. The barotropic model was applied to a blocking 
situation over Eastern Atlantic and Western Europe. 
The computed soo mb tendencies showed certain 
significant discrepancies. Tendency computations 
were also performed on the maps of height averaged 
with respect to pressure, but there were still con- 
siderable errors. The incorporation of the mean 
thermal vorticity advection improved the tendencies 
somewhat. 

Interpreting the errors as a combined effect of 
horizontal divergence and vertical vorticity trans- 
port, the vertical motion necessary to explain the 
local vorticity changes was computed at 850, 700, 
500 and 300 mb. The result is in excellent qualitative 
agreement with the vertical motion inferred from 
observed weather phenomena. 

Fil. lic. B. Bolin: A study of the three dimensional 
field of vertical motion in a selected case and some remarks 
on baroclinic models of the atmosphere. The vertical 
velocity was computed by the same method as 
reported by Eliassen and Hubert. The importance of 
the term expressing a transformation of vorticity 
around the horizontal axis into vorticity around a 
vertical axis was demonstrated in a synoptic example 
with a pronounced frontal zone. Some of the basic 
parameters for incorporating the baroclinic structure 
of the atmosphere in numerical forecasting proce- 
dures were discussed. 

Dr. C. W. Newton? (University of Chicago): 
Comments on the analysis and data problems associated 
with numerical forecasting. Since present numerical 
prognosis methods involve multiple differentiation 
of the .sobar.c contour field, the analys.s used must 
be more accurate than those in ordinary forecasting 


ZuSess note, PHL Ale 
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practices. Smoothness and space- and time-con- 
tinuity of the vorticity field are essential. Subjective 
analysis methods must be used at present since 
objective analysis of any but the larger-scale features 
of the circulanon 1s very inaccurate due to the in- 
homogeneity and inaccuracy of radiosonde observa- 
tions. For improvements in either subjective or 
objective analysis, an increase in the density of the 
radio wind-finding network is the most urgent need. 

Dr. S. Reisch (Royal Institute of Technology, 
Stockholm): A proposed analog type of computer for 
use in numercal forecasting. The principles for a 
machine for solving partial differential equations by 
relaxation methods were explained. By using a 
“spider” with a series of adjustable heads the infor- 
mation required from the storage fields (magnetic 
drum or magnetic tape) during the process of relaxa- 
tion are available simultaneously, whereby all transfer 
problems are simplified considerably. The compu- 
ational units also become simpler than those at 
present in use in electronic computers. 

May 15 

Prof. J. Van Mieghem? (University of Brussels): 
Energy conversions in the atmosphere on the scale of the 
general circulation. The equations of the mechanical 
energy of the zonally averaged motion, of the 
kinetic energy of the large scale eddies and of the 
internal energy are deduced from the equations of 
motion and from the first law of thermodynamics. 
The fluxes and the rates of production of these energy 
forms are given. The two most important energy 
conversions on the scale of the general circulation 
are: 1) For the maintenance of the large scale dis- 
turbances, the conversion of potential energy into 
the kinetic energy of the large scale eddies as a 
result of the positive correlation between the fluctu- 
ations of temperature and vertical velcity ; 2) For the 
maintenance of the mean zonal motion, the con- 
version of kinetic energy of the largescale eddies 
into kinetic energy of the mean motion as a result 
of the signs of the lateral Jeffreys stress and of the 
meridional mean zonal wind shear. 


Met. K. Hinkelmann (Weather Bureau, Bad 
Kissingen): A report on preliminary tests of a three- 
dimensional model in Bad Kissingen, Germany. 

From the quasi-static and geostrophic approx- 
imations, a linear, second-order differential equa- 
tion in the height tendency is derived, using 
pressure as vertical coordinate. The boundary 
conditions at the top and the bottom are used 
in an integrated form to avoid non-centered 
differences when a finite grid is used. The com- 


NOTE 


putation of the coefficients of the equation from 
the synoptic data is discussed. 

The equation is solved by relaxation for a 
standard atmosphere with the inhomogencous 
term equal to unity in one gridpoint, and zero 
in all other points. Three Green’s functions for 
the tendency field computed in this way, and 
with the disturbance located at 900, soo and 
200 mb respectively, are demonstrated. 


Met. S. Smebye? (Norwegian Meteorological 
Institute): A numerical test of a simplified atmospheric 
model. It is possible to take the thermal structure of 
the atmosphere into account in a crude manner by 
expressing the vertical temperature distribution as 
a function of the deviation of the 1,000—soo mb 
thickness from some standard value. The function 
is chosen so that a warm troposphere gives a cold 
stratosphere and vice versa. The height of an isobaric 
surface and the thickness are the variables to be 
evaluated. Surface tendencies may then be com- 
puted. In the first case computed the calculated 
tendencies are an improvement over those calculated 
with the barotropic model. 

Fil. lic. ©. Lônnquist (Swedish Meteorological 
and Hydrological Institute): Present methods in 
common use for forecasting upper flow patterns. Particular 
emphasis was placed on methods in use for preparing 
the prognostic soo-mb chart. It was stressed that 
the results of present methods must be used as a 
standard in judging the value of numerical forecast- 
ing techniques. The average of the correlation coef- 
ficients between forecast and observed 24-hour 
changes of the soo mb surface for two years at 
Swedish Weather Bureau was 0.72. The correlation 
coefficients were highest when pure extrapolation 
methods were successful. 


Mr. Ph. Clapp? (U. S. Weather Bureau): Applica- 
tion of the barotropic model to medium and long range 
forecasting. The simple barotropic vorticity equation 
has been used to compute height tendencies on a 
few 5-day mean soo mb charts. The results are 
approximately as good as similar computations made 
on synoptic charts despite the fact the time intervals 
between observed charts are 5 times as great. The 
method may be applied to forecasts of mean monthly 
circulation patterns provided some account is taken 
of the effects of transformation of potential into 
kinetic energy, surface friction, and non-adiabatic 
heating. A first approximation to these factors can 
be obtained by setting them equal to their climato- 
logical or normal values. 


2 See note p. I4I. 
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Dr. F. A. Berson (Swedish Meteorological and 
Hydrological Institute): Barotropic tendency com- 
putations on motions of different time scales. The 
response of s-day mean soo mb flow to barotropic 
redistribution of vorticity is tested in a summer 
blocking situation. The average linear coefficient 
of correlation between observed contour height 
changes in three days and those obtained from the 
barotropic tendencies is as low as 0.49 but increases 
to 0.71 by using a constant correction equal to the 
over all “error” for the entire period. A large part 
of the “errors” could be ascribed to two main 
effects: 

1) standing perturbation possibly set up topograph- 
ically by the Rocky Mountains; 

2) horizontal transfer of vorticity by eddies of 
comparatively small size. The coefficient of eddy 
diffusivity is evaluated for a period of 14 days 
and turns out to be negative. 

May 16 

Met U. Czapski? (University of Hamburg): A 
comparative study of barotropically computed 500 mb 
tendencies for different grid positions and grid sizes. 
Considerable differences between vorticity fields 
were obtained if different grid points and grid sizes 
were used in barotropic computations and similarly 
for tendency fields. An improvement of the forecast 
was obtained if all points were included in the com- 
putations instead of every other as was done in the 
tests reported by Bolin. It is important that com- 
putations are done for every grid point in the grid 
used. The most feasible grid size needs careful 
consideration. 

Dr. L. Vuorela? (University of Helsinki): A 
synoptic study of deformation fields. The horizontal de- 
formation has been computed for a situation where 
the barotropic tendency computation failed as well 
as for several other situations. The maximum values 
obtained meant a horizontal dilatation of the order 
of three times in six hours. 

Fil. lic. G. Arnason (University of Stockholm): 
Presentation of a baroclinic three dimensional atmos- 
pheric model. The model used consists of a troposphere 
with constant lapse rate and an isothermal strato- 
sphere. Four basic parameters are used to describe 
the temperature, pressure, and wind fields; 1) Surface 
temperature, 2) surface pressure, 3) tropospheric 
lapse rate and 4) height of tropopause. The math- 
ematical problem of solving the dynamic equations 
using this model is reduced to solving two partial 
differential equations in two dimensions. The 
surface tendencies and average vertical velocities are 
obtained. 
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A large number of atmospheric models have been 
suggested for numerical forecasting. It seems that 
the general approach will turn out to be fruitful, 
but we can not judge about the detailed path to be 
followed until a series of complete forecasts have 
been made using these various models. A few of the 
unanswered questions which received particular 
attention at the conference are as follows. In the 
first place is a point representation of the variables 
most feasible or is a representation by functions 
better? This question is intimately connected with 
the best way of solving the various differential 
equations. It also brings up for discussion the 
question of how to take care of inaccuracies that 
are introduced by the finite difference method. The 
magnitude of the deformation field indicates that 
features which at one instant are caught by the grid 
are sometimes lost on a later occasion. Furthermore, 
how many layers or parameters are needed in the 
baroclinic models in order to account for the major 
effects of the baroclinicity of the atmosphere? Many 
of these questions must wait for solution until the 
high speed computers are available to perform the 
complete forecasts. 

Until these computers are available for routine 
computations of this kind it is believed that the 
simple graphical method presented by Fjortoft will 
turn out to be very useful. The method is very well 
suited for being immediately incorporated into the 
daily routine of forecast centers and tried out in 
combination with the present methods in use. 

In this connection it should be stressed, however, 
that these methods merely give the three dimensional 
field of motion as a function of time. In order to 
prepare an ordinary forecast this field has to be 
interpreted in terms of weather. This problem repre- 
sents a major part of the problem of weather fore- 
casting. Here the large number of statistical 
investigations carried out at forecast centers and 
elsewhere all over the world should be utilized. 
However, the vertical velocity field, that is ob- 
tained explicitly as a part of the computations 
using any of the baroclinic models proposed, will 
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simplify this process of translating the flow pat- 
tern into weather phenomena. 

It was repeatedly stressed that an accurate and 
careful analysis of the atmospheric flow patterns is 
of primary importance. At present it is not possible 
to accomplish this except over limited areas. For 
example the number of upper air observations over 
the Atlantic Ocean is hardly sufficient for the needs 
asked for here. 

Even if the problem of numerical forecasting of 
the upper flow patterns is solved successfully, the 
process of preparing input data for the machine 
must be made more efficient. Here the communi- 
cation problem as well as the problem of objective 
analysis perhaps represent the two most important 
problems. With regard to the objective analysis the 
discussion brought out the need for a method of 
obtaining a more accurate analysis in the vicinity 
of a single station by consulting the time record at 
the station and making use of some dynamic 
principles. 

Some of the results reported here are very tenta- 
tive and should perhaps in many cases be interpret- 
ed as indications of where future research should 
be concentrated. It is quite obvious that the most 
important problem at present is to get high speed 
computors available for all different tests that could 
answer several of the questions formulated above. 
However it seems as important that a closer col- 
laboration be established between the mathematicians 
familiar with the machines and the meteorologists 
so that the former understand the basic problems 
of meteorology and the latter can think in terms of 
the machines. It may be of some interest to mention 
that an attempt to initiate such a closer collaboration 
will be started this fall at the Institute of Mete- 
orology, University of Stockholm in cooperation 
with the Computer Project in Stockholm. A series 
of lectures will be given dealing with the basic 
theory of the computers, programming and coding 
problems. 


Bert Bolin, Harriet Newton. 
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